OSMとturf.jsで、館林の最寄りバス停発見データを作ってみる (最寄りポリゴン生成編)

const fs = require('fs');
const voronoi = require('@turf/voronoi');
const buffer = require('@turf/buffer');
const booleanPointInPolygon = require('@turf/boolean-point-in-polygon').default;
const intersect = require('@turf/intersect').default;

const bus = JSON.parse(fs.readFileSync('./bus.geojson'));
const train = JSON.parse(fs.readFileSync('./train.geojson'));
const boundary = JSON.parse(fs.readFileSync('./boundary.geojson'));
// 市町境から2kmのバッファポリゴン生成
const bufferBoundary = buffer(boundary.features[0], 2);

// バス停データを名寄せ処理
bus.features = bus.features.reduce((prev, curr, index, arr) => {
  const name = curr.properties.name;
  if (!prev[name]) prev[name] = [];
  prev[name].push(curr);
  if (index !== arr.length - 1) return prev;
  return Object.keys(prev).map((key) => {
    return prev[key].reduce((prev, curr, index, arr) => {
      if (!prev) {
        prev = curr;
        prev.num = 1;
      } else {
        prev.geometry.coordinates[0] += curr.geometry.coordinates[0];
        prev.geometry.coordinates[1] += curr.geometry.coordinates[1];
        prev.num++;
      }
      if (index !== arr.length - 1) return prev;
      prev.geometry.coordinates[0] /= prev.num;
      prev.geometry.coordinates[1] /= prev.num;
      delete prev.num;
      return prev;
    }, null);
  });
}, {}).filter((item) => { // バッファポリゴンでフィルタ
  return booleanPointInPolygon(item, bufferBoundary);
});

// 鉄道駅データをフィルタしたうえで、名寄せ済みバス停データにマージ
train.features.filter((item) => {
  return booleanPointInPolygon(item, bufferBoundary);
}).map((item) => {
  bus.features.push(item);
});

// ボロノイ図生成
const voronoiObj = voronoi(bus);
// ボロノイ図をバッファポリゴンでくり抜き
voronoiObj.features = voronoiObj.features.map((item) => {
  return intersect(item, bufferBoundary);
}).filter((item) => {
  return item ? true : false;
});

// くり抜かれたボロノイ図を voronoi.geojsonの名前でセーブ
fs.writeFileSync('./voronoi.geojson', JSON.stringify(voronoiObj));