How can I apply Robinson projection in a GeoPandas plot?
import geopandas as gpd
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
world.plot(column='pop_est')
Robinson projection in GeoPandas
coordinate systemepsggeopandasplotpython
Related Solutions
if you are looking for description of Robinson here it is: http://en.wikipedia.org/wiki/Robinson_projection and algorithm in JS you can download from here: https://github.com/afar/robinson_projection and script descrition is here: http://www.flashmap.org/robinson-projection-in-as3-gpl/
You can do this using cartopy.
#!/usr/bin/env python
import cartopy
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
from pylab import imread
im = imread('Robinson-projection.jpg')
ax = plt.axes(projection=ccrs.Robinson())
plt.imshow(im, origin='upper', extent=[-17005833.330525, 17005833.330525, -8622512.772008, 8622512.772008], interpolation='nearest')
ax.coastlines(resolution='110m', color='yellow', linewidth=1, alpha=0.7)
plt.plot(-122.4194155, 37.7749295, marker='o', color='red', transform=ccrs.Geodetic())
plt.show()
To test I used a map in the Robinson projection downloaded from Wikipedia.
http://en.wikipedia.org/wiki/List_of_cartographers#mediaviewer/File:Robinson-projection.jpg
The result is shown below.
To convert coordinates between lnglat and WGS84 using pyproj:
#!/usr/bin/env python
import pyproj
crs_from = pyproj.Proj(init='EPSG:4326')
crs_to = pyproj.Proj('+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs')
x, y = pyproj.transform(crs_from, crs_to, -122.4194155, 37.7749295)
Best Answer
Use the
cartopylibrary to bring in the Robinson details.More details on the GeoPandas Docs.