PostGIS  2.5.0dev-r@@SVN_REVISION@@
double ptarray_length_spheroid ( const POINTARRAY pa,
const SPHEROID s 
)

Definition at line 3143 of file lwgeodetic.c.

References SPHEROID::a, SPHEROID::b, POINTARRAY::flags, FLAGS_GET_Z, geographic_point_init(), getPoint4d_p(), LW_FALSE, POINTARRAY::npoints, SPHEROID::radius, sphere_distance(), spheroid_distance(), POINT4D::x, POINT4D::y, and POINT4D::z.

Referenced by lwgeom_length_spheroid().

3144 {
3145  GEOGRAPHIC_POINT a, b;
3146  double za = 0.0, zb = 0.0;
3147  POINT4D p;
3148  uint32_t i;
3149  int hasz = LW_FALSE;
3150  double length = 0.0;
3151  double seglength = 0.0;
3152 
3153  /* Return zero on non-sensical inputs */
3154  if ( ! pa || pa->npoints < 2 )
3155  return 0.0;
3156 
3157  /* See if we have a third dimension */
3158  hasz = FLAGS_GET_Z(pa->flags);
3159 
3160  /* Initialize first point */
3161  getPoint4d_p(pa, 0, &p);
3162  geographic_point_init(p.x, p.y, &a);
3163  if ( hasz )
3164  za = p.z;
3165 
3166  /* Loop and sum the length for each segment */
3167  for ( i = 1; i < pa->npoints; i++ )
3168  {
3169  seglength = 0.0;
3170  getPoint4d_p(pa, i, &p);
3171  geographic_point_init(p.x, p.y, &b);
3172  if ( hasz )
3173  zb = p.z;
3174 
3175  /* Special sphere case */
3176  if ( s->a == s->b )
3177  seglength = s->radius * sphere_distance(&a, &b);
3178  /* Spheroid case */
3179  else
3180  seglength = spheroid_distance(&a, &b, s);
3181 
3182  /* Add in the vertical displacement if we're in 3D */
3183  if ( hasz )
3184  seglength = sqrt( (zb-za)*(zb-za) + seglength*seglength );
3185 
3186  /* Add this segment length to the total */
3187  length += seglength;
3188 
3189  /* B gets incremented in the next loop, so we save the value here */
3190  a = b;
3191  za = zb;
3192  }
3193  return length;
3194 }
double x
Definition: liblwgeom.h:351
double sphere_distance(const GEOGRAPHIC_POINT *s, const GEOGRAPHIC_POINT *e)
Given two points on a unit sphere, calculate their distance apart in radians.
Definition: lwgeodetic.c:913
double b
Definition: liblwgeom.h:313
double radius
Definition: liblwgeom.h:317
Point in spherical coordinates on the world.
Definition: lwgeodetic.h:47
unsigned int uint32_t
Definition: uthash.h:78
#define LW_FALSE
Definition: liblwgeom.h:76
uint8_t flags
Definition: liblwgeom.h:368
#define FLAGS_GET_Z(flags)
Macros for manipulating the 'flags' byte.
Definition: liblwgeom.h:139
double z
Definition: liblwgeom.h:351
int getPoint4d_p(const POINTARRAY *pa, uint32_t n, POINT4D *point)
Definition: lwgeom_api.c:113
double a
Definition: liblwgeom.h:312
void geographic_point_init(double lon, double lat, GEOGRAPHIC_POINT *g)
Initialize a geographic point.
Definition: lwgeodetic.c:171
double spheroid_distance(const GEOGRAPHIC_POINT *a, const GEOGRAPHIC_POINT *b, const SPHEROID *spheroid)
Computes the shortest distance along the surface of the spheroid between two points.
Definition: lwspheroid.c:186
double y
Definition: liblwgeom.h:351
uint32_t npoints
Definition: liblwgeom.h:370

Here is the call graph for this function:

Here is the caller graph for this function: