Develop

How can I delete the black part of the image so I can get just the NDVI in my AOI?

Hi, I’m trying to use the python module to acquire the S-2 NDMI for a specific area, but when I save them to my local drive they have some sort of glitch, like it was a result of merging two images. I checked the EO Browser and they seem fine. It only happens using the Python module. Image is attached.

When looking at this link - https://docs.sentinel-hub.com/api/latest/data/sentinel-1-grd/examples/#s1grd-orthorectified-decibel-gamma0-vh-between--20-db-and-0-db-png I can see examples of process API for Sentinel-1. When I look at this example: S1GRD orthorectified linear gamma0 VV, ascending orbit direction, GeoTIFF in EPSG:32648 (UTM zone 48N) I can see for example the property: “backCoeff”: “GAMMA0_ELLIPSOID”. Is there any documentation where I can see which other options I can use for the property “backCoeff”? or any other property for that matter (e.g., orbitDirection

Hi team. I was trying to download Sentinel-2 L2A files. I used this code. request = AwsTileRequest( tile=tile_name, time=date, aws_index=aws_idx, data_folder=f"./{year}", data_collection=DataCollection.SENTINEL2_L2A, metafiles=["tileInfo", "qi/MSK_CLOUDS_B00"] ) This code gave me image with cloud gml file before 2022. However, it gave me “SHRuntimeWarning: File in location s3://sentinel-s2-l2a/tiles/16/S/CG/2022/7/31/0/qi/MSK_CLOUDS_B00.gml is missing” error messages when I tried to download L2A image in 2022. How can this problme be solved? Thanks a lot.

Hi everyone,I have been trying to find out if it is possible to access Sentinel-3 Bottom of Atmosphere Reflectances using the Sentinel Hub services but haven’t been successful in finding either a positive or negative answer so far. Does anyone know if this is possible?I could also download L1 radiances, convert to reflectances and then apply an atmospheric correction with e.g. the Rayleigh correction within the ESA SNAP software. It looks like I could access Sentinel-3 OLCI L1B data but apparently already processed and not in the native netCDF file format (in .SEN3 folders) which I would need in order to use the data with the SNAP software.Thank you a lot in advance for your help and advice!best

I have been working on obtaining cloudless mosaic images but have been struggling with strong orbit-wise artifacts in some areas (see my post). As an alternative to seamless mosaics I am interested in obtaining orbit-wise cloudless images. The cloud-filling results in my tests so far look really good but the stripes introduced by the mosaicking process are going to complicate downstream applications of the imagery. Is there a way to specify a request to only pull images from a particular orbit rather than sourcing images from multiple orbits? For example, in the image below I would only want to pull raw images from the center orbit (relative orbit 54) and avoid pulling images from relative orbits 97 and 11.

Hello. We are trying to calculate the Processing unit per hectare. Here is the example which we use for our calculations. Please, can you give me some feedback if we calculate it correctly for NDVI?Your example on website:Output size: 0,01 (20x20 px)Input bands : 2/3Output format: 1Number of data sample : 1Orthorecification : 1PU = 0,01 x 2/3 x 1 x 1 x 1 = 0,0067 PUThe example on your website is designed for 20x20 px = 4h a = that is why Output size / multiplication factor : 0,01.Also you define that the minimum value of multiplication factor is 0,01 which corresponds to an are of 0,25km2 = 25 ha.Does it mean that the processing unit for 25ha in same example is the same ? 0,0067 PU ?If I am correct, then processing unit for 250 ha would be : 0,0067 x 10 = 0,067 PU. Is that correct ?Thank you for answering

