{"action":"opened","issue":{"url":"https://api.github.com/repos/geopandas/geopandas/issues/858","repository_url":"https://api.github.com/repos/geopandas/geopandas","labels_url":"https://api.github.com/repos/geopandas/geopandas/issues/858/labels{/name}","comments_url":"https://api.github.com/repos/geopandas/geopandas/issues/858/comments","events_url":"https://api.github.com/repos/geopandas/geopandas/issues/858/events","html_url":"https://github.com/geopandas/geopandas/issues/858","id":376634440,"node_id":"MDU6SXNzdWUzNzY2MzQ0NDA=","number":858,"title":"Plot a kde with epsg 3857 projection","user":{"login":"joanbn","id":12449412,"node_id":"MDQ6VXNlcjEyNDQ5NDEy","avatar_url":"https://avatars1.githubusercontent.com/u/12449412?v=4","gravatar_id":"","url":"https://api.github.com/users/joanbn","html_url":"https://github.com/joanbn","followers_url":"https://api.github.com/users/joanbn/followers","following_url":"https://api.github.com/users/joanbn/following{/other_user}","gists_url":"https://api.github.com/users/joanbn/gists{/gist_id}","starred_url":"https://api.github.com/users/joanbn/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/joanbn/subscriptions","organizations_url":"https://api.github.com/users/joanbn/orgs","repos_url":"https://api.github.com/users/joanbn/repos","events_url":"https://api.github.com/users/joanbn/events{/privacy}","received_events_url":"https://api.github.com/users/joanbn/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":0,"created_at":"2018-11-02T01:31:58Z","updated_at":"2018-11-02T01:31:58Z","closed_at":null,"author_association":"NONE","body":"I can't project KDE and base map tiles together, since map tiles need 3857 projection and I cannot find it in Geopandas projection.\r\n\r\nThis is a code example :\r\n\r\ncrs = {'init': 'epsg:4326'}\r\ngeo_df = GeoDataFrame(df, crs=crs, geometry=geometry)\r\ngeo_df = geo_df.to_crs(epsg=3857)\r\n\r\ndef add_basemap(ax, zoom, url='http://tile.stamen.com/terrain/tileZ/tileX/tileY.png'):\r\n    xmin, xmax, ymin, ymax = ax.axis()\r\n    basemap, extent = ctx.bounds2img(xmin, ymin, xmax, ymax, zoom=zoom, url=url)\r\n    ax.imshow(basemap, extent=extent, interpolation='bilinear')\r\n    # restore original x/y limits\r\n    ax.axis((xmin, xmax, ymin, ymax))\r\n\r\n\r\nax = plt.subplot(121,  projection=ccrs.epsg('3857'))\r\ngplt.kdeplot(geo_df, projection=gcrs.AlbersEqualArea(), shade=True, shade_lowest=False,ax=ax)\r\nadd_basemap(ax, zoom=12)    \r\n\r\nI can't find the geoplot projection that corresponds to EPSG:3857\r\n\r\nHow it could be solved?\r\n"}}
{"actor":{"display_login":"joanbn"}}