docs: def execute(self, eopatch):
          """ Execute function which adds new vector layer to the EOPatch
  
          :param eopatch: input EOPatch
          :type eopatch: EOPatch
          :return: New EOPatch with added vector layer
          :rtype: EOPatch
          """
          bbox_map = self._get_submap(eopatch)
          height, width = self._get_shape(eopatch)
          data_transform = rasterio.transform.from_bounds(*eopatch.bbox, width=width, height=height)
  
          if self.feature_name in eopatch[self.feature_type]:
              raster = eopatch[self.feature_type][self.feature_name].squeeze()
          else:
              raster = np.ones((height, width), dtype=self.raster_dtype) * self.no_data_value
  
          if not bbox_map.empty:
              rasterio.features.rasterize([(bbox_map.cascaded_union.buffer(0), self.raster_value)], out=raster,
                                          transform=data_transform, dtype=self.raster_dtype)
  
          eopatch[self.feature_type][self.feature_name] = raster[..., np.newaxis]
  
          return eopatch