使用 query-overpass 和 turf.js 将 geojson 多边形转换为带节点的点

Question

我使用节点模块 "query-overpass" 进行查询以从 openstreetmaps 获取农场商店。我想将所有多边形转换为此脚本中的点。我使用 turf.js 来获取这些多边形的质心，但我无法以永久方式更改对象。到目前为止，这是我的代码：

const query_overpass = require("query-overpass");
const turf = require ("turf");
const fs = require("fs")
let test
let filename = "data/test.js"
let bbox = "48.91821286473131,8.309097290039062,49.0610446187357,8.520584106445312";
console.log('starting query for ' +filename)
console.log('bbox: ' +bbox)
let query = ` 
  [out:json][timeout:250];
  // gather results
  (
  // query part for: “vending=milk”
  node["vending"="milk"](${bbox});
  way["vending"="milk"](${bbox});
  relation["vending"="milk"](${bbox});
  // query part for: “shop=farm”
  node["shop"="farm"](${bbox});
  way["shop"="farm"](${bbox});
  relation["shop"="farm"](${bbox});
  // query part for: “vending=food”
  node["vending"="food"](${bbox});
  way["vending"="food"](${bbox});
  relation["vending"="food"](${bbox});
  );
  // print results
  out body;
  >;
  out skel qt;
`;

// query overpass, write result to file
query_overpass(query, (error, data)  => {
    data = JSON.stringify(data , null, 1)
    console.log(data)
    test = JSON.parse(data)


//create centroids for every polyon and save them as a point
    for (var i = 0; i < test.features.length; i++) { 
      console.log("Log: " +test.features[i].geometry.type)
      console.log("Log: " +test.features[i].properties.name)
      if (test.features[i].geometry.type === "Polygon"){
        console.log("polygon detected")
        var centroid = turf.centroid(test.features[i]);
            var lon = centroid.geometry.coordinates[0];
            var lat = centroid.geometry.coordinates[1];
            console.log(" lon: " +lon +" lat: " +lat)

            test.features[i].geometry.type = 'Point'
            //delete Polygon structure and insert centroids as new points here
            console.log("polygon deleted and changed to point")
      }

  }
  console.log(test)
   fs.writeFile(filename, `var file = ${test};` , ["utf-8"], (error, data) => {if (error) {console.log(error)}})
 }, {flatProperties: true}
)

似乎我可以在for循环中改变一些东西，但是当稍后保存数据时它们不会出现。这基本上是一个如何正确编辑 json 对象的问题，但我不明白为什么这在这里根本不起作用。

所以基本上有两个问题：

为什么我不能覆盖上面示例中的 geometry.type？
如何删除旧多边形并向要素添加新点？

感谢您的帮助。

Answer 1

这很复杂...你为什么不让 Overpass API 来做这个工作并使用 out center; 而不是 out body;>;out skel qt; 到 return 的中心点所有节点、方式和关系。你可以先用overpass-turbo.eu试试这个

使用 query-overpass 和 turf.js 将 geojson 多边形转换为带节点的点

Convert geojson polygons into points with node using query-overpass and turf.js

node.js

openstreetmap

geojson

overpass-api

turfjs