Skip to main content

Request image with multipolygons

  • 26 April 2024
  • 5 replies
  • 4 views

Hello everyone, I hope you’re all doing well.


I’m currently working on sending a request that involves multiple polygons. However, I’m encountering an issue where there’s a red area appearing between the polygons, when visualizing the NDVI index


Here’s my script :


//VERSION=3

function setup() {
return {
input: "B04", "B08", "SCL", "dataMask"],
output:
{ id: "default", bands: 4 },
{ id: "index", bands: 1, sampleType: "FLOAT32" },
{ id: "eobrowserStats", bands: 2, sampleType: "FLOAT32" },
{ id: "dataMask", bands: 1 },
],
};
}

function evaluatePixel(samples) {
let val = index(samples.B08, samples.B04);
let imgVals = null;

const indexVal = samples.dataMask === 1 ? val : NaN;

if (val < -0.5) imgVals = i0.05, 0.05, 0.05, samples.dataMask];
else if (val < -0.2) imgVals = i0.75, 0.75, 0.75, samples.dataMask];
else if (val < -0.1) imgVals = i0.86, 0.86, 0.86, samples.dataMask];
else if (val < 0) imgVals = i0.92, 0.92, 0.92, samples.dataMask];
else if (val < 0.025) imgVals = i1, 0.98, 0.8, samples.dataMask];
else if (val < 0.05) imgVals = i0.93, 0.91, 0.71, samples.dataMask];
else if (val < 0.075) imgVals = i0.87, 0.85, 0.61, samples.dataMask];
else if (val < 0.1) imgVals = i0.8, 0.78, 0.51, samples.dataMask];
else if (val < 0.125) imgVals = i0.74, 0.72, 0.42, samples.dataMask];
else if (val < 0.15) imgVals = i0.69, 0.76, 0.38, samples.dataMask];
else if (val < 0.175) imgVals = i0.64, 0.8, 0.35, samples.dataMask];
else if (val < 0.2) imgVals = i0.57, 0.75, 0.32, samples.dataMask];
else if (val < 0.25) imgVals = i0.5, 0.7, 0.28, samples.dataMask];
else if (val < 0.3) imgVals = i0.44, 0.64, 0.25, samples.dataMask];
else if (val < 0.35) imgVals = i0.38, 0.59, 0.21, samples.dataMask];
else if (val < 0.4) imgVals = i0.31, 0.54, 0.18, samples.dataMask];
else if (val < 0.45) imgVals = i0.25, 0.49, 0.14, samples.dataMask];
else if (val < 0.5) imgVals = i0.19, 0.43, 0.1972, samples.dataMask];
else if (val < 0.55) imgVals = i0.13, 0.38, 0.07, samples.dataMask];
else if (val < 0.6) imgVals = i0.06, 0.33, 0.04, samples.dataMask];
else imgVals = i0, 0.27, 0, samples.dataMask];

return {
default: imgVals,
index: indexVal],
eobrowserStats: sval, isCloud(samples.SCL) ? 1 : 0],
dataMask: dsamples.dataMask],
};
}

function isCloud(scl) {
if (scl == 3) {
// SC_CLOUD_SHADOW
return false;
} else if (scl == 9) {
// SC_CLOUD_HIGH_PROBA
return true;
} else if (scl == 😎 {
// SC_CLOUD_MEDIUM_PROBA
return true;
} else if (scl == 7) {
// SC_CLOUD_LOW_PROBA
return false;
} else if (scl == 10) {
// SC_THIN_CIRRUS
return true;
} else if (scl == 11) {
// SC_SNOW_ICE
return false;
} else if (scl == 1) {
// SC_SATURATED_DEFECTIVE
return false;
} else if (scl == 2) {
// SC_DARK_FEATURE_SHADOW
return false;
}
return false;
}


Hi,

Could you try to set the format of output response to "image/png" in your request body? The transparent band only works for PNG.

"output": {
"width": 512,
"height": 296.512,
"responses": [
{
"identifier": "default",
"format": {
"type": "image/png"
}
}
]
}

I greatly appreciate it; the script now performs exactly as intended.

 


Hy, I am using similar logic but there is one edge case. If I have 2 Polygons and they are intersecting with each other than the process api through 400 error. “geometry Polygon rings are intersecting”. Any idea how to solve that. FYI I can’t update the data


The only way is to try to solve the geometry issues before passing the data to the Processing API


Hi,

The API only takes valid MultiPolygon and a valid MultiPolygon may not collect any overlapping polygons (defined in shapely).

I would recommend making the overlapping polygons into one polygon in this case.


Reply