经纬度计算
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
今天遇到计算经纬度的问题,百度到的,记录下。