一组多边形的外包络答案

【问题标题】：Outer envelope to a group of polygons一组多边形的外包络
【发布时间】：2018-02-13 13:24:01
【问题描述】：

我有一个包含多个多边形的 GeoJson 文件。像这样。

我使用 Leaflet 在网站中呈现此 GeoJson。
我想在包围所有多边形的多边形周围画一个轮廓。像这样的东西。

我正在使用的 GeoJSOn 格式：

{
"features": [
  {
    "geometry": {
      "coordinates": [
        [
          [
            1074.426,
            -1136.986
          ],
          [
            1088.241,
            -1123.171
          ]
        ]
      ],
      "type": "Polygon"
    },
    "properties": {
      "number": "2009",
      "type": "",
      "spaceid": null,
      "alias": null,
      "roomkey": "5/2009"
    },
    "type": "Feature"
  }
],
"bbox": [
  2445.578,
  2445.578
],
"crs": {
  "properties": {
    "name": "urn:ogc:def:crs:OGC:1.3:CRS84"
  },
  "type": "name"
},
"type": "FeatureCollection"

}

任何指针都会有所帮助:) 谢谢

【问题讨论】：

您希望包围的多边形是否共享边？如果没有，凸包是您的主要选择（如另一个答案中所述）。如果您的多边形共享边，则有一些潜在的解决方案可以根据未共享的边（以及组合多边形的外部限制）组装外部多边形（和内部环）。

标签： javascript leaflet polygon geojson turfjs

【解决方案1】：

我设法做到了，但使用的是 TopoJSON。 I post a video here.

首先我尝试在 tuffjs 中使用“联合”功能。但是它非常慢（我有一个非常详细的geojson）。

所以我切换到了 topojson。首先，您需要将 geojson 转换为 topojson。然后使用 topojson-client 库中的“合并”功能。合并返回多面体，但它比原始几何体简单得多。

然后你需要在代码中做额外的处理，移除一些属于其他多边形的多边形。

【讨论】：

【解决方案2】：

@Shaswat，关于凸包缺少内点的说法是正确的。所以我尝试了 turf.union：

const turf = require('@turf/turf')

const originGeojson = require('./SECCIONES_13_geo.json')

const totalUnion = originGeojson.features.reduce((union, feature, index) => {
  if (!union) {
    return turf.polygon(feature.geometry.coordinates)
  }

  try {
    return turf.union(union, turf.polygon(feature.geometry.coordinates))
  } catch (err) {
    return union
  }
})

console.log(JSON.stringify(totalUnion))

但它会产生这样的东西，里面有很多洞。

代码本身是不正确的，对于catch块，这只是一种遍历整个列表的方法。 catch 中的错误是：

错误：多边形的每个线性环必须有 4 个或更多位置。

如果有人能分享解决此问题的正确方法，我将不胜感激。

【讨论】：

【解决方案3】：

您正在寻找“凸包”：

在数学中，在欧几里得平面或欧几里得空间（或更一般地，在实数上的仿射空间）中的一组点 X 的凸包或凸包络是包含 X 的最小凸集。

参考：https://en.wikipedia.org/wiki/Convex_hull

您可以使用 Turf.js convex 方法做到这一点：

采用 Feature 或 FeatureCollection 并返回凸包 Polygon。

参考：http://turfjs.org/docs/#convex

例子：

var map = new L.Map('leaflet', {center: [0, 0], zoom: 0});

var collection = turf.featureCollection([
    turf.polygon([[[-80,-80],[-40,-80],[-40,-40],[-80,-40],[-80,-80]]]),
    turf.polygon([[[80,80],[40,80],[40,40],[80,40],[80,80]]])
]);

new L.GeoJSON(collection, {color: 'red'}).addTo(map);

var polygon = turf.convex(collection);

new L.GeoJSON(polygon, {color: 'black', 'fill': false }).addTo(map);

body {
    margin: 0;
}

html, body, #leaflet {
    height: 100%;
}

<!DOCTYPE html>
<html>
  <head>
    <title>Leaflet 1.2.0</title>
    <meta charset="utf-8" />
    <meta name="viewport" content="width=device-width, initial-scale=1.0" />
    <link type="text/css" rel="stylesheet" href="//unpkg.com/leaflet@1.2.0/dist/leaflet.css" />
  </head>
  <body>
    <div id="leaflet"></div>
    <script type="application/javascript" src="//unpkg.com/leaflet@1.2.0/dist/leaflet.js"></script>
    <script type="application/javascript" src="//npmcdn.com/@turf/turf/turf.min.js"></script>
</script>
  </body>
</html>

【讨论】：

但是凸包错过了内部点。如果您考虑 L 形角.. 将选择前两个点如果您看到此图像 raw.githubusercontent.com/AndriiHeonia/hull/master/readme-imgs/… 我不想要那些空白。边界必须沿着边缘延伸。