方位角算法:
$xa,$ya A点经纬度
$xb,$yb B点经纬度
$angle=sin(($xa-$xb) * $pi/ 180.0) / (cos($yb*$pi/180.0)*tan($ya*$pi/180.0)-sin($yb*$pi/180.0)*cos(($xa-$xb)*$pi/180.0));
if($angle > 0)
{
$angle = atan($angle) * 180.0 / $pi;
if( $xb < $xa) $angle = 180 + $angle;
}
else
{
$angle = atan($angle) * 180.0 / $pi;
if( $xb < $xa) $angle = 360 + $angle;
else $angle = 180 + $angle;
}