Ecosyste.ms: Awesome

An open API service indexing awesome lists of open source software.

Awesome Lists | Featured Topics | Projects

https://github.com/aidatorajiro/takasa

OpenStreetMap 高低差を考えた最短（？）経路問題ソルバー
https://github.com/aidatorajiro/takasa

numpy openstreetmap scipy shortest-path

Last synced: 5 days ago
JSON representation

OpenStreetMap 高低差を考えた最短（？）経路問題ソルバー

Awesome Lists containing this project

README 

        
# OpenStreetMap 高低差を考えた最短（？）経路問題ソルバー

![takasa wo kangaeru](takasa.png)



OpenStreetMapの地図データと国土地理院の標高データを組み合わせて、高低差が出来るだけ小さくなり、かつ経路も出来るだけ短くなるような経路をかんがえます。



## 実行方法



まず国土地理院のサイトにいき、探索範囲の「数値標高モデル」をダウンロードしましょう。（要会員登録）



zipファイルを全て解凍し、ダウンロードした大量のxmlファイルを実行ディレクトリの`DL`下に配置します。



次にPythonファイルを作成し、以下のように目的地などの情報を設定して、solve関数を呼びます。



経路の探索は指定された緯度・経度からなる矩形の範囲内で行われます。緯度・経度は小数点表記で書いてください。



```python

import solve

import elevation

import mapdata



pathbound_small = [探索範囲の最小緯度, 探索範囲の最小経度]

pathbound_large = [探索範囲の最大緯度, 探索範囲の最大経度]



pathstart_nodeid = 出発地のOpenStreetMap上のNode ID

pathend_nodeid = 目的地のOpenStreetMap上のNode ID



mapdata.generate_mapdata(pathbound_small, pathbound_large)

elevation.generate_elevation()

solve.solve(pathstart_nodeid, pathend_nodeid)

```



Node IDはなどで検索しましょう。右クリックで「地物を検索(Query Features)」すると付近の道が出現します。道のリンクをクリックすると、その道を構成するノードが現れます。



実行後、経路のノードID一覧`nodeids.results`と、



それをOverpass APIが処理できるようにした`overpass_query.results`が生成されます。に入力すると、マップと照らし合わせることができます。



また、範囲内の標高データ・マップデータをPythonが扱いやすいように処理したファイルである`elevation.shel`と`mapdata`も生成されます。これらは範囲(`pathbound_small`と`pathbound_large`)が変わらない限り、使い回すことができます。



## 仕様



てきとうに考えました。



1. scipyの`shortest_path`関数を使って解きます。

1. 探索対象の道は`highway`パラメータの存在する全ての道です。つまりフェリーや電車の路線などでない、歩行者や自動車が通るような全ての道にマッチします。

1. 一方通行などの区別はないです。つまり、`shortest_path`関数に渡す行列は対称行列です。

1. 道の長さと高低差をスコア化し、そのスコアの合計値が最小になるような経路を選びます。スコアの計算方法にはさまざまなものを実装しておりますますので、スコアの詳細は`solve.py`の100行目あたりを読もう！