Adding a base map to census geometry in python

NHGIS
1

Hello!

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?

2

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.