Merge pull request #347 from azlinszkysinergise/sentinel2_magic_eyes

Agri Growth Stage for Quarterly Mosaics

Author: Guyon-Duifhuizen

sentinel-2/agriculture_growth_stage/README.md

@@ -5,6 +5,9 @@ grand_parent: Sentinel
 layout: script
 permalink: /sentinel-2/agriculture_growth_stage/
 nav_exclude: true
+scripts:
+  - [Copernicus Browser, script.js]
+  - [Sentinel-2 mosaic, quarterly_mosaics.js]
 examples:
 - zoom: '11'
   lat: '45.6246'

sentinel-2/agriculture_growth_stage/quarterly_mosaics.js

@@ -0,0 +1,75 @@
+//VERSION=3 (auto-converted from 1)
+/*
+Source: @HarelDan - https://github.com/hareldunn/GIS_Repo/blob/master/Multi-Temporal%20NDVI%20for%20Sentinel%20Hub%20Custom%20Scripts
+Visualizing NDVI multi-temporal trends in Sentinel-2 imagery.
+will take the current image as baseline and calculate average NDVI for the previous 2 months
+Based on:
+https://twitter.com/sentinel_hub/status/922813457145221121
+https://twitter.com/sentinel_hub/status/1020755996359225344
+Script requires multi-temporal processing so parameter TEMPORAL=true should be added to the request.
+*/
+
+function setup() {
+    return {
+        input: [
+            {
+                bands: ["B04", "B08"],
+            },
+        ],
+        output: { bands: 3 },
+        mosaicking: "ORBIT",
+    };
+}
+
+function calcNDVI(sample) {
+    var denom = sample.B04 + sample.B08;
+    return denom != 0 ? (sample.B08 - sample.B04) / denom : 0.0;
+}
+function stretch(val, min, max) {
+    return (val - min) / (max - min);
+}
+
+function evaluatePixel(samples, scenes) {
+    var avg1 = 0;
+    var count1 = 0;
+    var avg2 = 0;
+    var count2 = 0;
+    var avg3 = 0;
+    var count3 = 0;
+    var endMonth = scenes[0].date.getMonth();
+
+    for (var i = 0; i < samples.length; i++) {
+        var ndvi = calcNDVI(samples[i]);
+        if (scenes[i].date.getMonth() == endMonth) {
+            avg3 = avg3 + ndvi;
+            count3++;
+        } else if (scenes[i].date.getMonth() == endMonth - 1) {
+            avg2 = avg2 + ndvi;
+            count2++;
+        } else {
+            avg1 = avg1 + ndvi;
+            count1++;
+        }
+    }
+    avg1 = avg1 / count1;
+    avg2 = avg2 / count2;
+    avg3 = avg3 / count3;
+    avg1 = stretch(avg1, 0.1, 0.7);
+    avg2 = stretch(avg2, 0.1, 0.7);
+    avg3 = stretch(avg3, 0.1, 0.7);
+
+    return [avg1, avg2, avg3];
+}
+
+function preProcessScenes(collections) {
+    collections.scenes.orbits = collections.scenes.orbits.filter(function (
+        orbit
+    ) {
+        var orbitDateFrom = new Date(orbit.dateFrom);
+        return (
+            orbitDateFrom.getTime() >=
+            collections.to.getTime() - 3 * 31 * 24 * 3600 * 1000
+        );
+    });
+    return collections;
+}

sentinel-5p/sentinel-5p.md

@@ -41,6 +41,6 @@ Sentinel-5P provides atmospheric measurements, relating to air quality, climate
  - [Cloud Optical Thickness](/sentinel-5p/cloud-optical-thickness)
  - [Cloud Effective Radiometric Fraction](/sentinel-5p/cloud-radiometric-fraction)
 
- #### Additional Products
+#### Additional Products
 
  - [Nitrogen Dioxide monthly mean](/sentinel-5p/no2_monthly_mean)