PostGIS  2.5.0dev-r@@SVN_REVISION@@
LWPOINT* lwgeom_project_spheroid ( const LWPOINT r,
const SPHEROID spheroid,
double  distance,
double  azimuth 
)

Calculate the location of a point on a spheroid, give a start point, bearing and distance.

Calculate the location of a point on a spheroid, give a start point, bearing and distance.

Parameters
r- location of first point.
spheroid- spheroid definition.
distance- distance, in units of the spheroid def'n.
azimuth- azimuth in radians.
Returns
s - location of projected point.

Definition at line 2008 of file lwgeodetic.c.

References distance(), geographic_point_init(), GEOGRAPHIC_POINT::lat, latitude_radians_normalize(), GEOGRAPHIC_POINT::lon, longitude_radians_normalize(), LW_FAILURE, LW_TRUE, LWDEBUGF, lwerror(), lwgeom_set_geodetic(), lwpoint_as_lwgeom(), lwpoint_construct(), lwpoint_get_x(), lwpoint_get_y(), POINT4D::m, ptarray_construct(), ptarray_set_point4d(), rad2deg, SPHEROID::radius, spheroid_project(), LWPOINT::srid, pixval::x, POINT4D::x, pixval::y, POINT4D::y, and POINT4D::z.

Referenced by geography_project().

2009 {
2010  GEOGRAPHIC_POINT geo_source, geo_dest;
2011  POINT4D pt_dest;
2012  double x, y;
2013  POINTARRAY *pa;
2014  LWPOINT *lwp;
2015 
2016  /* Normalize distance to be positive*/
2017  if ( distance < 0.0 ) {
2018  distance = -distance;
2019  azimuth += M_PI;
2020  }
2021 
2022  /* Normalize azimuth */
2023  azimuth -= 2.0 * M_PI * floor(azimuth / (2.0 * M_PI));
2024 
2025  /* Check the distance validity */
2026  if ( distance > (M_PI * spheroid->radius) )
2027  {
2028  lwerror("Distance must not be greater than %g", M_PI * spheroid->radius);
2029  return NULL;
2030  }
2031 
2032  /* Convert to ta geodetic point */
2033  x = lwpoint_get_x(r);
2034  y = lwpoint_get_y(r);
2035  geographic_point_init(x, y, &geo_source);
2036 
2037  /* Try the projection */
2038  if( spheroid_project(&geo_source, spheroid, distance, azimuth, &geo_dest) == LW_FAILURE )
2039  {
2040  LWDEBUGF(3, "Unable to project from (%g %g) with azimuth %g and distance %g", x, y, azimuth, distance);
2041  lwerror("Unable to project from (%g %g) with azimuth %g and distance %g", x, y, azimuth, distance);
2042  return NULL;
2043  }
2044 
2045  /* Build the output LWPOINT */
2046  pa = ptarray_construct(0, 0, 1);
2047  pt_dest.x = rad2deg(longitude_radians_normalize(geo_dest.lon));
2048  pt_dest.y = rad2deg(latitude_radians_normalize(geo_dest.lat));
2049  pt_dest.z = pt_dest.m = 0.0;
2050  ptarray_set_point4d(pa, 0, &pt_dest);
2051  lwp = lwpoint_construct(r->srid, NULL, pa);
2053  return lwp;
2054 }
double x
Definition: liblwgeom.h:351
POINTARRAY * ptarray_construct(char hasz, char hasm, uint32_t npoints)
Construct an empty pointarray, allocating storage and setting the npoints, but not filling in any inf...
Definition: ptarray.c:62
double longitude_radians_normalize(double lon)
Convert a longitude to the range of -PI,PI.
Definition: lwgeodetic.c:41
double m
Definition: liblwgeom.h:351
void ptarray_set_point4d(POINTARRAY *pa, uint32_t n, const POINT4D *p4d)
Definition: lwgeom_api.c:425
double radius
Definition: liblwgeom.h:317
Point in spherical coordinates on the world.
Definition: lwgeodetic.h:47
int spheroid_project(const GEOGRAPHIC_POINT *r, const SPHEROID *spheroid, double distance, double azimuth, GEOGRAPHIC_POINT *g)
Given a location, an azimuth and a distance, computes the location of the projected point...
Definition: lwspheroid.c:359
#define LW_FAILURE
Definition: liblwgeom.h:78
double lwpoint_get_x(const LWPOINT *point)
Definition: lwpoint.c:63
#define rad2deg(r)
Definition: lwgeodetic.h:75
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:75
double z
Definition: liblwgeom.h:351
tuple x
Definition: pixval.py:53
Datum distance(PG_FUNCTION_ARGS)
void lwgeom_set_geodetic(LWGEOM *geom, int value)
Set the FLAGS geodetic bit on geometry an all sub-geometries and pointlists.
Definition: lwgeom.c:945
int32_t srid
Definition: liblwgeom.h:409
double latitude_radians_normalize(double lat)
Convert a latitude to the range of -PI/2,PI/2.
Definition: lwgeodetic.c:69
void geographic_point_init(double lon, double lat, GEOGRAPHIC_POINT *g)
Initialize a geographic point.
Definition: lwgeodetic.c:171
double lwpoint_get_y(const LWPOINT *point)
Definition: lwpoint.c:76
LWGEOM * lwpoint_as_lwgeom(const LWPOINT *obj)
Definition: lwgeom.c:334
LWPOINT * lwpoint_construct(int srid, GBOX *bbox, POINTARRAY *point)
Definition: lwpoint.c:129
double y
Definition: liblwgeom.h:351
#define LWDEBUGF(level, msg,...)
Definition: lwgeom_log.h:88
void lwerror(const char *fmt,...)
Write a notice out to the error handler.
Definition: lwutil.c:190
tuple y
Definition: pixval.py:54

Here is the call graph for this function:

Here is the caller graph for this function: