Adding a base map to census geometry in python


I am trying to map census data (downloaded with geometry from NHGIS) in python using matplotlib and adding a basemap using contextily.

Basemaps are in crs 3587, but when I convert the NHGIS geometry from its native format (crs 4257?) to 3587 in order to add the basemap, the geometry column turns into “inf inf”.

Any ideas?

NHGIS shapefiles use Esri’s USA Contiguous Albers Equal Area Conic projected coordinate system. The files come with a “.prj” file that contains the exact specification for this system. I don’t know what the CRS code is. This site suggests to me that the code is 102003.