docs: def plot(self, figsize=None, ax=None, **kwargs):
          '''
          geo.plot()
  
          Returns plot of raster data
          '''
          if ax is None:
              fig, ax = plt.subplots(figsize=figsize)
          ax.set_aspect('equal')
          ax.matshow(self.raster, **kwargs)
          plt.draw()
          return ax