vector_angle()

double vector_angle	(	const POINT3D *	v1,
		const POINT3D *	v2
	)

Angle between two unit vectors.

Definition at line 505 of file lwgeodetic.c.

{
        POINT3D v3, normal;
        double angle, x, y;
 
        cross_product(v1, v2, &normal);
        normalize(&normal);
        cross_product(&normal, v1, &v3);
 
        x = dot_product(v1, v2);
        y = dot_product(v2, &v3);
 
        angle = atan2(y, x);
        return angle;
}

References cross_product(), dot_product(), normalize(), pixval::x, and pixval::y.

Referenced by gbox_geocentric_get_gbox_cartesian(), and test_vector_angle().

Here is the call graph for this function:

Here is the caller graph for this function: