Hey, sorry for a bit late answer.

The problem is that using m.add_wms_layer() directly, the additional params are not used by the underlying ipyleaflet’s WMSLayer class.

The solution is to create a custom class which extends ipyleaflet’s WMSLayer class with additional params needed for Sentinel Hub WMS request.

Then we can create a layer from that class and add it to the map.

I followed an example in ipyleaflet’s documentation for WMS layer and added time and geometry. I also added maxcc (max cloud coverage) parameter, so it overrides the value that is set in the dashboard for the layer in case you would want to do that.

There is no need to set the bounding box manually because it gets set automatically and changes when you interact with the map (zoom / move).

There is no need to extend the default ipyleaflet’s WMSLayer with the crs parameter because it already accepts it - it just needs to be of correct type. (I followed this clue in the ipyleaflet’s documentation)

The code

import leafmap.leafmap as leafmap

import localtileserver

from ipyleaflet import WMSLayer, projections

from traitlets import Unicode



url = "https://services.sentinel-hub.com/ogc/wms/<instance>"

p = "POLYGON ((38.193 12.236, 38.193 12.249, 38.180 12.249, 38.180 12.236, 38.193 12.236))"



# custom class which extends ipyleaflet's WMSLayer class with

# additional params needed for Sentinel Hub WMS request

class SHWMSLayer(WMSLayer):

    time = Unicode('').tag(sync=True, o=True)

    geometry = Unicode('').tag(sync=True, o=True)

    maxcc=Unicode('').tag(sync=True, o=True)



sh_wms = SHWMSLayer(

    url=url,

    layers="1_TRUE_COLOR",

    format='image/png',

    transparent=True,

    geometry=p,

    time="2022-05-15/2022-05-20",

    crs=projections.EPSG4326,

    maxcc="100"

)



m = leafmap.Map(zoom=16, center=612.2425, 38.1865])

m.add_basemap("SATELLITE")

m.add_layer(sh_wms)

m

Hope this helps. If you have any additional questions, don’t be afraid to ask.

Cheers

Thanks very much, that works perfectly!

In case of interest for others, I have deployed the solution here:

GitHub - agroimpacts/labelreview: Notebook for reviewing labelling work

Notebook for reviewing labelling work. Contribute to agroimpacts/labelreview development by creating an account on GitHub.

agroimpacts/labelreview/blob/main/src/labelreview.py

import pandas as pd

import geopandas as gpd

import leafmap.leafmap as leafmap

from ipyleaflet import WMSLayer, projections

from traitlets import Unicode

import localtileserver

import psycopg

import yaml

import os

from pathlib import Path

import pandas.io.sql as sqlio

import sqlalchemy as sa

from shapely.geometry import Polygon, box



class SHWMSLayer(WMSLayer):

    """Custom class to enable subsetting of images read from SentinelHub in 

    leafmap/ipyleaflet, provided by Ziga Cernigoj of SentinelHub

    https://forum.sentinel-hub.com/t/\

        select-area-of-byoc-image-using-bbox-geometry-in-iypleaflet/8211/2

This topic was automatically closed 60 days after the last reply. New replies are no longer allowed.