1 files changed, 244 insertions, 0 deletions
diff --git a/scripts/update-osm-parks.py b/scripts/update-osm-parks.py
new file mode 100644
index 0000000..fffb2f3
--- /dev/null
+++ b/scripts/update-osm-parks.py
@@ -0,0 +1,244 @@
+#!/usr/bin/env python3
+"""小平市の公園ポリゴンをOpenStreetMapから取得して更新する。
+
+使い方:
+    py -X utf8 scripts/update-osm-parks.py
+
+出力:
+    src/data/kodaira-parks-osm.json
+    src/data/kodaira-boundary.json（小平市の境界線）
+
+OpenStreetMapのデータは ODbL 1.0。地図表示時は
+「© OpenStreetMap contributors」の帰属表示が必要。
+"""
+import json
+import os
+import sys
+import time
+import unicodedata
+import urllib.request
+
+CITY = "小平市"
+OURS = "src/data/kodaira-parks.json"
+OUT = "src/data/kodaira-parks-osm.json"
+BOUNDARY_OUT = "src/data/kodaira-boundary.json"
+OVERRIDES = "src/data/park-name-overrides.json"
+
+# OSMの位置補正（メートル）。地図側のLeafletペインに付いていた
+# 16pxの余白をCSSで解消したため、現在は0（補正なし）。
+# 実位置との差が残る場合だけ正の値を入れる。
+NORTH_SHIFT_M = 0
+LAT_PER_M = 1 / 111320
+
+ENDPOINTS = [
+    "https://overpass-api.de/api/interpreter",
+    "https://overpass.kumi.systems/api/interpreter",
+    "https://overpass.private.coffee/api/interpreter",
+]
+
+QUERY = f"""[out:json][timeout:90];
+area["name"="{CITY}"]["boundary"="administrative"]->.a;
+(
+  way["leisure"="park"](area.a);
+  relation["leisure"="park"](area.a);
+);
+out geom;
+"""
+
+BOUNDARY_QUERY = f"""[out:json][timeout:60];
+rel["name"="{CITY}"]["boundary"="administrative"];
+out geom;
+"""
+
+# 名称が一致しないが、住所・座標から同一と確認できた公園。
+ALIAS = {
+    "小平中央公園": "中央公園",
+    "じゃぶじゃぶ公園": "天神じゃぶじゃぶ公園",
+    "リサイクルセンター東公園": "リサイクルセンター南公園",
+    "小平市立薬用植物園前公園": "薬用植物園前公園",
+    "花南ポケットパーク": "花南ポケットパーク公園",
+    "花５四季公園": "花五四季公園",
+    "天神町のびのび広場": "天神町のびのび広場公園",
+}
+
+
+def norm(s):
+    s = unicodedata.normalize("NFKC", str(s)).replace(" ", "").replace("　", "")
+    return s.replace("ヶ", "ケ").replace("ツ", "ッ")
+
+
+def fetch(query):
+    last = None
+    for ep in ENDPOINTS:
+        try:
+            req = urllib.request.Request(
+                ep,
+                data=query.encode("utf-8"),
+                headers={"User-Agent": "yasutakeyohei.com park map updater"},
+            )
+            with urllib.request.urlopen(req, timeout=120) as r:
+                data = json.loads(r.read().decode("utf-8"))
+            if data.get("elements"):
+                return data
+            last = "0 elements"
+        except Exception as e:  # noqa: BLE001
+            last = e
+        print(f"  retry: {ep}: {last}", file=sys.stderr)
+        time.sleep(2)
+    raise SystemExit(f"Overpass API failed: {last}")
+
+
+def pt(g):
+    return (round(g["lon"], 6), round(g["lat"] + NORTH_SHIFT_M * LAT_PER_M, 6))
+
+
+def ring_from_ways(ways):
+    rings = []
+    remaining = [list(w) for w in ways if len(w) >= 2]
+    while remaining:
+        ring = remaining.pop(0)
+        while ring[0] != ring[-1] and remaining:
+            found = False
+            for i, w in enumerate(remaining):
+                if w[0] == ring[-1]:
+                    ring.extend(w[1:]); remaining.pop(i); found = True; break
+                if w[-1] == ring[-1]:
+                    ring.extend(list(reversed(w))[1:]); remaining.pop(i); found = True; break
+                if w[0] == ring[0]:
+                    ring = list(reversed(w))[:-1] + ring; remaining.pop(i); found = True; break
+                if w[-1] == ring[0]:
+                    ring = w[:-1] + ring; remaining.pop(i); found = True; break
+            if not found:
+                break
+        if ring[0] != ring[-1]:
+            ring.append(ring[0])
+        if len(ring) >= 4:
+            rings.append(ring)
+    return rings
+
+
+def build_boundary():
+    """小平市の行政境界をGeoJSONにする。"""
+    data = fetch(BOUNDARY_QUERY)
+    if not data.get("elements"):
+        print("  boundary relation not found", file=sys.stderr)
+        return None
+    rel = data["elements"][0]
+    ways = []
+    for m in rel.get("members", []):
+        if m.get("role") == "outer" and m.get("geometry"):
+            ways.append([pt(p) for p in m["geometry"]])
+    rings = ring_from_ways(ways)
+    if not rings:
+        print("  boundary ring not assembled", file=sys.stderr)
+        return None
+    if len(rings) == 1:
+        geom = {"type": "Polygon", "coordinates": [rings[0]]}
+    else:
+        geom = {"type": "MultiPolygon", "coordinates": [[r] for r in rings]}
+    props = {
+        "name": rel.get("tags", {}).get("name", CITY),
+        "boundary": "administrative",
+        "admin_level": rel.get("tags", {}).get("admin_level", ""),
+    }
+    fc = {"type": "FeatureCollection", "features": [{"type": "Feature", "properties": props, "geometry": geom}]}
+    return fc, len(rings), sum(len(r) for r in rings)
+
+
+def main():
+    ours = json.load(open(OURS, encoding="utf-8"))
+    our_names = {norm(r["公園名"]): r["公園名"] for r in ours}
+    alias_norm = {norm(k): v for k, v in ALIAS.items()}
+
+    # OSM ID単位の手動上書き（osm_type/osm_id -> 市の公園名）
+    overrides = {}
+    if os.path.exists(OVERRIDES):
+        raw = json.load(open(OVERRIDES, encoding="utf-8"))
+        for key, entry in raw.items():
+            city = entry.get("city", "") if isinstance(entry, dict) else str(entry)
+            city = str(city).strip()
+            if not city:
+                continue
+            canonical = our_names.get(norm(city))
+            if canonical:
+                overrides[key] = canonical
+            else:
+                print(f"  warning: override city not found: {key} -> {city}", file=sys.stderr)
+
+    def match(name):
+        k = norm(name)
+        if k in our_names:
+            return our_names[k]
+        if k in alias_norm:
+            return alias_norm[k]
+        # 「〇〇児童公園」→「〇〇公園」のような表記差を吸収
+        k2 = k.replace("児童", "")
+        if k2 in our_names:
+            return our_names[k2]
+        return ""
+
+    data = fetch(QUERY)
+    features = []
+    applied = 0
+    for e in data["elements"]:
+        tags = e.get("tags", {})
+        name = tags.get("name", "")
+        key = f"{e['type']}/{e['id']}"
+        our_name = overrides.get(key) or match(name)
+        if key in overrides:
+            applied += 1
+        props = {
+            "name": name,
+            "our_name": our_name,
+            "matched": bool(our_name),
+            "osm_type": e["type"],
+            "osm_id": e["id"],
+        }
+        for k in ["leisure", "addr:city", "addr:street", "addr:housenumber", "opening_hours"]:
+            if k in tags:
+                props[k] = tags[k]
+
+        if e["type"] == "way":
+            ring = [pt(p) for p in e.get("geometry", [])]
+            if len(ring) < 3:
+                continue
+            if ring[0] != ring[-1]:
+                ring.append(ring[0])
+            geom = {"type": "Polygon", "coordinates": [ring]}
+        else:
+            outers, inners = [], []
+            for m in e.get("members", []):
+                pts = [pt(p) for p in m.get("geometry", [])]
+                if len(pts) < 2:
+                    continue
+                (inners if m.get("role") == "inner" else outers).append(pts)
+            outer_rings = ring_from_ways(outers)
+            inner_rings = ring_from_ways(inners)
+            if not outer_rings:
+                continue
+            if len(outer_rings) == 1:
+                geom = {"type": "Polygon", "coordinates": [outer_rings[0]] + inner_rings}
+            else:
+                geom = {"type": "MultiPolygon", "coordinates": [[r] + inner_rings for r in outer_rings]}
+
+        features.append({"type": "Feature", "properties": props, "geometry": geom})
+
+    features.sort(key=lambda f: (f["properties"].get("name", ""), f["properties"]["osm_type"], f["properties"]["osm_id"]))
+    out = {"type": "FeatureCollection", "features": features}
+    json.dump(out, open(OUT, "w", encoding="utf-8"), ensure_ascii=False, separators=(",", ":"))
+
+    matched = sum(1 for f in features if f["properties"]["matched"])
+    named = sum(1 for f in features if f["properties"].get("name"))
+    print(f"features {len(features)} (named {named}, matched {matched}) -> {OUT}")
+    if overrides:
+        print(f"  overrides applied {applied} / {len(overrides)}")
+
+    boundary = build_boundary()
+    if boundary:
+        boundary_data, ring_count, point_count = boundary
+        json.dump(boundary_data, open(BOUNDARY_OUT, "w", encoding="utf-8"), ensure_ascii=False, separators=(",", ":"))
+        print(f"boundary rings {ring_count}, points {point_count} -> {BOUNDARY_OUT}")
+
+
+if __name__ == "__main__":
+    main()