PostGIS  2.4.9dev-r@@SVN_REVISION@@

◆ lwgeom_segmentize_sphere()

LWGEOM* lwgeom_segmentize_sphere ( const LWGEOM lwg_in,
double  max_seg_length 
)

Create a new, densified geometry where no segment is longer than max_seg_length.

Derive a new geometry with vertices added to ensure no vertex is more than max_seg_length (in radians) from any other vertex.

Input geometry is not altered, output geometry must be freed by caller.

Parameters
lwg_in= input geometry
max_seg_length= maximum segment length in radians

Definition at line 1739 of file lwgeodetic.c.

References COLLECTIONTYPE, LWCOLLECTION::geoms, LINETYPE, lwcollection_add_lwgeom(), lwcollection_as_lwgeom(), lwcollection_construct_empty(), lwerror(), lwgeom_as_lwcollection(), lwgeom_as_lwline(), lwgeom_as_lwpoly(), lwgeom_clone(), lwgeom_clone_deep(), lwgeom_has_m(), lwgeom_has_z(), lwgeom_is_empty(), lwline_as_lwgeom(), lwline_construct(), lwpoly_add_ring(), lwpoly_as_lwgeom(), lwpoly_construct_empty(), lwtype_name(), MULTILINETYPE, MULTIPOINTTYPE, MULTIPOLYGONTYPE, LWCOLLECTION::ngeoms, LWPOLY::nrings, LWLINE::points, POINTTYPE, POLYGONTYPE, ptarray_segmentize_sphere(), LWPOLY::rings, LWGEOM::srid, and LWGEOM::type.

Referenced by geography_segmentize(), and test_lwgeom_segmentize_sphere().

