iOS 地图根据经纬度计算距离的方法

11.经纬度计算距离

- (void)calculateTheDistance
{
    // 距离多少米
    NSString *otherlongitude = [NSString stringWithFormat:@"%@",_leftModel.longitude];
    NSString *otherlatitude = [NSString stringWithFormat:@"%@",_leftModel.latitude];
    LastLocation *lastLocation = [LocationManager lastLocation];
    double self_Lat = [lastLocation.latitude doubleValue];//32.065376; // 纬度
    double self_Lon = [lastLocation.longitude doubleValue];//118.78986; // 经度
    double distince = [self LantitudeLongitudeDist:[otherlongitude doubleValue] other_Lat:[otherlatitude doubleValue] self_Lon:self_Lon self_Lat:self_Lat];
    
    CGSize size1;
    NSString *nameStr;
    
    if (distince > 1000)
    {
        nameStr = [NSString stringWithFormat:@"%.2fkm",distince/1000];
        size1 = [nameStr sizeWithAttributes:[NSDictionary dictionaryWithObjectsAndKeys:[UIFont systemFontOfSize:14.0],NSFontAttributeName, nil]];
    }
    else{
        nameStr = [NSString stringWithFormat:@"%.0fm",distince];
        size1 = [nameStr sizeWithAttributes:[NSDictionary dictionaryWithObjectsAndKeys:[UIFont systemFontOfSize:14.0],NSFontAttributeName, nil]];
    }
    
    UILabel *distanceLab = [[UILabel alloc]initWithFrame:CGRectMake(SCREEN_SIZE.width-25-size1.width, 10,size1.width+20, size1.height)];
    distanceLab.text = [NSString stringWithFormat:@"%@",nameStr];
//    distanceLab.backgroundColor = [UIColor redColor];
    distanceLab.textAlignment = 1;
    distanceLab.font = [UIFont systemFontOfSize:15.0f];
    [_headerView addSubview:distanceLab];
}

#pragma mark - calculate distance 根据 2 个经纬度计算距离

#define PI 3.1415926
-(double) LantitudeLongitudeDist:(double)lon1 other_Lat:(double)lat1 self_Lon:(double)lon2 self_Lat:(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;///1000;
}