code: def blend(self, blend_function=stack):
          
          new_scn = Scene()
          common_datasets = self.shared_dataset_ids
          for ds_id in common_datasets:
              datasets = [scn[ds_id] for scn in self.scenes if ds_id in scn]
              new_scn[ds_id] = blend_function(datasets)
  
          return new_scn
docs: Blend the datasets into one scene.
  
          .. note::
  
              Blending is not currently optimized for generator-based
              MultiScene.