はじめに

久保研介研究会のサブゼミの資料です！限定公開です。（2020.08. 公開しました
。）RでのWeb スクレイピング入門の続篇になります。前回の資料で最後に出力したgas_naha.csvを使って、分析を進めていきます！

今日のゴール

ざっくりいうと、この2つです。

「位置情報データって役に立ちそう」と感じる。
ジオコーディングで重要な操作をRで実行できるようになる！

とりあえず、ジオコーディングやってみる．

緯度経度へ変換

最初なので気合を入れていきたいところですが、こちらは専用のサイトにファイルをアップロードするだけで完成します笑。今回は、東大の**ジオコーディングサイト**を使ってみます。手順は下の通りです。

住所を含むカラム番号を入力
- 今回は2行目ですね。csvファイルを開いて確認してみてください。
入力ファイルと出力ファイルの漢字コード
- gas_naha.csvを出力した際に、CP932という意味わからんエンコーディングを指定したので、今回はここをシフトJISと指定しています。
- 普通のcsvファイルでしたら、なにも指定しなくても大丈夫だと思います。もし文字化けしたら、エンコーディングを変えてみてください。
変換したいファイル名
- その名の通りです。

**送信**を押すと、ファイルが返ってくると思います。返ってきたファイルを**gas_naha_new.csv**という名前にしておきましょう。

データを綺麗にする

ジオコーディングしたデータをR studioに取り込みましょう。この時、fileEncodingを指定することに注意してください。Windowsの人は、fileEncoding = utf8とすると、文字化けがなくなるかもです。

gas_naha <- read.csv("gas_naha_new.csv", fileEncoding = "shift-jis")

おなじみのtidyverseを使って、データを綺麗にしていきましょう。selectで必要な変数だけを指定して、renameで変数の名前をわかりやすくしています。

library(tidyverse)
gas_naha <- gas_naha %>% 
  select(store, address,company, fX, fY) %>% # 必要な変数だけ取得
  rename(longitude = fX, latitude = fY) # 変数の名前を変更

おまけ① 地図上にプロットしてみる

早いですが、1つ目のおまけです。前にもやりましたが、地図にガソリンスタンドをプロットしてみます。**leaflet**というライブラリを使います。
下のコードを実行してみてください。mapとうって、世界地図が出てくることを確認してください。

## leafletというライブラリを使用します。 
library(leaflet)

## 地図をRstudio上に用意します。
map <- leaflet(gas_naha) %>% 
  addTiles()

この世界地図上に、沖縄県那覇市のガソリンスタンドをプロットしてみます。addCirclesの代わりにaddMarkersでも良いかもです。集合体恐怖症の方は注意してください。

map %>% 
  addCircles(lng = ~ longitude, lat = ~latitude)

やっぱり、国道沿いに多いですね〜。

おわりに

次回に続く…。今回作ったデータにメッシュ統計データを組み合わせて分析してみます。

地理情報講座リスト