I’m wondering if you could point me in the right direction with re to the best way to download S1 and S2 data - I have a list of S1 and S2 tiles I want to download and I would like to pass on that list to query what’s in the archive and move to my S3 bucket. I’ve seen the jupyter notebook on batch processing heregithub.com sentinel-hub/sentinelhub-py/blob/master/examples/batch_processing.ipynb { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Sentinel Hub Batch Processing\n", "\n", "A tutorial about [Large area utilities](./large_area_utilities.ipynb) shows how to split a lar

Hi everyone, I have I problem. I was downloading data using dhusget. Now it starts to be the past. So I want to download native .nc file of Sentinel5P L2 products. Lets say CO. I have accounts credentials with client id and secret key and still don’t know how to do this. I have a following code wich actually find the products but I don’t know that to do next. here is the code #!/usr/bin/env python3 import requests url = "https://creodias.sentinel-hub.com/api/v1/catalog/1.0.0/search" headers = { "Content-Type": "application/json", "Authorization": "Bearer <mytoken>" } data = { "collections": [ "sentinel-5p-l2" ], "datetime": "2020-04-01T00:00:00Z/2020-04-01T23:59:59Z", "bbox": [ -9.037709, 38.822591, 35.272283, 67.609221 ], "limit": 10, "filter": { "op": "and", "args": [ { "op": "=", "args": [ { "property": "s5p:type" }, "C

Hello, I am trying to compute the processing units that will be spent based on user request while using process and statistical APIs as per the processing units definition page (https://docs.sentinel-hub.com/api/latest/api/overview/processing-unit/). I wish to calculate the processing units that will be spent on user request. In this regard, I wish to compute the size of the requested image based on the bounding box input by the user. The request builder option on sentinel hub shows the output size of the image in width and height. I wish to compute the same using the user defined bbox in my code for a fixed resolution of the image, so that I can compute the processing units that will be spent on the basis of the size of the image. Could anybody please guide me in this regard to compute the size of the image (no. of pixels in width and height) using the bounding box and thereby in computing the processing units. I have attached screenshots of the output th

Dear members in the forum, In the description of the ARI1 (Anthocyanin Reflectance Index), it is written that its values ranges from 0 to 0.2 (see attached screenshot). But when I used its given formula [(1/R550)-(1/R700)], I get values from -3 to 17. I checked images from different dates but the values are not in the range from 0 to 0.2. What could be possible cause of wrong range, does it require any conversions?

Been downloading Sentinel 2 imagery on the Copernicus Open Access Hub for a while now. However since about 3 days ago, every request for retrieval after adding my imagery to the cart says “status failed”. Any idea why and what I can do about it?

Hi, I’m using WCS request in Python to query PlanetScope data that I recently ordered. I can see the data in EO Browser, however, when I run a request in Python, I receive zero arrays as a response. I can, however, iterate via timestamps, so I think the query goes thorough (so the polygon intersect with purchased imagery), however, it won’t actually return correct values. As an input, I feed correct instance and collection IDs, as well as am referring to a pre-defined layer (True Color). Any idea of what might be the issue here? Thanks! Edit: I’ve just ran a similar request for SPOT data where I only substituted instance and collection IDs, and it went through without issues, so I would assume there is some sort of issue with PlanetScope? Edit2: Below is an example of WCS request I’m using: config_ps = SHConfig() INSTANCE_ID_ps = ‘d46e3de6-2ee5-masked’ config_ps.instance_id = INSTANCE_ID_ps collection_id_ps= ‘2

Hi Team! I kindly need your guidance on something that regards Sentinel Hub. I’m totally new to this and I wanted to develop a small web GIS platform that will be later used for crop monitoring. Right now, I just needed to develop a portion of it that can allow me to either import shapefile or json and also be able to draw a polygon in my area of interest. I also wanted to be able to extract the NDVI of the area of interest and at least be able to plot the NDVI values say over time. I will appreciate any assistance. I would want to know what tools I can use to develop the said portion Thanks!

What is the most efficient method to create a WMS service showing the change between two dates? I have successfully set up a NDVI WMS service for a single date, and understand the “evaluate pixel” evaluates the pixel for the date. But is there a way to evaluate the pixel for two dates to get the change of the NDVI, and serve that as the layer in the WMS? If this is not possible, then what would be the best approach to analyze the change between two dates? My thought would be getting the array from both dates and doing the analysis once data is received from the service, but before I do so I thought I’d pose the question first. Due to some inexperience in Sentinel Hub Python API usage, I thought it may be best to ask on this forum 😃 Thank you in advance.

Hello everyone, I have some troubles with the QGIS plug-in. I need to download maps from different ranges with different combinations of variables (2-3-4, 3-4-8, 5-6-7, 8a-11-12) to study a forest during growing season. Unfortunately, it ain’t work : the app says “The bounding box area is too large! Please zoom in.” That doesn’t really makes sense since I select a very small area compared to some I already selected before. I use a shp with geolocated points as reference. I also tried with the OpenStreetMap background map and zooming, which gave no different result. I even tried to use a custom bounding box but that way either it ain’t work. My subscription is active, and I also tried to download a simple NDVI card to check : it doesn’t work either. If anyone has a solution I will be grateful, I’m totally stucked for now ! I apologize if it’s a very easy to deal with mistake, I’m just beginning with sci-hub. Best regards !

Is it possible to download cloud probabilities/cloud masks with the sentinelhub.py API ?

I’m trying to use sentinelhub package on jupyterlab notebook and for some reason I don’t manage to be able to acess images. When I use my credentials in the configuration at the beginning it works and pronts the configuration file location, but when I try to acess images as descrribed here I get error message: DownloadFailedException: Failed to download from: long link of sentinel-hub with HTTPError: 404 Client Error: Not Found for urlanother long link of sentinelhub Server response: “Invalid UUID string: my credentials” This is how I have tried to acess the images: !sentinelhub.config --instance_id MySecretCredentials !sentinelhub.config --show #that works, print the configurations from sentinelhub import WebFeatureService,

while I request https://services.sentinel-hub.com/api/v1/statistics,it responses {“error”:{“status”:403,“reason”:“Forbidden”,“message”:“Invalid or expired account.”,“code”:“COMMON_INSUFFICIENT_PERMISSIONS”}},why?

I would like to get a Sentinel-1 GRD IW image with a resolution (pixel spacing) of 10 m. I followed the documentation here https://docs.sentinel-hub.com/api/latest/data/sentinel-1-grd/#resolution-pixel-spacing and added the following arguments: "resolution":"HIGH", "acquisitionMode":"IW", "polarization":"DV" }, "processing": { "orthorectify": "true", "backCoeff": "SIGMA0_ELLIPSOID", "demInstance":"COPERNICUS_30" But the image that I get have weird-looking pixels:

I am using the subscription API for subscribing AOI to Planet but getting BAD Request in response Request: POST https://services.sentinel-hub.com/api/v1/dataimport/subscriptions { "name": "Headland-HL1", "collectionId": "4f7d7518-9333-MASKED", "provider": "PLANET", "input": { "bounds": { "geometry": { "type": "Polygon", "coordinates": [ [ [ 118.5200492301355, -20.35358590810601, 0.0 ], [ 118.525520243307, -20.35208216090831, 0.0 ], [ 118.5315806188978, -20.35010044578434, 0.0 ], [ 118.5367336293801, -20.34862457841812, 0.0 ],

Hello, I wanted to use the B10 band from Sentinel-2 L2A AWS datastore, but I cannot find an image file for this band. Base on this article, the B10 band should be in 60m resolution, but there is no such image when I list the images from AWS. $ aws s3 ls --request-payer requester s3://sentinel-s2-l2a/tiles/22/K/GV/2022/2/28/0/R60m/ 2022-02-28 21:30:02 580081 AOT.jp2 2022-02-28 21:30:02 3749625 B01.jp2 2022-02-28 21:30:02 3736007 B02.jp2 2022-02-28 21:30:02 3752443 B03.jp2 2022-02-28 21:30:02 3773814 B04.jp2 2022-02-28 21:30:02 3753745 B05.jp2 2022-02-28 21:30:02 3742935 B06.jp2 2022-02-28 21:30:02 3753147 B07.jp2 2022-02-28 21:30:02 4398952 B08.jp2 2022-02-28 21:30:02 3763914 B09.jp2 2022-02-28 21:30:02 3739951 B11.jp2 2022-02-28 21:30:02 3752529 B12.jp2 2022-02-28 21:30:0

Hi there, I know that the major difference between L2A and L1C data is that the former is atmoshperically corrected. When I howeverd download L1C data from Sentinelhub using a WMS request I can also get L1C data atmospherically corrected. Simply by using custom_url_params={CustomUrlParam.ATMFILTER: ‘ATMCOR’}. So what is then is the difference of requesting L1C data with CustomUrlParam.ATMFILTER: ‘ATMCOR’ and requesting L2A data?

I have list of polygons that are shapely type. These polygons are different boudning boxes that I have slected for certain process: I want to have these polygons as bbox type. For that I have tried to use the get_bbox_from_shape function, but that doesn’t work as I understand it needs to have sentinelhub geometry: def get_bbox_from_shape(shape,resolution): minx, miny, maxx, maxy = shape.geometry.total_bounds bbox_coords_wgs84=[minx, miny, maxx, maxy] bbox = BBox(bbox=bbox_coords_wgs84, crs=CRS.WGS84) bbox_size = bbox_to_dimensions(bbox, resolution=resolution) return bbox_size,bbox,bbox_coords_wgs84 get_bbox_from_shape(bbox_list[0],10) AttributeError Traceback (most recent call last) in ----> 1 get_bbox_

Hi, I’m trying to use SentinelHub in QGIS and it is working both with WMS and WCS. However, I do not know how to get QGIS to manipulate the TIME parameter of the web-service requests, such that I can control the date of which images I will show in QGIS. Do anyone know how to handle this in QGIS?