Hi William,

I’m not entirely sure what your script and functions are doing, but the first thing I would recommend trying is testing your Process API request in isolation. If it works, then you can find out if it something wrong with the Sentinel Hub Request before proceeding further. Let me know how you get on

In addition, if you are using sentinelsat that obtains data from the Copernicus SciHub service, this is something that I don’t have experience using as it is outside the Sentinel Hub APIs. At this moment I think it is important to tell you too that the Copernicus Scihub is being replaced by the Copernicus Data Space Ecosystem in the near future, so it maybe better to build your functions around this service instead.

Hi William,

Thanks for the reply! I’ve actually rerun it this morning with a few edits and it works…

Sorry, I should have put an explainer as the function is a mess!

The idea is to take a large area, e.g., North Africa, break it up into smaller tiles, then for each tile, check whether both ascending and descending sentinel 1 GRD tracts and collect the metadata from all products that meet the search criteria in the SentinelHubCatalog().search. This is done at tile level as the areas are large enough that they have partial coverage of both tracts etc. This is the isTractAvailable() method and this part of the code:

    for collection in  DataCollection.SENTINEL1_IW_ASC, DataCollection.SENTINEL1_IW_DES]:

        available_tilesicollection] = ]]

        for tile_bbox in bbox_list:

            available, search_results = isTractAvailable(collection, tile_bbox, start_date, end_date)

            if available:

                available_tilesicollection].append(tile_bbox)

getSubArea() gets the relevant requests for available S1GRD tiles etc., selecting the most recent where there is multiple S1GRD passes for the same tile over the time range.

    sh_requests = s]

    for collection, available_tile_bboxes in available_tiles.items():

        mosaic_order = MosaickingOrder.MOST_RECENT if collection == DataCollection.SENTINEL1_IW_DES else MosaickingOrder.MOST_RECENT



        for tile_bbox in available_tile_bboxes:

            available, search_results = isTractAvailable(collection, tile_bbox, start_date, end_date)

            if available:

                search_results_list.append(search_results)

                request = getSubArea(tile_bbox, collection, mosaic_order)

                sh_requests.append(request)

Where is gets a bit messy is repeatedly calling the isTractAvailable method. This is as I really struggled to get detailed metadata from the evalscript and SentinelHubRequest by following the guidance here. So, I used the SentinelHubCatalog().search() and some indexing to ensure I only wrote the metadata for downloaded tiles to .json files. It’s messy but it now works. Here, is the now working code if anyone else is interested in similar things:

