Question

我在geoJSON文件中有一个用多边形描述的城市列表。

我想在多边形的内部获取一个样本。

基本数学说，重心在多边形内部，足以将所有经度和所有纬度加在一起，然后将其除以点数。

Full file to process（可视化在GitHub上可用）

{
  "type": "FeatureCollection",
  "features": [
  {
    "type": "Feature",
    "geometry": {
      "type": "Polygon",
      "coordinates": [[[2.41101, 48.72605], [2.41554, 48.72656], [2.41718, 48.72791], [2.4211, 48.72953], [2.42603, 48.72824], [2.42756, 48.72865], [2.42922, 48.72723], [2.43133, 48.72646], [2.43404, 48.72665], [2.43513, 48.72409], [2.42554, 48.7227], [2.42072, 48.72105], [2.41426, 48.71782], [2.41327, 48.71869], [2.41582, 48.72086], [2.41238, 48.72193], [2.41136, 48.72325], [2.41101, 48.72605]]]
    },
    "properties": {
      "code": "94001",
      "nom": "Ablon-sur-Seine"
    }
  },
  {
    "type": "Feature",
    "geometry": {
      "type": "Polygon",
      "coordinates": [[[2.41959, 48.81691], [2.4159, 48.81633], [2.40936, 48.81667], [2.40787, 48.81746
    },
    "properties": {
      "code": "94018",
      "nom": "Charenton-le-Pont"
    }
  },
  ...
  ]
}

我已经有一个计算多边形顶点长度的命令。

$ curl -s https://raw.githubusercontent.com/gregoiredavid/france-geojson/master/departements/94-val-de-marne/communes-94-val-de-marne.geojson \
> | jq '.features[0].geometry.coordinates[0][][0]' \
> | jq -s 'add/length'
2.4206944444444445

请参见https://unix.stackexchange.com/questions/13731/

如何使用jq和简单的bash命令，如何计算经度和和纬度的总和，并在另一个geoJSON文件的properties字段中重新注入重心？

谢谢。

Answer 1

因此，如果我理解正确，那么您正在尝试获取第一组坐标的平均值，然后更新属性以存储结果。

.features[] |= (
    (.geometry.coordinates[0] | length as $len | reduce .[] as [$x, $y] ([0,0];
        [.[0] + $x, .[1] + $y]
    ) | map(. / $len)) as $barrycenter |
    .properties.barycenter = $barrycenter
)

Answer 2

警告

多边形的“重心”通常与根据其顶点的x和y坐标的平均值定义的点不同。参见例如https://math.stackexchange.com/questions/3177/why-doesnt-a-simple-mean-give-the-position-of-a-centroid-in-a-polygon

纬度和经度的平均值

这是一个jq过滤器，它将只需一次调用jq且没有冗余，就可以计算每个“特征”的纬度和经度的平均值：

.features[].geometry.coordinates[0]
| [ [.[][0]], [.[][1]] ]
| map(add/length)

使用-c命令行选项，将生成一个数组流，每个“功能”一个数组。流开始：

[2.4206944444444445,48.724651111111115]
[2.407614,48.82250133333333]
...

当然还有其他选择，但是请注意，不需要使用字符串插值法进行分组，因此通常不需要像您的第一个版本那样使用tr。

使用jq计算geoJSON多边形的重心

2 个答案:

警告

纬度和经度的平均值