php根据经纬度计算距离和方向--摘录自http://haotushu.sinaapp.com/post-520.html

define('EARTH_RADIUS', 6367000);//需定义的静态变量

function getRadian($d)



{



return $d * M_PI / 180;



}



function getDistance ($lat1, $lng1, $lat2, $lng2)



{



$lat1 =getRadian($lat1);



$lat2 = getRadian($lat2);



$a = $lat1 - $lat2;



$b = getRadian($lng1) - getRadian($lng2);



$v = 2 * asin(sqrt(pow(sin($a/2),2) +



cos($lat1) * cos($lat2) * pow(sin($b/2),2)));



$v = round(EARTH_RADIUS * $v * 10000) / 10000;



return $v;



}

//根据经纬度计算方向



function getDirection($lat1, $lng1, $lat2, $lng2){



$k1 = $lng2-$lng1;



$k2 = $lat2-$lat1;



if( 0 == $k1){



if($k2>0){



$str="聊友在您的正北方 ";



}



else if( $k2<0){



$str ="聊友在您的正南方 ";



}



else if( $k2 == 0){



$str="聊友正在您的附近 ";



}



}else if( 0 == $k2){



if($k1>0){



$str="聊友在您的正东方 ";



}



else if( $k1<0){



$str="聊友在您的正西方 ";



}



}else{



$k=$k2/$k1;



if($k2>0){



if($k1>0){



$angle = 180*atan($k)/M_PI;



$str="聊友在您的东偏北 $angle 度方向  ";



}else if($k1<0){



$angle = 180*atan(-$k)/M_PI;



$str= "聊友在您的西偏北 $angle 度方向  ";



}



}else if($k2<0){



if($k1<0){



$angle = 180*atan($k)/M_PI;



$str = "聊友在您的西偏南 $angle 度方向  ";



}



else if($k1>0){



$angle = 180*atan($k)/M_PI;



$str="聊友在您的东偏南 $angle 度方向  ";



}



}



}



return $str;



}

echo getDirection(40.770623,-73.964367,40.758224,-73.917404);



exit;

结果