经纬度计算

define(EARTH_RADIUS, 6371);//地球半径，平均半径为6371km

 /**

 *计算某个经纬度的周围某段距离的正方形的四个点

 *

 *@param lng float 经度

 *@param lat float 纬度

 *@param distance float 该点所在圆的半径，该圆与此正方形内切，默认值为0.5千米

 *@return array 正方形的四个点的经纬度坐标

 */



 $lng = 117.30876;

 $lat = 31.8636;



 function returnSquarePoint($lng, $lat,$distance = 0.5){



    $dlng =  2 * asin(sin($distance / (2 * EARTH_RADIUS)) / cos(deg2rad($lat)));

    $dlng = rad2deg($dlng);



    $dlat = $distance/EARTH_RADIUS;

    $dlat = rad2deg($dlat);



    return array(

                'left-top'=>array('lat'=>$lat + $dlat,'lng'=>$lng-$dlng),

                'right-top'=>array('lat'=>$lat + $dlat, 'lng'=>$lng + $dlng),

                'left-bottom'=>array('lat'=>$lat - $dlat, 'lng'=>$lng - $dlng),

                'right-bottom'=>array('lat'=>$lat - $dlat, 'lng'=>$lng + $dlng)

                );

 }

//使用此函数计算得到结果后，带入sql查询。

$squares = returnSquarePoint($lng, $lat);



print_r($squares);

//$info_sql = "select id,locateinfo,lat,lng from `lbs_info` where lat<>0 and lat>{$squares['right-bottom']['lat']} and lat<{$squares['left-top']['lat']} and lng>{$squares['left-top']['lng']} and lng<{$squares['right-bottom']['lng']} ";

SELECT store_id,map_lng_b,map_lat_b, ROUND(6378.138*2*ASIN(SQRT(POW(SIN((31.86486929872*PI()/180-map_lat_b*PI()/180)/2),2)+COS(31.86486929872*PI()/180)*COS(map_lat_b*PI()/180)*POW(SIN((117.30853228912*PI()/180-map_lng_b*PI()/180)/2),2)))*1000) AS juli

FROM hfecm_store

where map_lat_b >0 AND map_lng_b > 0 

ORDER BY juli ASC

LIMIT 20

 

今天遇到计算经纬度的问题，百度到的，记录下。