code: def get_overlay_bounds(cls, overlay):
          
          if all(el.bounds==overlay.get(0).bounds for el in overlay):
              return overlay.get(0).bounds
          else:
              raise ValueError("Extents across the overlay are inconsistent")
docs: Returns the extents if all the elements of an overlay agree on
          a consistent extents, otherwise raises an exception.