0

I need to assess TSM in a river using Sentinel-2. Which is why I need to add the function of TSM as a new band. I used Adding multiple indices as bands to image as a guide but it is not working for my script. It doesn't show any error, but it also doesn't add the band after I map it to the collection.

Below you can find my code:

//Programar SST usando una expresion 
var SST = function (image){
           var addSST =image.expression(
          '-229.34 * ((GREEN/NIR) ** 3 )+ 1001.65 * ((GREEN/NIR) ** 2) - 1422.7 * (GREEN/NIR) +665.7',{
            'GREEN': image.select('B3'),
            'NIR':image.select('B8')
            
          }).rename('SST');
          
        return image.addBands(SST);
  
};

//SE DEBE AÑADIR LA PROPIEDAD SYSTEM:TIME_START
SST = function (image){
        return image.copyProperties(image,['system:time_start']);
  
};


var SST_IC = S2Clipped.map(SST);
print (SST_IC,'Colección Sólidos Suspendidos Totales');

1 Answer 1

1

There are a couple of issues with your script. First, you declare SST, assign it to a function. Then you reassign SST to another function. Only the second function will be used when you're mapping over the function.

Then there are problem with the functions themselves. The first one (in terms of syntax at least) correctly calculates the SST, assigns it to the variable addSST. This is a confusing naming. You typically want to name the function a verb, like addSST, and the variable a noun, like SST. A function does something, a variable is something. The problem is image.addBands(SST), where you're not adding the calculated SST to your image, but the function.

The second function doesn't actually do anything. It takes an image and reassigns the system:time_start property from the same image.

Here's how you could implement this:

var addSST = function(image) {
  var SST = image.expression(
    '-229.34 * ((GREEN/NIR) ** 3 )+ 1001.65 * ((GREEN/NIR) ** 2) - 1422.7 * (GREEN/NIR) +665.7', {
      'GREEN': image.select('B3'),
      'NIR': image.select('B8')

    }).rename('SST');

  return image.addBands(SST);
};

var SST_IC = S2Clipped.map(addSST);

https://code.earthengine.google.com/12958271968ce90c83898ed20812276d

1
  • Thank you very much for the explanation and advice, the code runs perfectly and now I understand better! Apr 18, 2022 at 21:11

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.