docs: def allstats(self, approximate=False):
          """
          Compute some basic raster statistics
  
          Parameters
          ----------
          approximate: bool
              approximate statistics from overviews or a subset of all tiles?
  
          Returns
          -------
          list of dicts
              a list with a dictionary of statistics for each band. Keys: `min`, `max`, `mean`, `sdev`.
              See :osgeo:meth:`gdal.Band.ComputeStatistics`.
          """
          statcollect = []
          for x in self.layers():
              try:
                  stats = x.ComputeStatistics(approximate)
              except RuntimeError:
                  stats = None
              stats = dict(zip(['min', 'max', 'mean', 'sdev'], stats))
              statcollect.append(stats)
          return statcollect