IOS根据两个经纬度计算相距距离

//第一种苹果自带的

1. CLLocation *orig=[[[CLLocation alloc] initWithLatitude:[mainDelegate.latitude_self doubleValue]  longitude:[mainDelegate.longitude_self doubleValue]] autorelease];
2. CLLocation* dist=[[[CLLocation alloc] initWithLatitude:[tmpNewsModel.latitude doubleValue] longitude:[tmpNewsModel.longitude doubleValue] ] autorelease];
3. CLLocationDistance kilometers=[orig distanceFromLocation:dist]/1000;
4. NSLog(@"距离:",kilometers);

//第二种，手动计算

1. #pragma mark - calculate distance  根据2个经纬度计算距离
2. #define PI 3.1415926
3. +(double) LantitudeLongitudeDist:(double)lon1 other_Lat:(double)lat1 self_Lon:(double)lon2 self_Lat:(double)lat2{
4. double er = 6378137; // 6378700.0f;
5. //ave. radius = 6371.315 (someone said more accurate is 6366.707)
6. //equatorial radius = 6378.388
7. //nautical mile = 1.15078
8. double radlat1 = PI*lat1/180.0f;
9. double radlat2 = PI*lat2/180.0f;
10. //now long.
11. double radlong1 = PI*lon1/180.0f;
12. double radlong2 = PI*lon2/180.0f;
13. if( radlat1 < 0 ) radlat1 = PI/2 + fabs(radlat1);// south
14. if( radlat1 > 0 ) radlat1 = PI/2 - fabs(radlat1);// north
15. if( radlong1 < 0 ) radlong1 = PI*2 - fabs(radlong1);//west
16. if( radlat2 < 0 ) radlat2 = PI/2 + fabs(radlat2);// south
17. if( radlat2 > 0 ) radlat2 = PI/2 - fabs(radlat2);// north
18. if( radlong2 < 0 ) radlong2 = PI*2 - fabs(radlong2);// west
19. //spherical coordinates x=r*cos(ag)sin(at), y=r*sin(ag)*sin(at), z=r*cos(at)
20. //zero ag is up so reverse lat
21. double x1 = er * cos(radlong1) * sin(radlat1);
22. double y1 = er * sin(radlong1) * sin(radlat1);
23. double z1 = er * cos(radlat1);
24. double x2 = er * cos(radlong2) * sin(radlat2);
25. double y2 = er * sin(radlong2) * sin(radlat2);
26. double z2 = er * cos(radlat2);
27. double d = sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2));
28. //side, side, side, law of cosines and arccos
29. double theta = acos((er*er+er*er-d*d)/(2*er*er));
30. double dist  = theta*er;
31. return dist;
32. }