计算某经纬度1公里内的对象

从具体的某站到其他基站的距离,可以简单的遍历一次 得到其他站点和目标站点的距离.再做判断就可以了,已知2点坐标,获得他们的距离的方法:

2点分别为P1(Lat1, Lng1) , P2(Lat2, Lng2)


公式的单位都是弧度。

  
/** 【简注】  根据已知两点经纬度求距离.   经纬度输入以度为单位
*                     使用此函数,需要包含头文件
*     #include <math.h>
*    
* 【参数】 double lat1 第一个点的纬度
* 【参数】 double lng1 第一个点的经度
* 【参数】 double   lat2 第二个点的纬度
* 【参数】 double lng2 第二个点的经度
* 【返回】 double dm     两点间的距离,单位为千米。
*/ 
double CJwd2Doc::GetDistance(double lat1, double lng1,   double lat2, double lng2)
{  
double radlat1 = lat1*3.1415926/180.0;
double radlat2 =   lat2*3.1415926/180.0;
double radlng1 =   lng1*3.1415926/180.0;
double radlng2 = lng2*3.1415926/180.0;  

double ff = (radlat1+radlat2)/2.0;
double gg =   (radlat1-radlat2)/2.0;
double ll = (radlng1-radlng2)/2.0;  
double ss =   pow((sin(gg)),2)*pow((cos(ll)),2)+pow((cos(ff)),2)*pow((sin(ll)),2);
double   cc =   pow((cos(gg)),2)*pow((cos(ll)),2)+pow((sin(ff)),2)*pow((sin(ll)),2);
double   ww = atan(sqrt(ss/cc));  
double banjin = 6378.135;  //地球半径km
double dist =   2*ww*banjin;
double rr = sqrt(ss*cc)/ww;
double h1 =   (3*rr-1)/(2*cc);
double h2 = (3*rr+1)/(2*ss);
double xx =   1/298.257223543; //修正率
double dm;  
if(ww!=0)
dm =   dist*(1+xx*h1*pow((sin(ff)),2)*pow((cos(gg)),2)-xx*h2*pow((cos(ff)),2)*pow((sin(gg)),2));
else
dm   = 0;  
return dm;
}
阅读更多
个人分类: 算法
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭
关闭