iOS：从起点给定距离和方位的目标点答案

【问题标题】：iOS: Destination point given distance and bearing from start pointiOS：从起点给定距离和方位的目标点
【发布时间】：2011-04-12 15:52:28
【问题描述】：

给定起点、初始方位和距离，这将计算沿（最短距离）大圆弧行进的目的地点和最终方位：

var lat2 = 

Math.asin( Math.sin(lat1)*Math.cos(d/R) +                       
Math.cos(lat1)*Math.sin(d/R)*Math.cos(brng) );

var lon2 = 

lon1 + Math.atan2(Math.sin(brng)*Math.sin(d/R)*Math.cos(lat1),  
       Math.cos(d/R)-Math.sin(lat1)*Math.sin(lat2));

此代码使用 JavaScript。我希望 iOS 也能做同样的事情，所以在 Objective-c 中。

任何人都知道一个可以解决问题的类吗？

【问题讨论】：

亚马逊书店里有一些关于 C 和 Objective-C 的好书。我建议买一个并学习语言。除此之外，看看 math.h

标签： javascript objective-c ios geolocation

【解决方案1】：

我不是在翻译你的 Javascript，这只是我在某处为我的一个项目找到的一个例程，它做同样的事情：

- (CLLocationCoordinate2D) NewLocationFrom:(CLLocationCoordinate2D)startingPoint 
                         atDistanceInMiles:(float)distanceInMiles 
                     alongBearingInDegrees:(double)bearingInDegrees {

    double lat1 = DEG2RAD(startingPoint.latitude);
    double lon1 = DEG2RAD(startingPoint.longitude);

    double a = 6378137, b = 6356752.3142, f = 1/298.257223563;  // WGS-84 ellipsiod
    double s = distanceInMiles * 1.61 * 1000;  // Convert to meters
    double alpha1 = DEG2RAD(bearingInDegrees);
    double sinAlpha1 = sin(alpha1);
    double cosAlpha1 = cos(alpha1);

    double tanU1 = (1 - f) * tan(lat1);
    double cosU1 = 1 / sqrt((1 + tanU1 * tanU1));
    double sinU1 = tanU1 * cosU1;
    double sigma1 = atan2(tanU1, cosAlpha1);
    double sinAlpha = cosU1 * sinAlpha1;
    double cosSqAlpha = 1 - sinAlpha * sinAlpha;
    double uSq = cosSqAlpha * (a * a - b * b) / (b * b);
    double A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq)));
    double B = uSq / 1024 * (256 + uSq * (-128 + uSq * (74 - 47 * uSq)));

    double sigma = s / (b * A);
    double sigmaP = 2 * kPi;

    double cos2SigmaM;
    double sinSigma;
    double cosSigma;

    while (abs(sigma - sigmaP) > 1e-12) {
        cos2SigmaM = cos(2 * sigma1 + sigma);
        sinSigma = sin(sigma);
        cosSigma = cos(sigma);
        double deltaSigma = B * sinSigma * (cos2SigmaM + B / 4 * (cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM) - B / 6 * cos2SigmaM * (-3 + 4 * sinSigma * sinSigma) * (-3 + 4 * cos2SigmaM * cos2SigmaM)));
        sigmaP = sigma;
        sigma = s / (b * A) + deltaSigma;
    }

    double tmp = sinU1 * sinSigma - cosU1 * cosSigma * cosAlpha1;
    double lat2 = atan2(sinU1 * cosSigma + cosU1 * sinSigma * cosAlpha1, (1 - f) * sqrt(sinAlpha * sinAlpha + tmp * tmp));
    double lambda = atan2(sinSigma * sinAlpha1, cosU1 * cosSigma - sinU1 * sinSigma * cosAlpha1);
    double C = f / 16 * cosSqAlpha * (4 + f * (4 - 3 * cosSqAlpha));
    double L = lambda - (1 - C) * f * sinAlpha * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM)));

    double lon2 = lon1 + L;

    // Create a new CLLocationCoordinate2D for this point
    CLLocationCoordinate2D edgePoint = CLLocationCoordinate2DMake(RAD2DEG(lat2), RAD2DEG(lon2));

    return edgePoint;
}

【讨论】：

这个答案中的计算比 OP 中给出的公式更复杂，但它们似乎确实有效。知道是否有关于此功能背后逻辑的文章吗？
@D-Nice 我已经保存了一些有用的网站，我曾经用它来完成这个，包括一个以伪代码提供它并解释了它的大部分工作原理的网站，但不幸的是，我似乎已经删除了它们。如果我追踪到它们，我会发布它们，但可能性不大（从备份中检索它们比我现在可以回答的工作要多一些）。
这个计算是 Vincenty 公式，它会比 OP 中的 Haversine 公式更准确。 Vincenty 公式精确到 0.5 毫米（这非常疯狂），Haversine 不太准确但更快。但是，既然你不是手工做的，不妨选择 Vincenty。[链接]en.wikipedia.org/wiki/Vincenty%27s_formulae
迄今为止我给予的最当之无愧的支持。干得好，先生，干得好！

【解决方案2】：

你让事情变得非常困难:)

试试我的代码：

LatDistance = Cos(BearingInDegrees)*distance;
LongDistance = Sin(BearingInDegrees)*distance;

CLLocationDegrees *destinationLat = CurrentLatitude + (LatDistance*0.00001);
CLLocationDegrees *destinationLong = CurrentLongitude + (LongDistance*0.00001);

就是这样。很简单。

【讨论】：