-1

地球の曲率を考慮して、半径 X (メートル単位) の円を作成する必要があります。OpenLayers 2 ではdestinationVincenty、角度を考慮してポイントから目的地を計算するのに役立つ がありました。ただし、これは OpenLayers 5 では欠落しているようです。

ジオメトリをテストしましたCricleが、不正確すぎます (画像を確認してください)。

OpenLayers 5 でこれを行うにはどうすればよいでしょうか?

青い円は半径 245km で描かれていますが、実際の半径 245km は Sedalia 付近の赤い線のあたりまでです。

青い円は半径 245km で描かれていますが、実際の半径 245km は Sedalia 付近の赤い線のあたりまでです。

4

1 に答える 1

0

destinationVincenty 関数が必要な場合は、OL2 ソースから簡単にコピーできます。3 つの明らかな依存関係を追加し、入力の lonlat と戻り値を OL3/4/5 型の座標配列に変更しましたが、何もテストしていません。

destinationVincenty = function(lonlat, brng, dist) {
    //var u = OpenLayers.Util;
    var u = {};
    u.rad = function(x) {return x*Math.PI/180;};
    u.deg = function(x) {return x*180/Math.PI;};
    u.VincentyConstants = {
        a: 6378137,
        b: 6356752.3142,
        f: 1/298.257223563
    };

    var ct = u.VincentyConstants;
    var a = ct.a, b = ct.b, f = ct.f;

    var lon1 = lonlat[0]; //lonlat.lon;
    var lat1 = lonlat[1]; //lonlat.lat;

    var s = dist;
    var alpha1 = u.rad(brng);
    var sinAlpha1 = Math.sin(alpha1);
    var cosAlpha1 = Math.cos(alpha1);

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

    var sigma = s / (b*A), sigmaP = 2*Math.PI;
    while (Math.abs(sigma-sigmaP) > 1e-12) {
        var cos2SigmaM = Math.cos(2*sigma1 + sigma);
        var sinSigma = Math.sin(sigma);
        var cosSigma = Math.cos(sigma);
        var 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;
    }

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

    var revAz = Math.atan2(sinAlpha, -tmp);  // final bearing

    //return new OpenLayers.LonLat(lon1+u.deg(L), u.deg(lat2));
    return [lon1+u.deg(L), u.deg(lat2)];
};

利用可能な最も近い方法は OL3/4/5 です ol.geom.Polygon.circular()

于 2018-12-11T20:24:32.257 に答える