版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/qq_15642411/article/details/85116897
according latitude and longtitude compute distance
1、最大圆距离
const double earth_radius=6371;
const double pi=3.1415926;
double degtorad(double deg)
{
return ( deg *pi/ 180);
}
double GreatCircleDistance(double lat1,double lon1,double lat2, double lon2)
{
lat1 = degtorad(lat1);
lon1 = degtorad(lon1);
lat2 = degtorad(lat2);
lon2 = degtorad(lon2);
double delta_lon = lon2 - lon1;
double central_ang = acos ( sin(lat1) *sin(lat2) + cos(lat1) *cos(lat2) * cos(delta_lon) );
return earth_radius * central_ang;
}
2、Haversine 距离
double HaversineDistace(double lat1,double lon1,double lat2, double lon2)
{
lat1 = degtorad(lat1);
lon1 = degtorad(lon1);
lat2 = degtorad(lat2);
lon2 = degtorad(lon2);
double dlon=lon2-lon1;
double dlat=lat2-lat1;
double a=pow(sin(dlon/2),2)*cos(lat1)*cos(lat2)+pow(sin(dlat/2),2);
double c=2*atan2(sqrt(a),sqrt(1-a));
return earth_radius*c;
}