update to ghs urban layer

code/end_to_end_gee/01-Earth-Engine-Basics/06a_Import

@@ -1,6 +1,7 @@
 // Let's import some data to Earth Engine
 
-// Download the 'Urban Areas' from Natural Earth
-// https://www.naturalearthdata.com/downloads/10m-cultural-vectors/
+// We will upload the 'GHS Urban Centers' database from JRC
+// https://ghsl.jrc.ec.europa.eu/ghs_stat_ucdb2015mt_r2019a.php
 
-// Unzip and upload the ne_10m_urban_areas shapefile
+// Download the shapefile from https://bit.ly/ghs-ucdb-shapefile
+// Unzip and upload

code/end_to_end_gee/01-Earth-Engine-Basics/06b_Import_(complete)

@@ -1,8 +1,13 @@
 // Let's import some data to Earth Engine
-// Upload the ne_10m_urban_areas.shp shapefile
+
+// Upload the 'GHS Urban Centers' database from JRC
+// https://ghsl.jrc.ec.europa.eu/ghs_stat_ucdb2015mt_r2019a.php
+
+// Download the shapefile from https://bit.ly/ghs-ucdb-shapefile
+// Unzip and upload
 
 // Import the collection
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
 // Visualize the collection
 Map.addLayer(urban, {color: 'blue'}, 'Urban Areas');

code/end_to_end_gee/01-Earth-Engine-Basics/06c_Import_(exercise)

@@ -1,5 +1,10 @@
-// Apply a filter to select only large urban areas
-// Use the property 'area_sqkm' and select all features 
-// with area > 400 sq.km
+// Apply a filter to select only large urban centers
+// in your country and display it on the map.
 
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+// Select all urban centers in your country with
+// a population greater than 1000000
+
+// Use the property 'CTR_MN_NM' containing country names
+// Use the property 'P15' containing 2015 Population
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
+print(urban.first())

code/end_to_end_gee/01-Earth-Engine-Basics/07a_Clipping

@@ -1,7 +1,7 @@
 /**** Start of imports. If edited, may not auto-convert in the playground. ****/
 var geometry = /* color: #d63000 */ee.Geometry.Point([77.5732612484775, 12.974168627832677]),
