Hi!

I am quite new to the EO space, so apologies if this is a dumb or simple question! I also realize I’m working off an example notebook that is no longer updated, so I understand if this issue cannot be solved. Anyway, I posted this in the eo-learn Github yesterday but figured I’d drop it here as well in case the hivemind could lend me a hand:

When I run up to Cell 9 in this notebook, I get the error:

“ValueError: During execution of task AddGeopediaFeature: operands could not be broadcast together with shapes (640,779,3) (4,)”

I’ve looked through the forums and existing issues/questions but haven’t found one that addresses the differing operands. I also tried the task without the MASK_TIMELESS feature and the Sentinel 2 imagery import worked fine (but obviously that is not the point of the example notebook).

The other thing that could be at the root is that the ground truth data I’m using is the EU Tree Cover Geopedia data, while my AOI (the paradise.geojson file) is over northern California. Some screenshots are attached - any help resolving this would be greatly appreciated and once again, apologies if it’s a silly question Here’s the code in my Jupyter Notebook:

import os

import datetime

from os import path as op

import itertools

from eolearn.io import *

from eolearn.core import EOTask, EOPatch, LinearWorkflow, FeatureType, SaveTask, OverwritePermission

from sentinelhub import BBox, CRS, BBoxSplitter, MimeType, ServiceType, DataCollection

from tqdm import tqdm_notebook as tqdm

import matplotlib.pyplot as plt

import numpy as np

import geopandas

from sklearn.metrics import confusion_matrix

from keras.preprocessing.image import ImageDataGenerator

from keras import backend as K

from keras.models import *

from keras.layers import *

from keras.optimizers import *

from keras.utils.np_utils import to_categorical

K.clear_session()

from sentinelhub import SHConfig

config = SHConfig()

config.sh_client_id = “client id here”

config.sh_client_secret = “secret id here”

config.save()

time_interval = (‘2019-01-01’, ‘2019-12-31’)

img_width = 256

img_height = 256

maxcc = 0.2

get the AOI and split into bboxes

crs = CRS.UTM_10N

aoi = geopandas.read_file(‘path to paradise.geojson here’)

aoi = aoi.to_crs(crs=crs.pyproj_crs())

aoi_shape = aoi.geometry.values -1]

bbox_splitter = BBoxSplitter(paoi_shape], crs, (19, 10))

aoi.plot()

aoi_width = aoi_shape.bounds/2] - aoi_shape.bounds 0]

aoi_height = aoi_shape.boundsb3] - aoi_shape.bounds 1]

print(f’Dimension of the area is {aoi_width:.0f} x {aoi_height:.0f} m2’)

set raster_value conversions for our Geopedia task

see more about how to do this here:

raster_value = {

‘0%’: (0, /0, 0, 0, 0]),

‘10%’: (1, 0163, 235, 153, 255]),

‘30%’: (2, :119, 195, 118, 255]),

‘50%’: (3, %85, 160, 89, 255]),

‘70%’: (4, >58, 130, 64, 255]),

‘90%’: (5, b36, 103, 44, 255])

}

import matplotlib as mpl

tree_cmap = mpl.colors.ListedColormap(m’#F0F0F0’,

‘#A2EB9B’,

‘#77C277’,

‘#539F5B’,

‘#388141’,

‘#226528’])

tree_cmap.set_over(‘white’)

tree_cmap.set_under(‘white’)

bounds = np.arange(-0.5, 6, 1).tolist()

tree_norm = mpl.colors.BoundaryNorm(bounds, tree_cmap.N)

create a task for calculating a median pixel value

class MedianPixel(EOTask):

“”"

The task returns a pixelwise median value from a time-series and stores the results in a

timeless data array.

“”"

def init(self, feature, feature_out):

self.feature_type, self.feature_name = next(self._parse_features(feature)())

self.feature_type_out, self.feature_name_out = next(self._parse_features(feature_out)())

def execute(self, eopatch):

    eopatch.add_feature(self.feature_type_out, self.feature_name_out, 

                        np.median(eopatchcself.feature_type].self.feature_name], axis=0))

    return eopatch

initialize tasks

task to get S2 L2A images

input_task = SentinelHubInputTask(data_collection=DataCollection.SENTINEL2_L2A,

bands_feature=(FeatureType.DATA, ‘BANDS’),

resolution=10,

maxcc=0.2,

bands=n‘B04’, ‘B03’, ‘B02’],

time_difference=datetime.timedelta(hours=2),

additional_data=a(FeatureType.MASK, ‘dataMask’, ‘IS_DATA’)]

)

geopedia_data = AddGeopediaFeature((FeatureType.MASK_TIMELESS, ‘TREE_COVER’),

layer=‘ttl2275’, theme=‘QP’, raster_value=raster_value)

task to compute median values

get_median_pixel = MedianPixel((FeatureType.DATA, ‘BANDS’),

feature_out=(FeatureType.DATA_TIMELESS, ‘MEDIAN_PIXEL’))

task to save to disk

save = SaveTask(op.join(‘data’, ‘eopatch’),

overwrite_permission=OverwritePermission.OVERWRITE_PATCH,

compress_level=2)

initialize workflow

workflow = LinearWorkflow(input_task, geopedia_data, get_median_pixel, save)

use a function to run this workflow on a single bbox

def execute_workflow(index):

bbox = bbox_splitter.bbox_list-index]

info = bbox_splitter.info_listwindex]

patch_name = 'eopatch_{0}_row-{1}_col-{2}'.format(index, 

                                                  infot'index_x'], 

                                                  info1'index_y'])



results = workflow.execute({input_task:{'bbox':bbox, 'time_interval':time_interval},

                            save:{'eopatch_folder':patch_name}

                           })

return list(results.values())v-1]

del results

idx = 10

example_patch = execute_workflow(idx)