code: def overlay(top, bottom, maxval=None):
      """Blending two layers.
  
      from: https://docs.gimp.org/en/gimp-concepts-layer-modes.html
      """
      if maxval is None:
          maxval = np.maximum(top.max(), bottom.max())
  
      res = (
          ((2 * top / maxval - 1) * bottom + 2 * top) * bottom / maxval
      )
      return res.clip(min=0)
  
imports: import logging
  import numpy
  from satpy.composites import GenericCompositor
  from satpy.dataset import combine_metadata
label_desc: Use image node_id 1 for calling a global function with example usage: overlay(top, bottom, maxval) and returns: res