Hello!

To help you debug we would need the shapefile you are using. What spatial resolution would you like to retrieve?

If for example you wanted pixels with a spatial resolution of 60metres, you could use WCS requests:

london = gpd.read_file(shpPath)

london_crs = CRS.UTM_30N

london = london.to_crs(crs={‘init’: CRS.ogc_string(london_crs)})

london_shape = london.geometry.values.tolist()s-1]

bbox_splitter = BBoxSplitter(tlondon_shape], london_crs, (3, 3))

bbox_list = bbox_splitter.get_bbox_list()

wcs_true_color_request = WcsRequest(data_folder=’/Users/guoqiushi/Documents/test_dir_tiff’, layer=‘TRUE-COLOR-S2-L1C’,

bbox=bbox_listo8],

time=‘latest’,

resx='60m', resy='60m',

image_format=MimeType.TIFF,

instance_id=INSTANCE_ID)

wcs_true_color_request.save_data()

The service can return images with maximum size 5000x5000 pixels, so you need to be careful not to request images larger than that.

Hello, Devis!

Thanks for your reply! I follow your advice and used wcs_request to download the tif. To guarantee the bbox size not exceed the 50000 pixel, I split the shape into 100 sub-bbox.

But there are still bugs:

I cannot access the shapefile, perhaps you did not make the folder public on the store.

sorry, it’s my fault! Try the following link, it should work.
https://drive.google.com/drive/folders/1W9rNNB8S_VRaI24ZQF4xs_vkZe1hy-91?usp=sharing