PostGIS  2.5.1dev-r@@SVN_REVISION@@

◆ robust_cross_product()

void robust_cross_product ( const GEOGRAPHIC_POINT p,
const GEOGRAPHIC_POINT q,
POINT3D a 
)

Computes the cross product of two vectors using their lat, lng representations.

Good even for small distances between p and q.

Definition at line 603 of file lwgeodetic.c.

References GEOGRAPHIC_POINT::lat, GEOGRAPHIC_POINT::lon, POINT3D::x, POINT3D::y, and POINT3D::z.

Referenced by clairaut_geographic(), edge_distance_to_point(), edge_intersection(), edge_point_side(), sphere_angle(), and test_sphere_project().

604 {
605  double lon_qpp = (q->lon + p->lon) / -2.0;
606  double lon_qmp = (q->lon - p->lon) / 2.0;
607  double sin_p_lat_minus_q_lat = sin(p->lat-q->lat);
608  double sin_p_lat_plus_q_lat = sin(p->lat+q->lat);
609  double sin_lon_qpp = sin(lon_qpp);
610  double sin_lon_qmp = sin(lon_qmp);
611  double cos_lon_qpp = cos(lon_qpp);
612  double cos_lon_qmp = cos(lon_qmp);
613  a->x = sin_p_lat_minus_q_lat * sin_lon_qpp * cos_lon_qmp -
614  sin_p_lat_plus_q_lat * cos_lon_qpp * sin_lon_qmp;
615  a->y = sin_p_lat_minus_q_lat * cos_lon_qpp * cos_lon_qmp +
616  sin_p_lat_plus_q_lat * sin_lon_qpp * sin_lon_qmp;
617  a->z = cos(p->lat) * cos(q->lat) * sin(q->lon-p->lon);
618 }
double y
Definition: liblwgeom.h:342
double x
Definition: liblwgeom.h:342
double z
Definition: liblwgeom.h:342
Here is the caller graph for this function: