Urban Heat Island Mapping with Landsat and DeepSeek

// Define the region of interest (Fez, Morocco)
var fez = ee.Geometry.Point([ -5.0069, 34.0372]);

// Define the date range
var startDate = '2023-07-01';
var endDate = '2023-08-31';

// Load Landsat 8 SR data
var landsat8 = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR')
  .filterBounds(fez)
  .filterDate(startDate, endDate)
  .filter(ee.Filter.lt('CLOUD_COVER', 20));

// Function to calculate NDVI
var calculateNDVI = function(image) {
  return image.normalizedDifference(['B5', 'B4']).rename('NDVI');
};

// Apply NDVI calculation to the Landsat 8 collection
var ndviCollection = landsat8.map(calculateNDVI);

// Composite the NDVI collection (e.g., median)
var ndvi = ndviCollection.median();

//Display NDVI
Map.addLayer(ndvi, {min: -1, max: 1, palette: ['blue', 'white', 'green']}, 'NDVI');

// Function to calculate LST (simplified)
var calculateLST = function(image) {
  var thermalBand = image.select('B10');
  var kelvin = thermalBand.multiply(0.1).add(273.15); // Convert to Kelvin
  return kelvin.rename('LST');
};

// Apply LST calculation to the Landsat 8 collection
var lstCollection = landsat8.map(calculateLST);

// Composite the LST collection (e.g., median)
var lst = lstCollection.median();

//Display LST
Map.addLayer(lst, {min: 290, max: 320, palette: ['blue', 'yellow', 'red']}, 'LST');

// Define a threshold for UHI (e.g., LST > 308 K)
var uhiThreshold = 308;

// Create a UHI mask
var uhiMask = lst.gt(uhiThreshold);

// Display the UHI mask
Map.addLayer(uhiMask.updateMask(uhiMask), {palette: 'red'}, 'UHI Mask');

// Export the LST image to Google Drive
Export.image.toDrive({
  image: lst,
  description: 'LST_Fez',
  folder: 'GEE_Exports',
  scale: 30,
  region: fez,
  maxPixels: 1e9
});