Maximum size of Image

Apps / Services Python SDK
1

Good evening,

I want to ask if it is posible to know how can I have a maximum of size (default size is 256 x 256).
I want to know how can I get maximum please.
(size=[512, 512])

Thanks

2

Hi @antonio.santos.rodri ,

the maximum size of the image you can get with Process API is 2500 * 2500 pixel. Note that there’s also a limitation on the size of a pixel, which is 1500 * 1500 meter per pixel.

3

Good Morning,

I have some few question:

  1. I want the max resolution. I tried to calculate but the result of the resolution was very small:
    My question is how can I get max resolution using x and y.
resolucion=5
al=(caja[2]-caja[0])*111.2
an=(caja[3]-caja[1])*111.2
w=int(an/resolucion)
h=int(al/resolucion)
  1. Another question that I have I used index to get my image.
    TCARIOSAVI->IDB - Index: TCARI/OSAVI
    NDII-> IDB - Information for Sensor and Index
    NDRE-> Band of rededge
    MSAVI → (0.5)((2800nm+1)−sqrt(2800nm+1)2−8(800nm−670nm))
    NDWI2 → (857nm−1241nm)/(857nm+1241nm)

My question is that I don’t know equivalent band that is in planet (I use collection of fincas Bustos).

  1. Another question I used script of python to create subscription but I didn’t find Information to create new collection. How can I create new one?
  1. Last question I get the image using script of python the result the image is square. I want just polygon.

#!/usr/bin/env python3

-- coding: utf-8 --

“”"
Created on Wed Aug 10 15:15:02 2022

@author: oumnia
“”"

This is script may only work with sentinelhub.version >= ‘3.4.0’

from sentinelhub import SentinelHubRequest, DataCollection, MimeType, CRS, BBox, SHConfig, Geometry
import os
from PIL import Image
import csv

“”"
parametros poligono tienen qu
“”"

def buscarImg(poligono, fechaIni, fechaFin, nombreCarpeta, tipo):
evalscript = “”"
//VERSION=3
function setup() {
return {
input: [{“bands”: [“Blue”, “Green”, “Red”]}],
output: { bands: 3}
}
}
function evaluatePixel(sample) {
return [2.5 * sample.Red / 10000,
2.5 * sample.Green / 10000,
2.5 * sample.Blue / 10000]
}
“”"
caja=obtenerCaja(poligono)
“”"resolucion=5
al=(caja[2]-caja[0])*111.2
an=(caja[3]-caja[1])*111.2

w=int(an/resolucion)
h=int(al/resolucion)
"""

    
#bbox = BBox(bbox=[-5.7585989, 37.3886755, -5.7565442, 37.3910112], crs=CRS.WGS84)
#geometry = Geometry(geometry={"type":"Polygon","coordinates":[[[-5.7580946,37.3905082],[-5.7579981,37.3903292],[-5.7585989,37.3902355],[-5.7585453,37.3893489],[-5.7583629,37.3886755],[-5.7573811,37.3887778],[-5.7565442,37.3910112],[-5.7580946,37.3905082]]],"type":"Polygon"}, crs=CRS.WGS84)
bbox = BBox(bbox=caja, crs=CRS.WGS84)
geometry = Geometry(geometry={"type":"Polygon","coordinates":poligono,"type":"Polygon"}, crs=CRS.WGS84)

request = SentinelHubRequest(
    data_folder=nombreCarpeta,
    evalscript=evalscript,
    input_data=[
    SentinelHubRequest.input_data(
        data_collection=DataCollection.define_byoc('---'),          
        time_interval=('2021-09-30', '2021-10-10'),          
    ),
],
    responses=[
        SentinelHubRequest.output_response('default', MimeType.TIFF),
    ],
    bbox=bbox,
    geometry=geometry,
    size=[512, 512],
    config=config
)
return request
4

Hi @antonio.santos.rodri ,

  1. The provided script is calculating height and width using a resolution of 5 m. To get the height and width with a resolution of 3 m (highest resolution of PlanetScope), you could use the same script and change the resolution to 3.

  2. The band info of PlanetScope can be found here and here. You could determine the corresponding band by the wavelength. Here’s an example script for NDVI. Note that not all bands required for your interested indices are provided by PlanetScope.

  3. The BYOC collections section in the Bring Your Own COG example demonstrate how to create a new collection using Python scripts.

  4. Your script should return an image similar to Fig 1. The pixels outside the specified polygon should be masked. If you are trying to get a TIFF or a Shapefile having a boundary of your polygon, it is not supported unfortunately.

Screenshot 2022-09-07 at 10.46.39
Screenshot 2022-09-07 at 10.46.39868×900 62 KB

Fig 1

5

Thanks for answer me.

  1. If I don’t create collection they put that is optional so if I don’t create it they could create it automatically?
    1.1. I can have more than one subscription in the same collection?
  1. In my collection of finca Busto I have 4 band but I need 8 how can I update it to 8?
  1. Last question, Getting img TIFF for just polygon it isn’t supported? (jpeg,png it is supported)

Thanks
Regards

6

Hi @antonio.santos.rodri ,

  1. Could you please rephrase the question? A collection will be automatically created if you don’t specify a collection id when ordering the data (see Import data into existing BYOC collection).

    1.1. You can have more than one subscription in the same collection. In general we suggest to have one collection per data type (e.g. PlanetScope, Pleiades, SPOT, WorldView) and reuse it every time you order data of corresponding type.

  2. You order Item Type PSScene4Band so there’s only 4 band in your data. To have 8 band please order Item Type PSScene.

  3. You still get a rectangular image when specifying JPEG and PNG as output format (See the image below), but there’s no data outside the polygon. If this is what you desire in TIFF format, you can simply specify TIFF as output and get a TIFF that has data only inside your polygon.

    31f86d9e-5497-44fc-be16-2592de7ebda3
    31f86d9e-5497-44fc-be16-2592de7ebda3512×528 68.8 KB

7
  1. There is way to change Item Type [PSScene4Band] to Item Type [PSScene] in my collection FincaBustos? (I need 8 band to create some index)
8

Hi @antonio.santos.rodri ,

The easiest way is to reorder the PSScene data. There won’t be extra cost for the same AOI (see here)