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(\theta_{... ...) + sin(\theta_{1})sin(\theta_{2}) \right] + sin(\phi_{1})sin(\phi_{2}) \right)$$ (19)