### OSMnxとQGISで ### AEDエリアマップを ### 作ってみたよ Yasuto FURUKAWA --- ## あんた誰 - GISのたのしい匍匐前進業 - Python真面目にやり(直し)始めたのここ１年くらいのﾍﾀﾚ🐣 - 所属 - [MIERUNE](https://mierune.co.jp) - [Code for Japan](https://www.code4japan.org/) - [OSGeoJP](https://www.osgeo.jp) - Qiitaへの技術ポエムの一例 - [GeoPandasは100万行の夢をみるか](https://qiita.com/Yfuruchin/items/6b6f5b6e787e83ea2709) - [上司に「ひなたGISみたいなの午後までに作って」と言われた時のために](https://qiita.com/Yfuruchin/items/9f466d4bedcb5bda3038) --- ### UDC 2014 アプリ部門金賞  https://www.slideshare.net/Mayumit/ss-45428322 --- ### UDCの（サブ）マークつくりました <img src="https://i.imgur.com/jSuX6yY.png" width="60%"> --- ## 今回やったこと <!-- .slide: data-background="https://i.imgur.com/EPr07Ft.png"--> --- ### ゴール [日本AED財団](https://aed-zaidan.jp/knowledge/index.html)によると、国内では一日に200人もの人が心臓突然死で亡くなっています。このような悲劇をすこしでも軽減するためAED（自動体外式除細動器）の設置が各地ですすめられていまが、AEDがあったとしても、どれくらいの範囲の人々を助けれるのでしょうか？ ここでは、ケーススタディとして東京都大田区のAEDの位置情報をもとに駆け足（300m/min）で往復3分(900m)以内の経路をOSMnxというPythonパッケージで計算し、QGISで可視化しました。 <!-- .slide: data-background="https://i.imgur.com/EPr07Ft.png"--> --- ### 要するに 大田区AEDからの経路計算をして 可視化してみました --- ## いいたいこと - Pythonでざっくり経路検索できるよ - カッコイ良い図面はQGISにおまかせ - データやツール０円でこんな絵ができます - 大田区のAED情報… --- ## ながれ - 経路検索って - OSMnxとは - OSMnxデモ - 実際には - QGISでの表現 - QGISデモ - 大田区のAED情報… - issue --- ## 経路検索アツい（個人的に） - バッファ距離…空飛ばないし - 実際に歩く(走る/逆立ち)で移動する距離で考えましょうや --- ### バッファvs経路検索 （バッファのほうがOver-Estimate）  --- ## OSMnxとは - [OpenStreetMap](https://www.openstreetmap.org)から道路網を - 取得 - 構築 - 分析 - 視覚化 してくれるパッケージ --- ## OSMnx参考情報 - [翻訳解説](https://githubja.com/gboeing/osmnx) - [Example](https://github.com/gboeing/osmnx-examples/tree/master/notebooks) --- ## 注意 - OpenStreetMapの[OverpassAPI](https://wiki.openstreetmap.org/wiki/JA:Overpass_API)たたくのでネット接続必須 - 1日あたり10000クエリまたは5GBのデータくらいが上限 - 優しくね！ --- ## OSMnxインストール ``` brew install spatialindex pip3 install osmnx ``` Pandas 🐼 GeoPandas 🌏🐼 など目的によってinstall --- ## OSMnxの挙動 ３−５秒/ポイントくらい？  --- ## 実際にやったこと 1. 大田区AED情報から緯度経度取得 2. 緯度経度地点のリストつくる 3. 距離のリストつくる 4. ２と３からloopでOSMnxに投げる 5. 結果を **GeoDataFrame** に統合 6. QGISでお化粧 --- ## QGISでの表現 [地理院地図ダーク化](https://maps.gsi.go.jp/development/ichiran.html)＋色合い工夫  --- ## 大田区のAED情報… --- ## pdfェ…  --- ## 頑張ってスクレイピング！ 緯度経度っぽいデータある！ http://ota-wagamachi.denen.com/?zoom=18&lat= ***35.57783380665689&lng=139.70937109611376*** &theme_id=22&layer_id=191&type=roadmap&icon_size=32&language=ja --- ## ぐりぐり回して可視化だぜ！ --- ## おかしい！…妙にAED少ない！ --- ## ファイル名🤮  --- ### またか…orz - 整然化されたオープンデータ超大事です - 失望GeoPythonさんをもう増やさないで！ --- ## Issue --- ## Issue - OpenStreetMapのデータ精度に左右されるので過信してはいけない - 特にノードの接続関係とか - QGISの結果とずれてるので要検証 - ガチでやるならキレイなデータもってくる - レスポンスそんな早くない - GRASS, PgRoutingパイセン - Google Colaboratoryでやってみたい！ --- ## まとめ - Pythonで（ざっくりな）経路検索できるよ - カッコイ良い図面はQGISにおまかせ - データやツール０円でこんな絵ができます - 大田区のAED情報orz --- ## 参考：GeoPandas🌏🐼!  http://geopandas.org/ いろんなトライできて素敵だけど日本語の情報少ないので貢献者求む！ --- ## 参考:GeoPandas🌏🐼!  https://2019.foss4g.org/schedule/workshops/ --- ## 参考:GeoPandas🌏🐼!  https://qiita.com/Yfuruchin/items/6b6f5b6e787e83ea2709 --- ## 参考リンク - XYのリストからネットワークを計算して１つのGeoDataFrameにまとめてSHPにしてくれる雑い[サンプルコード](https://github.com/Yfuruchin/pyconjp2019_geopyrhon/blob/master/points2networkshp.ipynb) - [OSMnx](https://osmnx.readthedocs.io/en/stable/) - [QGIS](https://qgis.org/) - [大田区AED設置場所](http://www.city.ota.tokyo.jp/seikatsu/hoken/iryokikan_hoka/aed/setting_place/index.html)