-    urban = ee.FeatureCollection("users/ujavalgandhi/e2e/ne_10m_urban_areas"),
-    s2 = ee.ImageCollection("COPERNICUS/S2_HARMONIZED");
+    s2 = ee.ImageCollection("COPERNICUS/S2_HARMONIZED"),
+    urban = ee.FeatureCollection("users/ujavalgandhi/e2e/ghs_urban_centers");
 /***** End of imports. If edited, may not auto-convert in the playground. *****/
 var rgbVis = {
   min: 0.0,

code/end_to_end_gee/01-Earth-Engine-Basics/07b_Clipping_(complete)

@@ -1,8 +1,8 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection("users/ujavalgandhi/e2e/ne_10m_urban_areas");
-
-// Find the feature id by adding the layer to the map and using Inspector.
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
+// Find the name of the urban centre
+// by adding the layer to the map and using Inspector.
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 
 var geometry = filtered.geometry();
 
@@ -19,5 +19,6 @@ var image = filtered.median();
 
 var clipped = image.clip(geometry);
 
+Map.centerObject(geometry);
 Map.addLayer(clipped, rgbVis, 'Clipped');
 

code/end_to_end_gee/01-Earth-Engine-Basics/07c_Clipping_(exercise)

@@ -1,8 +1,11 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
 // Change the filter to your home city or any urban area of your choice
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+// Find the name of the urban centre
+// by adding the layer to the map and using Inspector.
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
+
 var geometry = filtered.geometry();
 
 var rgbVis = {

code/end_to_end_gee/01-Earth-Engine-Basics/08a_Export

@@ -1,7 +1,8 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
+
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
 var geometry = filtered.geometry();
 
 var rgbVis = {

code/end_to_end_gee/01-Earth-Engine-Basics/08b_Export_(complete)

@@ -1,7 +1,7 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 var geometry = filtered.geometry();
 
 var rgbVis = {
@@ -17,6 +17,7 @@ var image = filtered.median();
 
 var clipped = image.clip(geometry);
 
+Map.centerObject(geometry);
 Map.addLayer(clipped, rgbVis, 'Clipped');
 
 var exportImage = clipped.select('B.*');

code/end_to_end_gee/01-Earth-Engine-Basics/08c_Export_(exercise)

@@ -2,9 +2,9 @@
 // Write the export function to export the results
 
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 var geometry = filtered.geometry();
 
 var rgbVis = {
@@ -19,6 +19,7 @@ var filtered = s2.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
 var image = filtered.median(); 
 
 var clipped = image.clip(geometry);
+Map.centerObject(geometry);
 Map.addLayer(clipped, rgbVis, 'Clipped');
 
 var exportImage = clipped.select('B.*');

code/end_to_end_gee/03-Supervised-Classification/01c_Basic_Supervised_Classification_(exercise)

@@ -1,9 +1,12 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED');
-var urbanAreas = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urbanAreas = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
 // Perform supervised classification for your city
-// Find the feature id by adding the layer to the map and using Inspector.
-var city = urbanAreas.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+
+// Find the name of the urban centre
+// by adding the layer to the map and using Inspector.
+var city = urbanAreas.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
+
 var geometry = city.geometry();
 Map.centerObject(geometry);
 

code/end_to_end_gee/05-Earth-Engine-Apps/05a_Split_Panel_App

@@ -29,6 +29,7 @@ var scaleImage = function(image) {
 var filtered = s2
   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
   .filter(ee.Filter.bounds(geometry))
+  .filter(ee.Filter.date('2020-01-01', '2021-01-01'))
   .map(maskS2clouds)
   .map(scaleImage);
 

code/end_to_end_gee/05-Earth-Engine-Apps/05b_Split_Panel_App_(complete)

@@ -29,6 +29,7 @@ var scaleImage = function(image) {
 var filtered = s2
   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
   .filter(ee.Filter.bounds(geometry))
+  .filter(ee.Filter.date('2020-01-01', '2021-01-01'))
   .map(maskS2clouds)
   .map(scaleImage);
 

code/end_to_end_gee/05-Earth-Engine-Apps/05c_Split_Panel_App_(exercise)

@@ -29,6 +29,7 @@ var scaleImage = function(image) {
 var filtered = s2
   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
   .filter(ee.Filter.bounds(geometry))
+  .filter(ee.Filter.date('2020-01-01', '2021-01-01'))
   .map(maskS2clouds)
   .map(scaleImage);
 

code/end_to_end_gee/Supplement/Change_Detection/Gaza_Conflict_Change_Detection_Index

@@ -43,13 +43,21 @@ var nrgVis = {min: 0, max: 0.3, bands: ['B8', 'B4', 'B3']};
 Map.addLayer(beforeComposite.clip(geometry), {min: 0, max: 0.3, bands: ['B4', 'B3', 'B2']}, 'Before')
 Map.addLayer(afterComposite.clip(geometry), {min: 0, max: 0.4, bands: ['B4', 'B3', 'B2']}, 'After')
 
-var ndviBefore = beforeComposite.normalizedDifference(['B8', 'B4']).rename(['ndvi']);
-var ndviAfter = afterComposite.normalizedDifference(['B8', 'B4']).rename(['ndvi']);
+var addIndices = function(image) {
+  var nbr = image.normalizedDifference(['B8', 'B12']).rename(['nbr']);
+  var ndvi = image.normalizedDifference(['B8', 'B4']).rename(['ndvi'])
+  return image.addBands(nbr).addBands(ndvi);
+}
+
+// You can try change detection with either NBR or NDVI
+var selectedIndex = 'nbr';
+var beforeNbr = addIndices(beforeComposite).select(selectedIndex);
+var afterNbr = addIndices(afterComposite).select(selectedIndex);
 
-// Loss of NDVI
-var difference = ndviBefore.subtract(ndviAfter).clip(geometry)
+// Change in Index
+var difference = beforeNbr.subtract(afterNbr).clip(geometry)
 
-var indexThreshold = 0.12;
+var indexThreshold = 0.25;
 var change = difference.gt(indexThreshold)
 
 // Mask Waterbodies using WorldCover

code/end_to_end_gee/Supplement/Supervised_Classification/Classification_with_Migrated_Training_Samples

@@ -51,7 +51,7 @@ Map.addLayer(composite2020, rgbVis, 'Composite 2020');
 
 
 // Compute Spectral Distance between 2019 and 2020 images
-var distance = composite2019.spectralDistance(composite2020);
+var distance = composite2019.spectralDistance(composite2020, 'sam');
 // GCPs were collected on 2020 image
 // Find out which GCPs are still unchanged in 2021
 

code/end_to_end_gee/Supplement/Supervised_Classification/Time_Series_Modeling

@@ -203,7 +203,6 @@ var amplitude = harmonicTrendCoefficients.select('sin_1')
 var rgb = ee.Image.cat([phase, amplitude, ee.Image(1)]).hsvToRgb();
 Map.addLayer(rgb, {}, 'Phase (hue), Amplitude (sat)', false);
 
-
 // The Phase and Amplitude values will be very different
 // at farms following different cropping cycles
 // Let's plot and compare the fitted time series

docs/code/end_to_end_gee/01-Earth-Engine-Basics/06a_Import

@@ -1,6 +1,7 @@
 // Let's import some data to Earth Engine
 
-// Download the 'Urban Areas' from Natural Earth
-// https://www.naturalearthdata.com/downloads/10m-cultural-vectors/
+// We will upload the 'GHS Urban Centers' database from JRC
+// https://ghsl.jrc.ec.europa.eu/ghs_stat_ucdb2015mt_r2019a.php
 
-// Unzip and upload the ne_10m_urban_areas shapefile
+// Download the shapefile from https://bit.ly/ghs-ucdb-shapefile
+// Unzip and upload

docs/code/end_to_end_gee/01-Earth-Engine-Basics/06b_Import_(complete)

@@ -1,8 +1,13 @@
 // Let's import some data to Earth Engine
-// Upload the ne_10m_urban_areas.shp shapefile
+
+// Upload the 'GHS Urban Centers' database from JRC
+// https://ghsl.jrc.ec.europa.eu/ghs_stat_ucdb2015mt_r2019a.php
+
+// Download the shapefile from https://bit.ly/ghs-ucdb-shapefile
+// Unzip and upload
 
 // Import the collection
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
 // Visualize the collection
 Map.addLayer(urban, {color: 'blue'}, 'Urban Areas');

docs/code/end_to_end_gee/01-Earth-Engine-Basics/06c_Import_(exercise)

@@ -1,5 +1,10 @@
-// Apply a filter to select only large urban areas
-// Use the property 'area_sqkm' and select all features 
-// with area > 400 sq.km
+// Apply a filter to select only large urban centers
+// in your country and display it on the map.
 
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+// Select all urban centers in your country with
+// a population greater than 1000000
+
+// Use the property 'CTR_MN_NM' containing country names
+// Use the property 'P15' containing 2015 Population
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
+print(urban.first())

docs/code/end_to_end_gee/01-Earth-Engine-Basics/07a_Clipping

@@ -1,7 +1,7 @@
 /**** Start of imports. If edited, may not auto-convert in the playground. ****/
 var geometry = /* color: #d63000 */ee.Geometry.Point([77.5732612484775, 12.974168627832677]),
-    urban = ee.FeatureCollection("users/ujavalgandhi/e2e/ne_10m_urban_areas"),
-    s2 = ee.ImageCollection("COPERNICUS/S2_HARMONIZED");
+    s2 = ee.ImageCollection("COPERNICUS/S2_HARMONIZED"),
+    urban = ee.FeatureCollection("users/ujavalgandhi/e2e/ghs_urban_centers");
 /***** End of imports. If edited, may not auto-convert in the playground. *****/
 var rgbVis = {
   min: 0.0,

docs/code/end_to_end_gee/01-Earth-Engine-Basics/07b_Clipping_(complete)

@@ -1,8 +1,8 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection("users/ujavalgandhi/e2e/ne_10m_urban_areas");
-
-// Find the feature id by adding the layer to the map and using Inspector.
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
+// Find the name of the urban centre
+// by adding the layer to the map and using Inspector.
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 
 var geometry = filtered.geometry();
 
@@ -19,5 +19,6 @@ var image = filtered.median();
 
 var clipped = image.clip(geometry);
 
+Map.centerObject(geometry);
 Map.addLayer(clipped, rgbVis, 'Clipped');
 

docs/code/end_to_end_gee/01-Earth-Engine-Basics/07c_Clipping_(exercise)

@@ -1,8 +1,11 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
 // Change the filter to your home city or any urban area of your choice
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+// Find the name of the urban centre
+// by adding the layer to the map and using Inspector.
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
+
 var geometry = filtered.geometry();
 
 var rgbVis = {

docs/code/end_to_end_gee/01-Earth-Engine-Basics/08a_Export

@@ -1,7 +1,8 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
+
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
 var geometry = filtered.geometry();
 
 var rgbVis = {

docs/code/end_to_end_gee/01-Earth-Engine-Basics/08b_Export_(complete)

@@ -1,7 +1,7 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 var geometry = filtered.geometry();
 
 var rgbVis = {
@@ -17,6 +17,7 @@ var image = filtered.median();
 
 var clipped = image.clip(geometry);
 
+Map.centerObject(geometry);
 Map.addLayer(clipped, rgbVis, 'Clipped');
 
 var exportImage = clipped.select('B.*');

docs/code/end_to_end_gee/01-Earth-Engine-Basics/08c_Export_(exercise)

@@ -2,9 +2,9 @@
 // Write the export function to export the results
 
 var s2 = ee.ImageCollection('COPERNICUS/S2_HARMONIZED');
-var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urban = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
-var filtered = urban.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+var filtered = urban.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
 var geometry = filtered.geometry();
 
 var rgbVis = {
@@ -19,6 +19,7 @@ var filtered = s2.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
 var image = filtered.median(); 
 
 var clipped = image.clip(geometry);
+Map.centerObject(geometry);
 Map.addLayer(clipped, rgbVis, 'Clipped');
 
 var exportImage = clipped.select('B.*');

docs/code/end_to_end_gee/03-Supervised-Classification/01c_Basic_Supervised_Classification_(exercise)

@@ -1,9 +1,12 @@
 var s2 = ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED');
-var urbanAreas = ee.FeatureCollection('users/ujavalgandhi/e2e/ne_10m_urban_areas');
+var urbanAreas = ee.FeatureCollection('users/ujavalgandhi/e2e/ghs_urban_centers');
 
 // Perform supervised classification for your city
-// Find the feature id by adding the layer to the map and using Inspector.
-var city = urbanAreas.filter(ee.Filter.eq('system:index', '00000000000000002bf8'));
+
+// Find the name of the urban centre
+// by adding the layer to the map and using Inspector.
+var city = urbanAreas.filter(ee.Filter.eq('UC_NM_MN', 'Bengaluru'));
+
 var geometry = city.geometry();
 Map.centerObject(geometry);
 

docs/code/end_to_end_gee/05-Earth-Engine-Apps/05a_Split_Panel_App

@@ -29,6 +29,7 @@ var scaleImage = function(image) {
 var filtered = s2
   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
   .filter(ee.Filter.bounds(geometry))
+  .filter(ee.Filter.date('2020-01-01', '2021-01-01'))
   .map(maskS2clouds)
   .map(scaleImage);
 

docs/code/end_to_end_gee/05-Earth-Engine-Apps/05b_Split_Panel_App_(complete)

@@ -29,6 +29,7 @@ var scaleImage = function(image) {
 var filtered = s2
   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
   .filter(ee.Filter.bounds(geometry))
+  .filter(ee.Filter.date('2020-01-01', '2021-01-01'))
   .map(maskS2clouds)
   .map(scaleImage);
 

docs/code/end_to_end_gee/05-Earth-Engine-Apps/05c_Split_Panel_App_(exercise)

@@ -29,6 +29,7 @@ var scaleImage = function(image) {
 var filtered = s2
   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 30))
   .filter(ee.Filter.bounds(geometry))
+  .filter(ee.Filter.date('2020-01-01', '2021-01-01'))
   .map(maskS2clouds)
   .map(scaleImage);