未来のタクシー需要を先読みしよう!
mini_forest
・State Plane Coordinate Systemが初見だったので、緯度経度変換を試みましたのでメモです。・とりあえずfolium上に点を打ちたかったのでメモです。
import pandas as pd pd.options.display.max_columns = 100 pd.options.display.max_rows = 999 pd.options.display.float_format = '{:.6f}'.format import warnings warnings.filterwarnings('ignore') from pyproj import Proj, transform import folium
from google.colab import drive drive.mount('/content/drive')
Drive already mounted at /content/drive; to attempt to forcibly remount, call drive.mount("/content/drive", force_remount=True).
path = "/content/drive/MyDrive/tmp/7_probspace/3_taxi/" zone = pd.read_csv(path+'data/taxi_zones.csv')
# 緯度経度に変換 def sp2latlong(row): inProj = Proj(init='epsg:2263', preserve_units=True) outProj = Proj(init='epsg:4326') return transform(inProj, outProj, row["X"], row["Y"]) zone[["long","lat"]] = zone.apply(lambda x:sp2latlong(x), axis=1, result_type='expand') zone
map = folium.Map(location=[zone["lat"].mean(), zone["long"].mean()], tiles='cartodbpositron',zoom_start=11) for i, d in zone.iterrows(): folium.Marker(location=[d["lat"], d["long"]]).add_to(map) map
mini_forestさま
foliumやprojのような便利なツールがあるとは知りませんでした。ありがとうございました。
とんでもないです!!せっかくの位置情報なので何か有用なことに使えないかなと思い調べてみた内容の共有でした!