code: def cira_stretch(img, **kwargs):
      
      LOG.debug("Applying the cira-stretch")
  
      def func(band_data):
          log_root = np.log10(0.0223)
          denom = (1.0 - log_root) * 0.75
          band_data *= 0.01
          band_data = band_data.clip(np.finfo(float).eps)
          band_data = xu.log10(band_data)
          band_data -= log_root
          band_data /= denom
          return band_data
  
      return apply_enhancement(img.data, func)
docs: Logarithmic stretch adapted to human vision.
  
      Applicable only for visible channels.