PostGIS 3.0.6dev-r@@SVN_REVISION@@
Loading...
Searching...
No Matches

◆ ptarray_calculate_gbox_geodetic()

int ptarray_calculate_gbox_geodetic ( const POINTARRAY pa,
GBOX gbox 
)
extern

Calculate geodetic (x/y/z) box and add values to gbox.

Return LW_SUCCESS on success.

Definition at line 2889 of file lwgeodetic.c.

2890{
2891 uint32_t i;
2892 int first = LW_TRUE;
2893 const POINT2D *p;
2894 POINT3D A1, A2;
2895 GBOX edge_gbox;
2896
2897 assert(gbox);
2898 assert(pa);
2899
2900 gbox_init(&edge_gbox);
2901 edge_gbox.flags = gbox->flags;
2902
2903 if ( pa->npoints == 0 ) return LW_FAILURE;
2904
2905 if ( pa->npoints == 1 )
2906 {
2907 p = getPoint2d_cp(pa, 0);
2908 ll2cart(p, &A1);
2909 gbox->xmin = gbox->xmax = A1.x;
2910 gbox->ymin = gbox->ymax = A1.y;
2911 gbox->zmin = gbox->zmax = A1.z;
2912 return LW_SUCCESS;
2913 }
2914
2915 p = getPoint2d_cp(pa, 0);
2916 ll2cart(p, &A1);
2917
2918 for ( i = 1; i < pa->npoints; i++ )
2919 {
2920
2921 p = getPoint2d_cp(pa, i);
2922 ll2cart(p, &A2);
2923
2924 edge_calculate_gbox(&A1, &A2, &edge_gbox);
2925
2926 /* Initialize the box */
2927 if ( first )
2928 {
2929 gbox_duplicate(&edge_gbox, gbox);
2930 first = LW_FALSE;
2931 }
2932 /* Expand the box where necessary */
2933 else
2934 {
2935 gbox_merge(&edge_gbox, gbox);
2936 }
2937
2938 A1 = A2;
2939 }
2940
2941 return LW_SUCCESS;
2942}
int gbox_merge(const GBOX *new_box, GBOX *merge_box)
Update the merged GBOX to be large enough to include itself and the new box.
Definition gbox.c:257
void gbox_duplicate(const GBOX *original, GBOX *duplicate)
Copy the values of original GBOX into duplicate.
Definition gbox.c:433
void gbox_init(GBOX *gbox)
Zero out all the entries in the GBOX.
Definition gbox.c:40
#define LW_FALSE
Definition liblwgeom.h:108
#define LW_FAILURE
Definition liblwgeom.h:110
#define LW_SUCCESS
Definition liblwgeom.h:111
#define LW_TRUE
Return types for functions with status returns.
Definition liblwgeom.h:107
void ll2cart(const POINT2D *g, POINT3D *p)
Convert lon/lat coordinates to cartesian coordinates on unit sphere.
Definition lwgeodetic.c:423
int edge_calculate_gbox(const POINT3D *A1, const POINT3D *A2, GBOX *gbox)
The magic function, given an edge in spherical coordinates, calculate a 3D bounding box that fully co...
static const POINT2D * getPoint2d_cp(const POINTARRAY *pa, uint32_t n)
Returns a POINT2D pointer into the POINTARRAY serialized_ptlist, suitable for reading from.
Definition lwinline.h:91
double ymax
Definition liblwgeom.h:343
double zmax
Definition liblwgeom.h:345
double xmax
Definition liblwgeom.h:341
double zmin
Definition liblwgeom.h:344
double ymin
Definition liblwgeom.h:342
double xmin
Definition liblwgeom.h:340
lwflags_t flags
Definition liblwgeom.h:339
double z
Definition liblwgeom.h:388
double x
Definition liblwgeom.h:388
double y
Definition liblwgeom.h:388
uint32_t npoints
Definition liblwgeom.h:413

References edge_calculate_gbox(), GBOX::flags, gbox_duplicate(), gbox_init(), gbox_merge(), getPoint2d_cp(), ll2cart(), LW_FAILURE, LW_FALSE, LW_SUCCESS, LW_TRUE, POINTARRAY::npoints, POINT3D::x, GBOX::xmax, GBOX::xmin, POINT3D::y, GBOX::ymax, GBOX::ymin, POINT3D::z, GBOX::zmax, and GBOX::zmin.

Referenced by lwline_calculate_gbox_geodetic(), lwpoint_calculate_gbox_geodetic(), lwpolygon_calculate_gbox_geodetic(), and lwtriangle_calculate_gbox_geodetic().

Here is the call graph for this function:
Here is the caller graph for this function: