Function to calculate red-edgeNDVI in Google Earth Engine Unicorn Meta Zoo #1: Why another podcast? Announcing the arrival of Valued Associate #679: Cesar ManaraGEE cloud-free Sentinel2 and linear RegressionLinearFit with Google Earth EngineCalculate MSAVI (Modified Soil-adjusted Vegetation Index) in Google Earth EngineGoogle Earth Engine - Map.addLayerGoogle Earth Engine, how to distinguish between rivers/streams and ponds/lakes in a water maskMissing band in mapping function using Google Earth EngineEarth Engine - function for printing multiple histogramsEarth Engine create user defined function with default parametersEarth Engine how to avoid timeout for buffer calculation?Google Earth Engine formaTrend function

Map material from china not allowed to leave the country

A strange hotel

With indentation set to `0em`, when using a line break, there is still an indentation of a size of a space

Israeli soda type drink

Need of separate security plugins for both root and subfolder sites Wordpress?

What is /etc/mtab in Linux?

What is "leading note" and what does it mean to "raise a note"?

How to find the right literary agent in the USA?

Has a Nobel Peace laureate ever been accused of war crimes?

Are all CP/M-80 implementations binary compatible?

Suing a Police Officer Instead of the Police Department

What is ls Largest Number Formed by only moving two sticks in 508?

Will I lose my paid in full property

Raising a bilingual kid. When should we introduce the majority language?

Can I criticise the more senior developers around me for not writing clean code?

How do I check if a string is entirely made of the same substring?

Is a 5 watt UHF/VHF handheld considered QRP?

Identify story/novel: Tribe on colonized planet, not aware of this. "Taboo," altitude sickness, robot guardian (60s? Young Adult?)

Why did C use the -> operator instead of reusing the . operator?

Passing args from the bash script to the function in the script

Mistake in years of experience in resume?

What *exactly* is electrical current, voltage, and resistance?

What do you call the part of a novel that is not dialog?

Align column where each cell has two decimals with siunitx



Function to calculate red-edgeNDVI in Google Earth Engine



Unicorn Meta Zoo #1: Why another podcast?
Announcing the arrival of Valued Associate #679: Cesar ManaraGEE cloud-free Sentinel2 and linear RegressionLinearFit with Google Earth EngineCalculate MSAVI (Modified Soil-adjusted Vegetation Index) in Google Earth EngineGoogle Earth Engine - Map.addLayerGoogle Earth Engine, how to distinguish between rivers/streams and ponds/lakes in a water maskMissing band in mapping function using Google Earth EngineEarth Engine - function for printing multiple histogramsEarth Engine create user defined function with default parametersEarth Engine how to avoid timeout for buffer calculation?Google Earth Engine formaTrend function



.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;








1















I am trying to calculate the red-edgeNDVI ((NIR – red edge)/(NIR + red edge)) of an Image.Collection in Google Earth Engine. There is a built in function for NDVI but not for red-edge NDVI. The code I have is below but I get the error message "image.NIR is undefined". I'm sure this is a simple syntax issue related to the function but I've tried everything I could think of and I just can't get it to work. Any GEE/java experts who can tell me what I'm missing?



/**
* Function to mask clouds using the Sentinel-2 QA band
* @param ee.Image image Sentinel-2 image
* @return ee.Image cloud masked Sentinel-2 image
*/
function maskS2clouds(image)
var qa = image.select('QA60');

// Bits 10 and 11 are clouds and cirrus, respectively.
var cloudBitMask = 1 << 10;
var cirrusBitMask = 1 << 11;

// Both flags should be set to zero, indicating clear conditions.
var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
.and(qa.bitwiseAnd(cirrusBitMask).eq(0));

return image.updateMask(mask)
.divide(10000);


// Load Sentinel-2 TOA reflectance data.
var S2 = ee.ImageCollection('COPERNICUS/S2')
.filterDate('2017-06-01', '2017-09-30')
// Pre-filter to get less cloudy granules.
.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
//Select required bands only
.select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
//Apply cloud mask
.map(maskS2clouds);

//Create band variables
var redEdge = S2.select('B5');
var NIR = S2.select('B8');

//Function to calculate redEdgeNDVI
var add_reNDVI = function(image)
var redEdgeNDVI = image.NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
return image.addBands(reNDVI);
;









share|improve this question




























    1















    I am trying to calculate the red-edgeNDVI ((NIR – red edge)/(NIR + red edge)) of an Image.Collection in Google Earth Engine. There is a built in function for NDVI but not for red-edge NDVI. The code I have is below but I get the error message "image.NIR is undefined". I'm sure this is a simple syntax issue related to the function but I've tried everything I could think of and I just can't get it to work. Any GEE/java experts who can tell me what I'm missing?



    /**
    * Function to mask clouds using the Sentinel-2 QA band
    * @param ee.Image image Sentinel-2 image
    * @return ee.Image cloud masked Sentinel-2 image
    */
    function maskS2clouds(image)
    var qa = image.select('QA60');

    // Bits 10 and 11 are clouds and cirrus, respectively.
    var cloudBitMask = 1 << 10;
    var cirrusBitMask = 1 << 11;

    // Both flags should be set to zero, indicating clear conditions.
    var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
    .and(qa.bitwiseAnd(cirrusBitMask).eq(0));

    return image.updateMask(mask)
    .divide(10000);


    // Load Sentinel-2 TOA reflectance data.
    var S2 = ee.ImageCollection('COPERNICUS/S2')
    .filterDate('2017-06-01', '2017-09-30')
    // Pre-filter to get less cloudy granules.
    .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
    //Select required bands only
    .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
    //Apply cloud mask
    .map(maskS2clouds);

    //Create band variables
    var redEdge = S2.select('B5');
    var NIR = S2.select('B8');

    //Function to calculate redEdgeNDVI
    var add_reNDVI = function(image)
    var redEdgeNDVI = image.NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
    return image.addBands(reNDVI);
    ;









    share|improve this question
























      1












      1








      1








      I am trying to calculate the red-edgeNDVI ((NIR – red edge)/(NIR + red edge)) of an Image.Collection in Google Earth Engine. There is a built in function for NDVI but not for red-edge NDVI. The code I have is below but I get the error message "image.NIR is undefined". I'm sure this is a simple syntax issue related to the function but I've tried everything I could think of and I just can't get it to work. Any GEE/java experts who can tell me what I'm missing?



      /**
      * Function to mask clouds using the Sentinel-2 QA band
      * @param ee.Image image Sentinel-2 image
      * @return ee.Image cloud masked Sentinel-2 image
      */
      function maskS2clouds(image)
      var qa = image.select('QA60');

      // Bits 10 and 11 are clouds and cirrus, respectively.
      var cloudBitMask = 1 << 10;
      var cirrusBitMask = 1 << 11;

      // Both flags should be set to zero, indicating clear conditions.
      var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
      .and(qa.bitwiseAnd(cirrusBitMask).eq(0));

      return image.updateMask(mask)
      .divide(10000);


      // Load Sentinel-2 TOA reflectance data.
      var S2 = ee.ImageCollection('COPERNICUS/S2')
      .filterDate('2017-06-01', '2017-09-30')
      // Pre-filter to get less cloudy granules.
      .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
      //Select required bands only
      .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
      //Apply cloud mask
      .map(maskS2clouds);

      //Create band variables
      var redEdge = S2.select('B5');
      var NIR = S2.select('B8');

      //Function to calculate redEdgeNDVI
      var add_reNDVI = function(image)
      var redEdgeNDVI = image.NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
      return image.addBands(reNDVI);
      ;









      share|improve this question














      I am trying to calculate the red-edgeNDVI ((NIR – red edge)/(NIR + red edge)) of an Image.Collection in Google Earth Engine. There is a built in function for NDVI but not for red-edge NDVI. The code I have is below but I get the error message "image.NIR is undefined". I'm sure this is a simple syntax issue related to the function but I've tried everything I could think of and I just can't get it to work. Any GEE/java experts who can tell me what I'm missing?



      /**
      * Function to mask clouds using the Sentinel-2 QA band
      * @param ee.Image image Sentinel-2 image
      * @return ee.Image cloud masked Sentinel-2 image
      */
      function maskS2clouds(image)
      var qa = image.select('QA60');

      // Bits 10 and 11 are clouds and cirrus, respectively.
      var cloudBitMask = 1 << 10;
      var cirrusBitMask = 1 << 11;

      // Both flags should be set to zero, indicating clear conditions.
      var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
      .and(qa.bitwiseAnd(cirrusBitMask).eq(0));

      return image.updateMask(mask)
      .divide(10000);


      // Load Sentinel-2 TOA reflectance data.
      var S2 = ee.ImageCollection('COPERNICUS/S2')
      .filterDate('2017-06-01', '2017-09-30')
      // Pre-filter to get less cloudy granules.
      .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
      //Select required bands only
      .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
      //Apply cloud mask
      .map(maskS2clouds);

      //Create band variables
      var redEdge = S2.select('B5');
      var NIR = S2.select('B8');

      //Function to calculate redEdgeNDVI
      var add_reNDVI = function(image)
      var redEdgeNDVI = image.NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
      return image.addBands(reNDVI);
      ;






      google-earth-engine function vegetation-index






      share|improve this question













      share|improve this question











      share|improve this question




      share|improve this question










      asked 4 hours ago









      Adam GAdam G

      262




      262




















          1 Answer
          1






          active

          oldest

          votes


















          3














          In your code, S2 is an ImageCollection, so when you "create band variables" you're just getting ImageCollections in which every image inside has only the selected band, which is useful. As you well commented, add_reNDVI is a function that will take every image in the collection and calculate reNDVI. So you have to map that function over the collection to get what you want.



          // Load Sentinel-2 TOA reflectance data.
          var S2 = ee.ImageCollection('COPERNICUS/S2')
          .filterDate('2017-06-01', '2017-09-30')
          // Pre-filter to get less cloudy granules.
          .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
          //Select required bands only
          .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
          //Apply cloud mask
          .map(maskS2clouds);

          //Function to calculate redEdgeNDVI
          var add_reNDVI = function(image)
          //Create band variables
          var redEdge = image.select('B5');
          var NIR = image.select('B8');
          var redEdgeNDVI = NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
          return image.addBands(redEdgeNDVI);
          ;
          var s2_reNDVI = S2.map(add_reNDVI)
          Map.centerObject(s2_reNDVI.first())
          Map.addLayer(s2_reNDVI.first())


          I guess you'll need to filter by bounds at some point, that is why I just took the first image out of the collection to check if it's working






          share|improve this answer























            Your Answer








            StackExchange.ready(function()
            var channelOptions =
            tags: "".split(" "),
            id: "79"
            ;
            initTagRenderer("".split(" "), "".split(" "), channelOptions);

            StackExchange.using("externalEditor", function()
            // Have to fire editor after snippets, if snippets enabled
            if (StackExchange.settings.snippets.snippetsEnabled)
            StackExchange.using("snippets", function()
            createEditor();
            );

            else
            createEditor();

            );

            function createEditor()
            StackExchange.prepareEditor(
            heartbeatType: 'answer',
            autoActivateHeartbeat: false,
            convertImagesToLinks: false,
            noModals: true,
            showLowRepImageUploadWarning: true,
            reputationToPostImages: null,
            bindNavPrevention: true,
            postfix: "",
            imageUploader:
            brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
            contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
            allowUrls: true
            ,
            onDemand: true,
            discardSelector: ".discard-answer"
            ,immediatelyShowMarkdownHelp:true
            );



            );













            draft saved

            draft discarded


















            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f320766%2ffunction-to-calculate-red-edgendvi-in-google-earth-engine%23new-answer', 'question_page');

            );

            Post as a guest















            Required, but never shown

























            1 Answer
            1






            active

            oldest

            votes








            1 Answer
            1






            active

            oldest

            votes









            active

            oldest

            votes






            active

            oldest

            votes









            3














            In your code, S2 is an ImageCollection, so when you "create band variables" you're just getting ImageCollections in which every image inside has only the selected band, which is useful. As you well commented, add_reNDVI is a function that will take every image in the collection and calculate reNDVI. So you have to map that function over the collection to get what you want.



            // Load Sentinel-2 TOA reflectance data.
            var S2 = ee.ImageCollection('COPERNICUS/S2')
            .filterDate('2017-06-01', '2017-09-30')
            // Pre-filter to get less cloudy granules.
            .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
            //Select required bands only
            .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
            //Apply cloud mask
            .map(maskS2clouds);

            //Function to calculate redEdgeNDVI
            var add_reNDVI = function(image)
            //Create band variables
            var redEdge = image.select('B5');
            var NIR = image.select('B8');
            var redEdgeNDVI = NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
            return image.addBands(redEdgeNDVI);
            ;
            var s2_reNDVI = S2.map(add_reNDVI)
            Map.centerObject(s2_reNDVI.first())
            Map.addLayer(s2_reNDVI.first())


            I guess you'll need to filter by bounds at some point, that is why I just took the first image out of the collection to check if it's working






            share|improve this answer



























              3














              In your code, S2 is an ImageCollection, so when you "create band variables" you're just getting ImageCollections in which every image inside has only the selected band, which is useful. As you well commented, add_reNDVI is a function that will take every image in the collection and calculate reNDVI. So you have to map that function over the collection to get what you want.



              // Load Sentinel-2 TOA reflectance data.
              var S2 = ee.ImageCollection('COPERNICUS/S2')
              .filterDate('2017-06-01', '2017-09-30')
              // Pre-filter to get less cloudy granules.
              .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
              //Select required bands only
              .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
              //Apply cloud mask
              .map(maskS2clouds);

              //Function to calculate redEdgeNDVI
              var add_reNDVI = function(image)
              //Create band variables
              var redEdge = image.select('B5');
              var NIR = image.select('B8');
              var redEdgeNDVI = NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
              return image.addBands(redEdgeNDVI);
              ;
              var s2_reNDVI = S2.map(add_reNDVI)
              Map.centerObject(s2_reNDVI.first())
              Map.addLayer(s2_reNDVI.first())


              I guess you'll need to filter by bounds at some point, that is why I just took the first image out of the collection to check if it's working






              share|improve this answer

























                3












                3








                3







                In your code, S2 is an ImageCollection, so when you "create band variables" you're just getting ImageCollections in which every image inside has only the selected band, which is useful. As you well commented, add_reNDVI is a function that will take every image in the collection and calculate reNDVI. So you have to map that function over the collection to get what you want.



                // Load Sentinel-2 TOA reflectance data.
                var S2 = ee.ImageCollection('COPERNICUS/S2')
                .filterDate('2017-06-01', '2017-09-30')
                // Pre-filter to get less cloudy granules.
                .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
                //Select required bands only
                .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
                //Apply cloud mask
                .map(maskS2clouds);

                //Function to calculate redEdgeNDVI
                var add_reNDVI = function(image)
                //Create band variables
                var redEdge = image.select('B5');
                var NIR = image.select('B8');
                var redEdgeNDVI = NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
                return image.addBands(redEdgeNDVI);
                ;
                var s2_reNDVI = S2.map(add_reNDVI)
                Map.centerObject(s2_reNDVI.first())
                Map.addLayer(s2_reNDVI.first())


                I guess you'll need to filter by bounds at some point, that is why I just took the first image out of the collection to check if it's working






                share|improve this answer













                In your code, S2 is an ImageCollection, so when you "create band variables" you're just getting ImageCollections in which every image inside has only the selected band, which is useful. As you well commented, add_reNDVI is a function that will take every image in the collection and calculate reNDVI. So you have to map that function over the collection to get what you want.



                // Load Sentinel-2 TOA reflectance data.
                var S2 = ee.ImageCollection('COPERNICUS/S2')
                .filterDate('2017-06-01', '2017-09-30')
                // Pre-filter to get less cloudy granules.
                .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
                //Select required bands only
                .select('B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'QA60')
                //Apply cloud mask
                .map(maskS2clouds);

                //Function to calculate redEdgeNDVI
                var add_reNDVI = function(image)
                //Create band variables
                var redEdge = image.select('B5');
                var NIR = image.select('B8');
                var redEdgeNDVI = NIR.subtract(redEdge).divide(NIR.add(redEdge)).rename('reNDVI');
                return image.addBands(redEdgeNDVI);
                ;
                var s2_reNDVI = S2.map(add_reNDVI)
                Map.centerObject(s2_reNDVI.first())
                Map.addLayer(s2_reNDVI.first())


                I guess you'll need to filter by bounds at some point, that is why I just took the first image out of the collection to check if it's working







                share|improve this answer












                share|improve this answer



                share|improve this answer










                answered 3 hours ago









                Rodrigo E. PrincipeRodrigo E. Principe

                4,54111021




                4,54111021



























                    draft saved

                    draft discarded
















































                    Thanks for contributing an answer to Geographic Information Systems Stack Exchange!


                    • Please be sure to answer the question. Provide details and share your research!

                    But avoid


                    • Asking for help, clarification, or responding to other answers.

                    • Making statements based on opinion; back them up with references or personal experience.

                    To learn more, see our tips on writing great answers.




                    draft saved


                    draft discarded














                    StackExchange.ready(
                    function ()
                    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f320766%2ffunction-to-calculate-red-edgendvi-in-google-earth-engine%23new-answer', 'question_page');

                    );

                    Post as a guest















                    Required, but never shown





















































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown

































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown







                    Popular posts from this blog

                    Oświęcim Innehåll Historia | Källor | Externa länkar | Navigeringsmeny50°2′18″N 19°13′17″Ö / 50.03833°N 19.22139°Ö / 50.03833; 19.2213950°2′18″N 19°13′17″Ö / 50.03833°N 19.22139°Ö / 50.03833; 19.221393089658Nordisk familjebok, AuschwitzInsidan tro och existensJewish Community i OświęcimAuschwitz Jewish Center: MuseumAuschwitz Jewish Center

                    Valle di Casies Indice Geografia fisica | Origini del nome | Storia | Società | Amministrazione | Sport | Note | Bibliografia | Voci correlate | Altri progetti | Collegamenti esterni | Menu di navigazione46°46′N 12°11′E / 46.766667°N 12.183333°E46.766667; 12.183333 (Valle di Casies)46°46′N 12°11′E / 46.766667°N 12.183333°E46.766667; 12.183333 (Valle di Casies)Sito istituzionaleAstat Censimento della popolazione 2011 - Determinazione della consistenza dei tre gruppi linguistici della Provincia Autonoma di Bolzano-Alto Adige - giugno 2012Numeri e fattiValle di CasiesDato IstatTabella dei gradi/giorno dei Comuni italiani raggruppati per Regione e Provincia26 agosto 1993, n. 412Heraldry of the World: GsiesStatistiche I.StatValCasies.comWikimedia CommonsWikimedia CommonsValle di CasiesSito ufficialeValle di CasiesMM14870458910042978-6

                    Typsetting diagram chases (with TikZ?) Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)How to define the default vertical distance between nodes?Draw edge on arcNumerical conditional within tikz keys?TikZ: Drawing an arc from an intersection to an intersectionDrawing rectilinear curves in Tikz, aka an Etch-a-Sketch drawingLine up nested tikz enviroments or how to get rid of themHow to place nodes in an absolute coordinate system in tikzCommutative diagram with curve connecting between nodesTikz with standalone: pinning tikz coordinates to page cmDrawing a Decision Diagram with Tikz and layout manager