如何有效更改Geopandas数据框中的几何crs？

Question

我有一个shapefile，我正在将其转换为geopandas数据框。此shapefile的crs是epsg：4326，我想将其转换为epsg：3857。

import geopandas as gp

file_path = 'zip_poly.shp'
sf = gp.read_file(file_path)
sf = sf[['ZIP_CODE','PO_NAME','STATE','geometry']]
sf = sf.to_crs(epsg=3857)

您可以找到shapefile here。

shapefile包含约33,000个邮政编码，转换时间太长。有没有更好的方法来快速转换？