Is it possible to use sentinel-1-rtc with the Data API mosaic endpoint?

[scottyhq]

I tried modifying the data api example notebook to use display the following search as a mosaic on a folium map.

https://planetarycomputer.microsoft.com/explore?c=-121.1456%2C47.3259&z=7.77&v=2&d=sentinel-1-rtc&s=false%3A%3A100%3A%3Atrue&ae=0&sr=desc&m=cql%3A257c41fb13818eb6930e12a7bfcd82ed&r=VV%2C+VH+False-color+composite

Rendering single items as tile layers works fine. But I'm unable to get the mosaic to work. I see a lot of 404 and 500 errors in the console but am not noticing anything obviously wrong with the tiles URL:

'https://planetarycomputer.microsoft.com/api/data/v1/mosaic/tiles/5433eb3bbda6fade99420abc2d8f1aec/WebMercatorQuad/{z}/{x}/{y}@1x?collection=sentinel-1-rtc&expression=0.03+++log+%2810e-4+-+log+%280.05+%2F+%280.02+++2+%2A+vv%29%29%29%3B0.05+++exp+%280.25+%2A+%28log+%280.01+++2+%2A+vv%29+++log+%280.02+++5+%2A+vh%29%29%29%3B1+-+log+%280.05+%2F+%280.045+-+0.9+%2A+vv%29%29&rescale=0%2C.8000&rescale=0%2C1.000&rescale=0%2C1.000&asset_as_band=true'

Could this be an issue with the SAS token requirement for RTC compared to the example using sentinel-2 ?

[TomAugspurger]

Looks like a quoting issue. The URLs from the Explorer look like

In [5]: urllib.parse.parse_qs("https://planetarycomputer.microsoft.com/api/data/v1/mosaic/tiles/257c41fb13818eb6930e12a7bfcd82ed/WebMercatorQuad/9/83/179@2x?asset_as_band=true&expression=0.03+%2B+log+%2810e-4+-+log+%280.05+%2F+%280.02+%2B+2+%2A+vv%29%29%29%3B0.05+%2B+exp+%280.25+%2A+%28log+%280.01+%2B+2+%2A+vv%29+%2B+lo
   ...: g+%280.02+%2B+5+%2A+vh%29%29%29%3B1+-+log+%280.05+%2F+%280.045+-+0.9+%2A+vv%29%29&rescale=0%2C.8000&rescale=0%2C1.000&rescale=0%2C1.000&collection=sentinel-1-rtc&format=png")
Out[5]: 
{'https://planetarycomputer.microsoft.com/api/data/v1/mosaic/tiles/257c41fb13818eb6930e12a7bfcd82ed/WebMercatorQuad/9/83/179@2x?asset_as_band': ['true'],
 'expression': ['0.03 + log (10e-4 - log (0.05 / (0.02 + 2 * vv)));0.05 + exp (0.25 * (log (0.01 + 2 * vv) + log (0.02 + 5 * vh)));1 - log (0.05 / (0.045 - 0.9 * vv))'],
 'rescale': ['0,.8000', '0,1.000', '0,1.000'],
 'collection': ['sentinel-1-rtc'],
 'format': ['png']}

While the tilejson generated by the notebook is

In [6]: urllib.parse.parse_qs("https://planetarycomputer.microsoft.com/api/data/v1/mosaic/tiles/ebeaaa9747d8913c11bf077a9e717291/WebMercatorQuad/{z}/{x}/{y}@1x?collection=sentinel-2-l2a&expression=0.03+++log+%2810e-4+-+log+%280.05+%2F+%280.02+++2+%2A+vv%29%29%29%3B0.05+++exp+%280.25+%2A+%28log+%280.01+++2+%2A+vv%29+++lo
   ...: g+%280.02+++5+%2A+vh%29%29%29%3B1+-+log+%280.05+%2F+%280.045+-+0.9+%2A+vv%29%29&rescale=0%2C.8000&rescale=0%2C1.000&rescale=0%2C1.000&asset_as_band=true'")
Out[6]: 
{'https://planetarycomputer.microsoft.com/api/data/v1/mosaic/tiles/ebeaaa9747d8913c11bf077a9e717291/WebMercatorQuad/{z}/{x}/{y}@1x?collection': ['sentinel-2-l2a'],
 'expression': ['0.03   log (10e-4 - log (0.05 / (0.02   2 * vv)));0.05   exp (0.25 * (log (0.01   2 * vv)   log (0.02   5 * vh)));1 - log (0.05 / (0.045 - 0.9 * vv))'],
 'rescale': ['0,.8000', '0,1.000', '0,1.000'],
 'asset_as_band': ["true'"]}

Notice that the expression is missing the + operators. I'm checking on the correct way to quote / escape these.

Edit: This is a bit ugly, but seems to work:

import pystac_client
import requests
import itertools

def key(s):
    return s.split("=")[0]
COLLECTION = "sentinel-1-rtc"
catalog = pystac_client.Client.open(
    "https://planetarycomputer.microsoft.com/api/stac/v1"
)

aoi = {
  "type": "Polygon",
  "coordinates": [
    [
      [-123.60492153752948, 46.04968244060868],
      [-118.68923696942923, 46.04968244060868],
      [-118.68923696942923, 48.57164306792697],
      [-123.60492153752948, 48.57164306792697],
      [-123.60492153752948, 46.04968244060868]
    ]
  ]
}

# Define your temporal range
daterange = {"interval": ["2023-10-01T00:28:21Z", "2023-10-05T23:59:59Z"]}

# Define your search with CQL2 syntax
search = catalog.search(filter_lang="cql2-json", filter={
  "op": "and",
  "args": [
    {"op": "s_intersects", "args": [{"property": "geometry"}, aoi]},
    {"op": "anyinteracts", "args": [{"property": "datetime"}, daterange]},
    {"op": "=", "args": [{"property": "collection"}, COLLECTION]}
  ]
})

r_register = requests.post(
    "https://planetarycomputer.microsoft.com/api/data/v1/mosaic/register",
    json=search.get_parameters(),
)
registered = r_register.json()
registered

tilejson_url = registered["links"][1]["href"]

mosaic_info = requests.get(
    "https://planetarycomputer.microsoft.com/api/data/v1/mosaic/info",
    params=dict(collection=COLLECTION),
).json()
render_config = mosaic_info["renderOptions"][0]
render_config
params = {k: [x.split("=")[1] for x in v] for k, v in itertools.groupby(render_config["options"].split("&"), key=key)}

params["collection"] = COLLECTION

r_tiles = requests.get(
    tilejson_url,
    params=params
)
r_tiles.raise_for_status()
tiles = r_tiles.json()["tiles"][0]
tiles

I don't know why, but Python's urllib.parse.parse_qs is converting the +'s to empty spaces. This probably buggy piece of code seems to do what we need:

def key(s):
    return s.split("=")[0]
params = {k: [x.split("=")[1] for x in v] for k, v in itertools.groupby(render_config["options"].split("&"), key=key)}

[TomAugspurger]

@scottyhq I added a snippet that seems to do what we need. I'll get the notebook updated.

[scottyhq]

excellent, thanks @TomAugspurger !