angdist.m

Syntax:

psi = angdist(phi1,theta1,phi2,theta2)

This function calculates the angular distance between two points on a sphere. Given (latitude, longitude) for two points $(\phi_{1},\theta_{1})$ and $(\phi_{2},\theta_{2})$, the angle $\psi$ between them is:

$$\psi = \arccos\left( cos(\phi_{1})\cos(\phi_{2})\left[ \cos(\thet... ...\sin(\theta_{1})\sin(\theta_{2}) \right] + \sin(\phi_{1})\sin(\phi_{2}) \right)$$ (24)