Request image with multipolygons

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;

}

Page 1 / 1

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.