contextily | Context geo-tiles in Python | Map library

I'm doing a jointplot with a basemap, the problem is that when I add the basemap the main plot doesn't have the same size of the marginal plots. I've tried with different parameters without luck. Does anyone have an idea?

Using the matplotlib backend, is it possible to add a tile basemap such as OSM to a GeoViewsplot, e.g. by somehow calling contextily? Using the Bokeh backend, this is done via gv.tile_sources and then adding it to an overlay but is there a similar function for the mpl backend?

Adding a reproducible example assuming one is switching between backends, and using neighbourhood-level polygon gdfs in EPSG:4326.

What made me initially think adding a basemap was not possible was (1) not defining the WMTS zoom level (causing undecipherable pixelated text to be plotted instead of features), and, after reading James' answer, (2) adding the tiles layer to the layout last, not first, which caused tiles to cover the polygons layers (not an issue on the bokeh backend, but with matplotlib apparently it does matter).

I have a dataframe of individual addresses with lat/long that I converted to a GeoPandas dataframe. I set the projection to the Contextily EPSG:3857 projection.

I'm trying to plot a set of points on a basemap using contextily in python. I have the following points object generated by shapely.geometry, long first and then lat.

0 POINT (-122.44040676 37.7452892821)

1 POINT (-122.432337016 37.7814273807)

2 POINT (-122.409302275 37.7923551768)

3 POINT (-122.446989285 37.7963790178)

4 POINT (-122.412922159 37.8076612624)

I then converted it into a geopandas data and used to following to plot it.

Using Python 3.8, GeoPandas, and Contextily, I am plotting a lot of maps in different areas at different zoom levels (looping through list of points in a GeoDataFrame). Different zoom levels work for different global areas.

What is the best way to set the zoom to the max allowable, capturing the points that I am plotting? In the code below, I don't know how to:

1. Return max zoom allowable (I just used 13 in second plot because that was in error message)
2. Change the extents in the second plot-- the extents of the second plot are the same as the first; I would have expected that changing the zoom level would change the extents

I think I have correctly set the EPSG, which was the gist of the answer to a similar question here.

Code:

I have a point geometry of US locations contained in a GeoDataFrame. I want to plot this as a scatterplot over the US map. My code is:

I am trying to add a real map behind my plot in order to show the locations of my datapoints.

But I don't know how to add the real map...

I've tried several ways, like put another png behind the plot.

But it just not right.

Here is what I got.

I'm trying to add a contextily basemap to a Matplotlib figure containing a GeoPandas data frame. When I just plot the data frame using df.plot the map extent is calculated correctly.

However, when I try adding a contextily basemap the map extent (and zoom level) is calculated wrongly and the following warning is shown:

Please forgive me as this is my first go at a Python 'Project'.

A quick overview: I am trying to produce maps (MatPlotLib figures) by iterating through a GeoPandas dataframe of a GeoJSON file containing the boundaries of active subdivision phases in order to show the progress of construction of the individual lots in each subdivision phase. We normally did this manually in GIS, but I figured I'd take a stab at automating some, or all, of the process. I am also attempting to do this without using ESRI's python functionality and would prefer to keep it that way for stability for future use as ESRI can move around quite a bit.

I am able to iterate through the geodataframe and produce a map (figure) that is zoomed to the extent of the subdivision phase boundary, however, it is clipping the background aerial imagery basemap to the minimum binding box of the subdivision phase that I have been using to set the 'zoom' of the figure.

Example of what I am aiming to produce. This was made in ArcGIS Pro:

What I am able to make in Python:

I am unable to get the basemap to be a constant size that fills a standard landscape letter page whilst still being correctly zoomed in to the extent of the subdivision phase. My output jpeg is sized correctly, but the aerial imagery basemap is continuously croppped to the extent of the subdivision phase boundary leaving large borders around the figure.