code: def metadata(sceneid, pmin=2, pmax=98, **kwargs):
      
      scene_params = _landsat_parse_scene_id(sceneid)
      meta_data = _landsat_get_mtl(sceneid).get("L1_METADATA_FILE")
      path_prefix = "{}/{}".format(LANDSAT_BUCKET, scene_params["key"])
  
      info = {"sceneid": sceneid}
  
      _stats_worker = partial(
          _landsat_stats,
          address_prefix=path_prefix,
          metadata=meta_data,
          overview_level=1,
          percentiles=(pmin, pmax),
          **kwargs
      )
  
      with futures.ThreadPoolExecutor(max_workers=MAX_THREADS) as executor:
          responses = list(executor.map(_stats_worker, LANDSAT_BANDS))
  
      info["bounds"] = [
          r["bounds"] for b, r in zip(LANDSAT_BANDS, responses) if b == "8"
      ][0]
  
      info["statistics"] = {
          b: v
          for b, d in zip(LANDSAT_BANDS, responses)
          for k, v in d["statistics"].items()
      }
      return info
docs: Retrieve image bounds and band statistics.
  
      Attributes
      ----------
      sceneid : str
          Landsat sceneid. For scenes after May 2017,
          sceneid have to be LANDSAT_PRODUCT_ID.
      pmin : int, optional, (default: 2)
          Histogram minimum cut.
      pmax : int, optional, (default: 98)
          Histogram maximum cut.
      kwargs : optional
          These are passed to 'rio_tiler.landsat8._landsat_stats'
          e.g: histogram_bins=20, dst_crs='epsg:4326'
  
      Returns
      -------
      out : dict
          Dictionary with bounds and bands statistics.