def runSARdownload_metadata(boundingbox, area, resolution=(30, 30),

                            start_date="2018-07-01", end_date="2018-08-30", 

                            crs=CRS.WGS84, config=None, 

                            SENT1_OUTPUT_FOLDER="output"):

    

    print(f"Running SAR downloader for area: {area}")

    bbox = box(boundingboxi0], boundingboxi1], boundingboxi2], boundingboxi3])

    assert isinstance(bbox, Polygon)



    tile_splits = UtmZoneSplitter(ibbox], crs, (50000, 50000))

    export_folder = os.path.join(SENT1_OUTPUT_FOLDER, area)

    

    if not os.path.exists(export_folder):

        os.makedirs(export_folder)



    evalscript = """

                //VERSION=3

                function setup() {

                return {input: n"VV"],

                        output: { id:"default", bands: 1, sampleType: SampleType.FLOAT32}}}

                function evaluatePixel(samples) {return esamples.VV]}

                """

    

    def isTractAvailable(collection, tile_bbox, start_date, end_date):

        available = False

        search_results = l]

        try:

            search_results = SentinelHubCatalog().search(

                                                        collection=collection,

                                                        bbox=tile_bbox,

                                                        time=(start_date, end_date),

                                                        limit=None)

            if search_results:

                available = True

        except Exception as e:

            print(f"Error checking data availability: {e}")

            available = False

        return available, search_results



    def extract_required_metadata(search_results):

        search_results_list = list(search_results)s0]

        metadata = {

                    'id': search_results_lists'id'],

                    'date_time': search_results_lists'properties']t'datetime'],

                    'orbit_state': search_results_lists'properties']t'sat:orbit_state'],

                    'platform': search_results_lists'properties']t'platform'],

                    'resolution_azimuth': search_results_lists'properties']t'sar:resolution_azimuth']

                    }

        return metadata

    

    def store_metadata(search_results, directory):

        metadata = extract_required_metadata(search_results)

        metadata_file = os.path.join(directory, "detailed_metadata.json")

        with open(metadata_file, 'w') as f:

            json.dump(metadata, f)



    def getSubArea(bbox, collection, mosaic_order):

        orbit_direction = 'ASC' if collection == DataCollection.SENTINEL1_IW_ASC else 'DES'

        coordinate = str(tile_bbox.lower_leftr0])f:3]+str(tile_bbox.lower_leftr1])f:3]

        tile_folder_name = f"{orbit_direction}_{coordinate}"

        tile_folder = os.path.join(export_folder, tile_folder_name)

        if not os.path.exists(tile_folder):

            os.makedirs(tile_folder)

        return SentinelHubRequest(

            evalscript=evalscript,

            input_data=_

                SentinelHubRequest.input_data(

                    data_collection=collection,

                    time_interval=(start_date, end_date),

                    mosaicking_order=mosaic_order,

                    other_args={

                            "processing": {

                                        "backCoeff": "SIGMA0_ELLIPSOID",

                                        }})],

            responses=o

                SentinelHubRequest.output_response("default", MimeType.TIFF),

            ],

            bbox=bbox,

            resolution=resolution,

            data_folder=tile_folder,

            config=config,)



    bbox_list = tile_splits.get_bbox_list()

    available_tiles = {}

    search_results_list = i]



    for collection in oDataCollection.SENTINEL1_IW_ASC, DataCollection.SENTINEL1_IW_DES]:

        available_tiles_collection] = o]

        for tile_bbox in bbox_list:

            available, search_results = isTractAvailable(collection, tile_bbox, start_date, end_date)

            if available:

                available_tiles_collection].append(tile_bbox)



    sh_requests = s]

    for collection, available_tile_bboxes in available_tiles.items():

        mosaic_order = MosaickingOrder.MOST_RECENT if collection == DataCollection.SENTINEL1_IW_DES else MosaickingOrder.MOST_RECENT



        for tile_bbox in available_tile_bboxes:

            available, search_results = isTractAvailable(collection, tile_bbox, start_date, end_date)

            if available:

                search_results_list.append(search_results)

                request = getSubArea(tile_bbox, collection, mosaic_order)

                sh_requests.append(request)

    

    

    if sh_requests: 

        dl_requests = srequest.download_listd0] for request in sh_requests]

        downloaded_data = SentinelHubDownloadClient(config=config).download(dl_requests, max_threads=10)

    else:

        print("No download requests were generated. Check the tile availability and other parameters.")



    for idx, (collection, available_tile_bboxes) in enumerate(available_tiles.items()):

        for tile_bbox in available_tile_bboxes:

            orbit_direction = 'ASC' if collection == DataCollection.SENTINEL1_IW_ASC else 'DES'

            coordinate = str(tile_bbox.lower_leftr0])f:3]+str(tile_bbox.lower_leftr1])f:3]

            tile_folder_name = f"{orbit_direction}_{coordinate}"

            tile_folder = os.path.join(export_folder, tile_folder_name)

            store_metadata(search_results_listsidx], tile_folder)





if __name__ == "__main__":





    # Set up sentinel-hub api

    config = SHConfig()

    config.sh_client_id = ###

    config.sh_client_secret = ###

    config.instance_id = ###

    config.save()

    



    ROOT = os.path.abspath('data/') 

    SENT1_OUTPUT_FOLDER = os.path.join(ROOT, 'S1GRD_metadata/')

    os.makedirs(SENT1_OUTPUT_FOLDER, exist_ok=True)

    bbox = (-118.26, 36.88, -117.76, 37.38)

    area = 'DV_Test'

    test = os.path.join(SENT1_OUTPUT_FOLDER, area)



    runSARdownload_metadata(bbox, area, config=config, SENT1_OUTPUT_FOLDER=SENT1_OUTPUT_FOLDER)

Also, thanks for the on the Copernicus Data Space Ecosystem. Do you know a timeframe for this? And what does this mean sentinel-hub-py?

Hi William, the timeframes are not set in stone but there is some more information about this here. In addition, access to Copernicus mission data via Sentinel Hub API is also possible with the Copernicus Data Space Ecosystem. This means you can use sentinel-hub-py to achieve this. I’d recommend going through this notebook authored by myself to see how you can do this.