坐标三角法 - 计算飞行路径的弧中点

时间:2014-02-04 13:26:45

标签: javascript svg d3.js trigonometry bezier

我正在尝试使用SVG在地图上绘制航线。我在Leaflet上使用d3,但使用的框架不应该对我的问题产生影响 - 它是触发器。

http://fiddle.jshell.net/zw8TR/26

我尝试这样做的方法是创建一个quadratic bezier curve(如果你知道的话,我会对其他/更简单的方式开放)。我需要计算的是1个控制点,垂直于每条线的中点。此点应始终偏向比中点更高的y值/纬度,以创建看起来像飞行路径的弧。

在上面的演示中,我发现通过添加一些额外的临时点,更容易调试控制点的确切位置。正如您所看到的,我的一些控制点面朝下,没有一个看起来是正确垂直的。

Check out my previous question on this - with diagrams

我有一种感觉问题归结为这一行:

var theta = Math.atan2(t_area_y, t_area_x) * 180 / Math.PI;

我没有正确处理负坐标。我试图破解一套令人讨厌的大门来处理这个问题。

我试图很好地评论小提琴,以解释发生了什么。一旦我知道了这一点,就应该是创建custom interpolation in d3的简单案例。

1 个答案:

答案 0 :(得分:2)

如果您使用自定义行生成器,这实际上比您想象的要容易。您只需在计算路径时添加控制点,而不是将控制点添加到要素中。代码如下所示:

feature.attr("d", function(d) {
    var s, prev;
    d.geometry.coordinates.forEach(function(c) {
        var proj = map.latLngToLayerPoint(new L.LatLng(c[1], c[0]));
        if(s) {
            var length = Math.sqrt(Math.pow(proj.x - prev.x, 2), Math.pow(proj.y - prev.y, 2)),
                midx = prev.x + (proj.x - prev.x) / 2,
                midy = prev.y + (proj.y - prev.y) / 2 - length * 0.2 * (Math.abs(proj.x - prev.x) / length);
            s += "Q" + midx + "," + midy + " " + proj.x + "," + proj.y;
        } else {
            s = "M" + proj.x + "," + proj.y;
        }
        prev = proj;
    });
    return s;
});

让我们一步一步地完成它。最重要的是,我正在跟踪前一点的坐标,以便能够计算控制点。首先,s将为空并且采用else分支 - 只需移动到该点(起点)而不绘制线条。对于所有后续点,将进行实际计算。

首先,我们计算两点(先前和当前)之间的距离length。计算控制点的x坐标很简单,因为不需要偏移。 y坐标有点棘手 - 第一部分是相同的,然后添加偏移量。偏移的大小是此处路径长度的20%(为更长的路径制作更宽的弧),根据需要进行调整。这需要乘以x轴的角度余弦,但幸运的是我们不需要明确计算角度 - 它是由点之间的距离和差异之间的关系定义的。 x坐标(该角度的反余弦)。所以我们可以直接采用这种关系并乘以它。如果您希望弧指向(即负y偏移),我们将采用x坐标差异的绝对值。如果没有这个,一些控制点就会指向下方。

完整示例here