spheroid_project()

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.

Using the direct geodesic problem from GeographicLib (Karney 2013).

Parameters

r	- location of first point
distance	- distance in meters
azimuth	- azimuth in radians

Returns

g - location of projected point

Definition at line 128 of file lwspheroid.c.

{
        struct geod_geodesic gd;
        geod_init(&gd, spheroid->a, spheroid->f);
        double lat1 = r->lat * 180.0 / M_PI;
        double lon1 = r->lon * 180.0 / M_PI;
        double lat2, lon2; /* return projected position */
        geod_direct(&gd, lat1, lon1, azimuth * 180.0 / M_PI, distance, &lat2, &lon2, 0);
        g->lat = lat2 * M_PI / 180.0;
        g->lon = lon2 * M_PI / 180.0;
        return LW_SUCCESS;
}

References SPHEROID::a, distance(), SPHEROID::f, GEOGRAPHIC_POINT::lat, GEOGRAPHIC_POINT::lon, LW_SUCCESS, and r.

Referenced by lwgeom_project_spheroid().

Here is the call graph for this function:

Here is the caller graph for this function: