title: is it possible to make a join between a geotiff image with a geopandas dataframe?
body: <p>I have a shapfile and a geotiff image.</p>
  
  <pre><code>import geopandas as gpd
  import rasterio
  from rasterio.plot import show
  
  df = gpd.read_file('myShape.shp')
  fileI = 'myFile.tiff'
  
  data = rasterio.open(fileI)
  show((data), cmap='terrain', ax=ax)
  p1 = df.geometry.boundary.plot(color=None,edgecolor='red',linewidth = 2,ax=ax)
  ax.axis('off');
  </code></pre>
  
  <p><a href="https://i.stack.imgur.com/ocX70.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/ocX70.png" alt="enter image description here"></a></p>
  
  <p>I would like to assign the values of <code>data</code> only to the red region defined by <code>df</code>. Is it possible to crop the image around that region or to make a join between them?</p>