How do I select 10*10 pixels for a particular resolution centered on a certain coordinate?

[lizhuoq]

Take the example of downloading Sentinel 1 Radiometrically Terrain Corrected (RTC) data

[ghidalgo3]

@lizhuoq starting with the example notebook, you can do something like this:

import requests

point = [-112.4627, 36.3516]
point_feature = {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "coordinates": point,
        "type": "Point"
      }
    }
# this box surrounds the previous point,
# you will probably want to construct this box from the point coordinates
box = {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "coordinates": [
          [
            [
              -112.50775947420826,
              36.38271075702113
            ],
            [
              -112.50775947420826,
              36.3267430650056
            ],
            [
              -112.42430457762461,
              36.3267430650056
            ],
            [
              -112.42430457762461,
              36.38271075702113
            ],
            [
              -112.50775947420826,
              36.38271075702113
            ]
          ]
        ],
        "type": "Polygon"
      }
    }
search = catalog.search(
    collections=["sentinel-1-rtc"], datetime="2025-02-01/2025-02-19", max_items=5, intersects=point_feature
)

item = list(search.items())[0]
# 10, 10 is really small, I made it 100, 100 here to make it easier to see
width, height, format = 100, 100, "png"
response = requests.post(
    f"https://planetarycomputer.microsoft.com/api/data/v1/item/feature/{width}x{height}.{format}",
    params={
        "collection": item.collection_id,
        "item": item.id,
        "asset_as_band": True,
        "assets": ["vv", "vh"],
        "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,.80000", "0,1.0", "0,1.0"]
    },
    json=box
)
Image(response.content)

That will produce a 100x100 image:

You may find that other Planetary Computer Data APIs are more suitable for your use case. Good luck!

[lizhuoq]

@ghidalgo3 Thx