Change CRS of Pandana network

gistfile1.txt

import pyproj
import pandas as pd

s = pd.HDFStore('osm_bayarea.h5', "r")
s2 = pd.HDFStore('osm_bayarea4326.h5', "w")
s2["edges"] = s["edges"]

po = s["nodes"]

print po.describe()

x_name = "x"
y_name = "y"

p1 = pyproj.Proj('+init=epsg:%d' % 3740)
p2 = pyproj.Proj('+init=epsg:%d' % 4326)
x2, y2 = pyproj.transform(p1, p2,
                          po[x_name].values,
                          po[y_name].values)
po[x_name], po[y_name] = x2, y2

print po.describe()

s2["nodes"] = po