from cartography.proj import SpatialReference from cartography.proj.transform import transform_points from cartography.proj import ProjTransform #datum: #lat_1: Standard Parallel 1 #lat_2: Standard Parallel 2 #lat_0: Latitude of Origin #lon_0: Central Meridian #x_0: False Easting #y_0: False Northing #ellps: Ellipsoid #units: Meters custom = SpatialReference("""+proj=aea +datum=NAD83 +lat_1=24.0 +lat_2=31.5 +lat_0=24.00 +lon_0=-84.00 +x_0=400000 +y_0=0 +ellps=wgs84 +units=m """) wgs84 = SpatialReference(epsg=4326) t = ProjTransform(custom, wgs84) t.transform([(375365, 686613)]) print t