title: Make geopandas DataFrame plot on top a raster plot fit correctly
body: <p>I am trying to plot the Pumps.shp data on top of the OSMap.tif file from <a href="http://blog.rtwilson.com/john-snows-cholera-data-in-more-formats/" rel="nofollow noreferrer">this website</a> on the same figure. </p>
  
  <p>I tried using rasterio.plot() and geopandas.plot() methods, with matplotlibs subplots. </p>
  
  <p>The problem is that the plots don't match, the raster file gets plotted in the range(0,1000) for both axis and the shp gets plotted in the actual coordinates range(around 50000 on the x axis and around). </p>
  
  <p>The crs are equal in both objects and the coordinates are in the same range. Why is this? What am I doing wrong?</p>
  
  <p>Here is my code</p>
  
  <pre><code>   import rasterio as rast
     import rasterio.plot as rsplot
     import geopandas as gpd
     src=rast.open("OSMap.tif")
     data=gpd.read_file("Pumps.shp")
     fig,ax=plt.subplots()
     rsplot.show(src,ax=ax)
     data.plot(ax=ax)
     plt.show()
  </code></pre>
  
  <p><strong>This is the result of calling src.bounds:</strong></p>
  
  <p>BoundingBox(left=528765.0, bottom=180466.0, right=529934.0, top=181519.0)</p>
  
  <p><strong>This is the result of data.bounds</strong></p>
  
  <p>(528765.0, 180466.0, 529934.0, 181519.0)</p>
  
  <p><strong>This is crs of both:</strong></p>
  
  <p>CRS({'lon_0': -2, 'y_0': -100000, 'k': 0.9996012717, 'lat_0': 49, 'proj': 'tmerc', 'wktext': True, 'datum': 'OSGB36', 'no_defs': True, 'x_0': 400000, 'units': 'm'})  </p>