from geopandas import GeoDataFrame, GeoSeries, read_file, datasets
import numpy as np
import pandas as pd
from shapely.geometry import Point
class CRS:
def setup(self):
nybb = read_file(datasets.get_path('nybb'))
self.long_nybb = GeoDataFrame(pd.concat(10 * [nybb]),
crs=nybb.crs)
num_points = 20000
longitudes = np.random.rand(num_points) - 120
latitudes = np.random.rand(num_points) + 38
self.point_df = GeoSeries([Point(x, y) for (x, y)
in zip(longitudes, latitudes)])
self.point_df.crs = {"init": "epsg:4326"}
def time_transform_wgs84(self):
self.long_nybb.to_crs({"init": "epsg:4326"})
def time_transform_many_points(self):
self.point_df.to_crs({"init": "epsg:32610"})