5

我正在使用 Fortune 算法的 Javascript 实现来计算 voronoi 单元（https://github.com/gorhill/Javascript-Voronoi）。我要计算的站点是地图上的点（所以(lat,lng)）。我首先进行了投影(lat,lng) -> (x,y)，然后计算了 voronoi 细胞，并以另一种方式进行了半边的投影。
它工作正常，我使用传单显示结果，但我还需要做一件事。

我最初计算的每个站点都依赖于一个 ID，我按 ID 重新分类 voronoi 单元，最后，每个 ID 都具有如下所示的标准数据结构：

{
  "type": "FeatureCollection",
  "features": [
    {
      "type": "Feature",
      "geometry": {
        "type": "Polygon",
        "coordinates": [[
          [9.994812, 53.549487],
          [10.046997, 53.598209],
          [10.117721, 53.531737],
          [9.994812, 53.549487]
        ]]
      }
    }, {
      "type": "Feature",
      "geometry": {
        "type": "Polygon",
        "coordinates": [[
          [10.000991, 53.50418],
          [10.03807, 53.562539],
          [9.926834, 53.551731],
          [10.000991, 53.50418]
        ]]
      }
    }
  ]
};

给定 ID 的一组多边形（由 voronoi 单元的半边构成）。

我需要按 ID 合并这些多边形，我打算使用turf.merge()，但我有拓扑错误

 turf.min.js:13 Uncaught TopologyError: side location conflict

根据这篇文章（http://lists.refracions.net/pipermail/jts-devel/2009-March/002939.html），我尝试将(lat,lng)这对夫妇从 10^-14 到 10^-7 但它并没有真正奏效。在寻找扭结并尝试消除它们之前，我打印了一些数据样本，并且我知道我是否使用了来自 Fortune 算法的好数据。当我显示所有 ID 的所有多边形时，我得到了正确的图表，但是当我显示一个 ID 的所有多边形或一个 ID 的一些多边形时，我最终得到了不完整的图表：

完整图表的一部分

完整图表的一部分

在此处输入图像描述

一个 ID 的部分图表

在此处输入图像描述

给定 ID 的两个“多边形”

有谁知道如何合并共享至少一个公共顶点的多边形？为什么会出现拓扑错误？

编辑：多边形不是“不完整的”（我使用的是折线）

在此处输入图像描述

我还尝试了一个更简单的示例：

在此处输入图像描述

仍然得到错误：

Uncaught TopologyError: side location conflict [ (44.8220601, -0.5869532) ]

所以这不是（或至少不仅仅是）由于扭结

4

1 回答 1

2

在数据到达 Turf 之前，您的问题似乎已经发生。通过GeoJSON 验证器从您的GitHub 问题运行 GeoJSON 会发现两个错误。第一个是每个特征只包含一个对象，GeoJSON 要求所有特征也有一个对象，即使它是空的。其次，更重要的是，一个有效的 GeoJSON 多边形必须是一个闭环，第一个点和最后一个点的坐标相同。第二个问题似乎是导致 Turf 抛出错误的原因。一旦将第一组坐标复制到末端以闭合环，多边形将成功合并。geometryproperties

在地图上显示数据后，您的纬度和经度也很明显颠倒了。坐标应该lon,lat在 GeoJSON 中，因为你的坐标在 中lat,lon，所以多边形出现在印度洋的中间。一旦纠正，它们就会出现在正确的位置。

这是一个显示他们成功合并的小提琴：

http://fiddle.jshell.net/nathansnider/p7kfxvk7/

于 2016-03-30T09:06:33.647 回答