MapKit 地图上两点之间的距离(源码函数)
2014-08-01 10:57
459 查看
MapKit 地图上两点之间的距离(源码函数)
/***********************************************************
传入参数两点 及两对 latitude 、longitude,则返回以米为单位的 double 数值
/***********************************************************
#define PI 3.1415926
double LantitudeLongitudeDist(double lon1,double lat1,
double lon2,double lat2)
{
double er = 6378137; // 6378700.0f;
//ave. radius = 6371.315 (someone said more accurate is 6366.707)
//equatorial radius = 6378.388
//nautical mile = 1.15078
double radlat1 = PI*lat1/180.0f;
double radlat2 = PI*lat2/180.0f;
//now long.
double radlong1 = PI*lon1/180.0f;
double radlong2 = PI*lon2/180.0f;
if( radlat1 < 0 ) radlat1 = PI/2 + fabs(radlat1);// south
if( radlat1 > 0 ) radlat1 = PI/2 - fabs(radlat1);// north
if( radlong1 < 0 ) radlong1 = PI*2 - fabs(radlong1);//west
if( radlat2 < 0 ) radlat2 = PI/2 + fabs(radlat2);// south
if( radlat2 > 0 ) radlat2 = PI/2 - fabs(radlat2);// north
if( radlong2 < 0 ) radlong2 = PI*2 - fabs(radlong2);// west
//spherical coordinates x=r*cos(ag)sin(at), y=r*sin(ag)*sin(at), z=r*cos(at)
//zero ag is up so reverse lat
double x1 = er * cos(radlong1) * sin(radlat1);
double y1 = er * sin(radlong1) * sin(radlat1);
double z1 = er * cos(radlat1);
double x2 = er * cos(radlong2) * sin(radlat2);
double y2 = er * sin(radlong2) * sin(radlat2);
double z2 = er * cos(radlat2);
double d = sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2));
//side, side, side, law of cosines and arccos
double theta = acos((er*er+er*er-d*d)/(2*er*er));
double dist = theta*er;
return dist;//return 米
}
/***********************************************************
传入参数两点 及两对 latitude 、longitude,则返回以米为单位的 double 数值
/***********************************************************
#define PI 3.1415926
double LantitudeLongitudeDist(double lon1,double lat1,
double lon2,double lat2)
{
double er = 6378137; // 6378700.0f;
//ave. radius = 6371.315 (someone said more accurate is 6366.707)
//equatorial radius = 6378.388
//nautical mile = 1.15078
double radlat1 = PI*lat1/180.0f;
double radlat2 = PI*lat2/180.0f;
//now long.
double radlong1 = PI*lon1/180.0f;
double radlong2 = PI*lon2/180.0f;
if( radlat1 < 0 ) radlat1 = PI/2 + fabs(radlat1);// south
if( radlat1 > 0 ) radlat1 = PI/2 - fabs(radlat1);// north
if( radlong1 < 0 ) radlong1 = PI*2 - fabs(radlong1);//west
if( radlat2 < 0 ) radlat2 = PI/2 + fabs(radlat2);// south
if( radlat2 > 0 ) radlat2 = PI/2 - fabs(radlat2);// north
if( radlong2 < 0 ) radlong2 = PI*2 - fabs(radlong2);// west
//spherical coordinates x=r*cos(ag)sin(at), y=r*sin(ag)*sin(at), z=r*cos(at)
//zero ag is up so reverse lat
double x1 = er * cos(radlong1) * sin(radlat1);
double y1 = er * sin(radlong1) * sin(radlat1);
double z1 = er * cos(radlat1);
double x2 = er * cos(radlong2) * sin(radlat2);
double y2 = er * sin(radlong2) * sin(radlat2);
double z2 = er * cos(radlat2);
double d = sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2));
//side, side, side, law of cosines and arccos
double theta = acos((er*er+er*er-d*d)/(2*er*er));
double dist = theta*er;
return dist;//return 米
}
相关文章推荐
- iOS开发MapKit地图两点之间的距离及线路的绘制
- 计算指定两点经纬度之间的距离(附PHP源码)
- 计算地图上两点之间的距离
- 计算地图上两点之间的距离
- 第七周实验报告(二)利用成员函数、友元函数和一般函数求两点之间距离
- 根据经纬度计算两点之间的距离的公式推导过程以及google.maps的测距函数
- 高德地图设置中心点和缩放比例,获取两点之间距离
- 用友元函数设计点类求两点之间的距离
- 第七周项目 用友元函数求两点之间的距离
- 如何计算地图上两点之间的距离
- 地图两点之间获取距离公式
- 算法竞赛入门经典: 第四章 函数与递归 4.1求两点之间距离
- 用成员函数设计点类求两点之间的距离
- 高德地图计算两点之间的距离java+html
- 第七周项目 用一般函数,成员函数,友元函数求两点之间的距离
- 第七周项目 求两点之间的距离 成员函数
- 根据高德地图得出的坐标算出两点之间的距离 java
- 用一般函数设计点类求两点之间的距离
- 地图-获取两点之间的距离(经纬度)
- 第七周项目 用一般函数求两点之间的距离