1740 {
1741  POINTARRAY *pa_out;
1742  LWLINE *lwline;
1743  LWPOLY *lwpoly_in, *lwpoly_out;
1744  LWCOLLECTION *lwcol_in, *lwcol_out;
1745  int i;
1746 
1747  /* Reflect NULL */
1748  if ( ! lwg_in )
1749  return NULL;
1750 
1751  /* Clone empty */
1752  if ( lwgeom_is_empty(lwg_in) )
1753  return lwgeom_clone(lwg_in);
1754 
1755  switch (lwg_in->type)
1756  {
1757  case MULTIPOINTTYPE:
1758  case POINTTYPE:
1759  return lwgeom_clone_deep(lwg_in);
1760  break;
1761  case LINETYPE:
1762  lwline = lwgeom_as_lwline(lwg_in);
1763  pa_out = ptarray_segmentize_sphere(lwline->points, max_seg_length);
1764  return lwline_as_lwgeom(lwline_construct(lwg_in->srid, NULL, pa_out));
1765  break;
1766  case POLYGONTYPE:
1767  lwpoly_in = lwgeom_as_lwpoly(lwg_in);
1768  lwpoly_out = lwpoly_construct_empty(lwg_in->srid, lwgeom_has_z(lwg_in), lwgeom_has_m(lwg_in));
1769  for ( i = 0; i < lwpoly_in->nrings; i++ )
1770  {
1771  pa_out = ptarray_segmentize_sphere(lwpoly_in->rings[i], max_seg_length);
1772  lwpoly_add_ring(lwpoly_out, pa_out);
1773  }
1774  return lwpoly_as_lwgeom(lwpoly_out);
1775  break;
1776  case MULTILINETYPE:
1777  case MULTIPOLYGONTYPE:
1778  case COLLECTIONTYPE:
1779  lwcol_in = lwgeom_as_lwcollection(lwg_in);
1780  lwcol_out = lwcollection_construct_empty(lwg_in->type, lwg_in->srid, lwgeom_has_z(lwg_in), lwgeom_has_m(lwg_in));
1781  for ( i = 0; i < lwcol_in->ngeoms; i++ )
1782  {
1783  lwcollection_add_lwgeom(lwcol_out, lwgeom_segmentize_sphere(lwcol_in->geoms[i], max_seg_length));
1784  }
1785  return lwcollection_as_lwgeom(lwcol_out);
1786  break;
1787  default:
1788  lwerror("lwgeom_segmentize_sphere: unsupported input geometry type: %d - %s",
1789  lwg_in->type, lwtype_name(lwg_in->type));
1790  break;
1791  }
1792 
1793  lwerror("lwgeom_segmentize_sphere got to the end of the function, should not happen");
1794  return NULL;
1795 }
#define LINETYPE
Definition: liblwgeom.h:86
#define POLYGONTYPE
Definition: liblwgeom.h:87
#define MULTIPOINTTYPE
Definition: liblwgeom.h:88
LWGEOM * lwgeom_clone_deep(const LWGEOM *lwgeom)
Deep clone an LWGEOM, everything is copied.
Definition: lwgeom.c:482
LWPOLY * lwgeom_as_lwpoly(const LWGEOM *lwgeom)
Definition: lwgeom.c:174
int lwgeom_has_z(const LWGEOM *geom)
Return LW_TRUE if geometry has Z ordinates.
Definition: lwgeom.c:885
int32_t srid
Definition: liblwgeom.h:399
LWGEOM * lwpoly_as_lwgeom(const LWPOLY *obj)
Definition: lwgeom.c:288
const char * lwtype_name(uint8_t type)
Return the type name string associated with a type number (e.g.
Definition: lwutil.c:218
LWGEOM * lwline_as_lwgeom(const LWLINE *obj)
Definition: lwgeom.c:298
LWLINE * lwline_construct(int srid, GBOX *bbox, POINTARRAY *points)
Definition: lwline.c:42
LWGEOM ** geoms
Definition: liblwgeom.h:509
POINTARRAY ** rings
Definition: liblwgeom.h:457
static POINTARRAY * ptarray_segmentize_sphere(const POINTARRAY *pa_in, double max_seg_length)
Create a new point array with no segment longer than the input segment length (expressed in radians!)...
Definition: lwgeodetic.c:1680
int nrings
Definition: liblwgeom.h:455
LWGEOM * lwgeom_clone(const LWGEOM *lwgeom)
Clone LWGEOM object.
Definition: lwgeom.c:444
LWLINE * lwgeom_as_lwline(const LWGEOM *lwgeom)
Definition: lwgeom.c:138
#define MULTIPOLYGONTYPE
Definition: liblwgeom.h:90
LWGEOM * lwgeom_segmentize_sphere(const LWGEOM *lwg_in, double max_seg_length)
Create a new, densified geometry where no segment is longer than max_seg_length.
Definition: lwgeodetic.c:1739
LWCOLLECTION * lwgeom_as_lwcollection(const LWGEOM *lwgeom)
Definition: lwgeom.c:192
#define POINTTYPE
LWTYPE numbers, used internally by PostGIS.
Definition: liblwgeom.h:85
LWPOLY * lwpoly_construct_empty(int srid, char hasz, char hasm)
Definition: lwpoly.c:161
uint8_t type
Definition: liblwgeom.h:396
int lwpoly_add_ring(LWPOLY *poly, POINTARRAY *pa)
Add a ring, allocating extra space if necessary.
Definition: lwpoly.c:249
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:1346
#define MULTILINETYPE
Definition: liblwgeom.h:89
LWCOLLECTION * lwcollection_construct_empty(uint8_t type, int srid, char hasz, char hasm)
Definition: lwcollection.c:94
int lwgeom_has_m(const LWGEOM *geom)
Return LW_TRUE if geometry has M ordinates.
Definition: lwgeom.c:892
LWCOLLECTION * lwcollection_add_lwgeom(LWCOLLECTION *col, const LWGEOM *geom)
Appends geom to the collection managed by col.
Definition: lwcollection.c:187
void lwerror(const char *fmt,...)
Write a notice out to the error handler.
Definition: lwutil.c:190
#define COLLECTIONTYPE
Definition: liblwgeom.h:91
POINTARRAY * points
Definition: liblwgeom.h:422
LWGEOM * lwcollection_as_lwgeom(const LWCOLLECTION *obj)
Definition: lwgeom.c:268
Here is the call graph for this function:
Here is the caller graph for this function: