code: def get_res_stats(ds_list, t_srs=None):
      
      if t_srs is None:
          t_srs = get_ds_srs(ds_list[0]) 
      res = np.array([get_res(ds, t_srs=t_srs) for ds in ds_list])
      #Check that all projections are identical
      #gt_array = np.array([ds.GetGeoTransform() for ds in args])
      #xres = gt_array[:,1]
      #yres = -gt_array[:,5]
      #if xres == yres:
      #res = np.concatenate((xres, yres))
      min = np.min(res)
      max = np.max(res)
      mean = np.mean(res)
      med = np.median(res)
      return (min, max, mean, med)
docs: Return resolution stats for an input dataset list