code: def all_axis_bin_centers(self, axis):
          
          # Arcane hack that seems to work, at least in 3d... hope
          axis = self.get_axis_number(axis)
          return np.meshgrid(*self.bin_centers(), indexing='ij')[axis]
docs: Return ndarray of same shape as histogram containing bin center value along axis at each point