docs: def apply_colormap(image, colormap, contig=True):
      """Return palette-colored image.
  
      The image values are used to index the colormap on axis 1. The returned
      image is of shape image.shape+colormap.shape[0] and dtype colormap.dtype.
  
      Parameters
      ----------
      image : numpy.ndarray
          Indexes into the colormap.
      colormap : numpy.ndarray
          RGB lookup table aka palette of shape (3, 2**bits_per_sample).
      contig : bool
          If True, return a contiguous array.
  
      Examples
      --------
      >>> image = numpy.arange(256, dtype='uint8')
      >>> colormap = numpy.vstack([image, image, image]).astype('uint16') * 256
      >>> apply_colormap(image, colormap)[-1]
      array([65280, 65280, 65280], dtype=uint16)
  
      """
      image = numpy.take(colormap, image, axis=1)
      image = numpy.rollaxis(image, 0, image.ndim)
      if contig:
          image = numpy.ascontiguousarray(image)
      return image