Scribd
Upload
13 views3 pages

Define The Custom Polygon For The

The document outlines a script for processing Landsat 8 imagery to analyze a custom Area of Interest (AOI) in Bangalore. It includes steps for applying scaling factors, masking clouds, calculating NDVI, and deriving Land Surface Temperature (LST) using emissivity. The processed imagery and results are visualized on a map with specified parameters.

Uploaded by

Ashwin Adve
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
13 views3 pages

Define The Custom Polygon For The

The document outlines a script for processing Landsat 8 imagery to analyze a custom Area of Interest (AOI) in Bangalore. It includes steps for applying scaling factors, masking clouds, calculating NDVI, and deriving Land Surface Temperature (LST) using emissivity. The processed imagery and results are visualized on a map with specified parameters.

Uploaded by

Ashwin Adve
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 3

// Define the custom polygon for the new AOI using the provided coordinates

var customPolygon = ee.Geometry.Polygon([


[
[77.61055588885255, 12.973863092132746],
[77.62995362444825, 12.973863092132746],
[77.62995362444825, 12.989921386728478],
[77.61055588885255, 12.989921386728478],
[77.61055588885255, 12.973863092132746]
]
]);

// Set the AOI to the custom polygon


var aoi = customPolygon;

// Add the AOI to the map


Map.addLayer(aoi, {color: 'red'}, 'Custom AOI - Precise Polygon');

// Center the map on the AOI


Map.centerObject(aoi, 14);

var startDate = '2024-04-20';


var endDate = '2024-05-10' ;
// Applies scaling factors.
function applyScaleFactors(image) {
// Scale and offset values for optical bands
var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);

// Scale and offset values for thermal bands


var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0);

// Add scaled bands to the original image


return image.addBands(opticalBands, null, true)
.addBands(thermalBands, null, true);
}

// Function to Mask Clouds and Cloud Shadows in Landsat 8 Imagery

function cloudMask(image) {
// Define cloud shadow and cloud bitmasks (Bits 3 and 5)
var cloudShadowBitmask = (1 << 3);
var cloudBitmask = (1 << 5);

// Select the Quality Assessment (QA) band for pixel quality information
var qa = image.select('QA_PIXEL');

// Create a binary mask to identify clear conditions (both cloud and cloud shadow
bits set to 0)
var mask = qa.bitwiseAnd(cloudShadowBitmask).eq(0)
.and(qa.bitwiseAnd(cloudBitmask).eq(0));

// Update the original image, masking out cloud and cloud shadow-affected pixels
return image.updateMask(mask);
}

// Import and preprocess Landsat 8 imagery


var image = ee.ImageCollection("LANDSAT/LC08/C02/T1_L2")
.filterBounds(aoi)
.filterDate(startDate, endDate)
.map(applyScaleFactors)
.map(cloudMask)
.median()
.clip(aoi);

// Define visualization parameters for True Color imagery (bands 4, 3, and 2)


var visualization = {
bands: ['SR_B4', 'SR_B3', 'SR_B2'],
min: 0.0,
max: 0.15,
};

// Add the processed image to the map with the specified visualization
Map.addLayer(image, visualization, 'True Color 432');

// Calculate Normalized Difference Vegetation Index (NDVI)


var ndvi = image.normalizedDifference(['SR_B5', 'SR_B4']).rename('NDVI');

// Define NDVI Visualization Parameters


var ndviPalette = {
min: -1,
max: 1,
palette: ['blue', 'white', 'green']
};

Map.addLayer(ndvi, ndviPalette, 'NDVI Bangalore')


// Calculate the minimum NDVI value within the AOI
var ndviMin = ee.Number(ndvi.reduceRegion({
reducer : ee.Reducer.min(),
geometry : aoi,
scale : 30,
maxPixels : 1e9
}).values().get(0));

// Calculate the maximum NDVI value within the AOI


var ndviMax = ee.Number(ndvi.reduceRegion({
reducer : ee.Reducer.max(),
geometry : aoi,
scale : 30,
maxPixels : 1e9
}).values().get(0));
// Fraction of Vegetation (FV) Calculation
// Formula: ((NDVI - NDVI_min) / (NDVI_max - NDVI_min))^2

// Calculate the Fraction of Vegetation (FV) using the NDVI values within the
specified range.
// NDVI_min represents the minimum NDVI value, and NDVI_max represents the maximum
NDVI value
var fv = ((ndvi.subtract(ndviMin)).divide(ndviMax.subtract(ndviMin)))
.pow(ee.Number(2))
.rename('FV');

// Emissivity Calculation
// Formula: 0.004 * FV + 0.986

// Calculate Land Surface Emissivity (EM) using the Fraction of Vegetation (FV).
// The 0.004 coefficient represents the emissivity variation due to vegetation,
// and the 0.986 represents the base emissivity for other surfaces.
var em = fv.multiply(ee.Number(0.004)).add(ee.Number(0.986)).rename('EM');

// Select Thermal Band (Band 10) and Rename It


var thermal = image.select('ST_B10').rename('thermal');
// Now, lets calculate the land surface temperature (LST)

// Formula: (TB / (1 + (λ * (TB / 1.438)) * ln(em))) - 273.15


var lst = thermal.expression(
'(TB / (1 + (0.00115 * (TB / 1.438)) * log(em))) - 273.15', {
'TB': thermal.select('thermal'), // Select the thermal band (TB)
'em': em // Assign emissivity (em)
}).rename('LST Bangalore 2024');

// Add the LST Layer to the Map with Custom Visualization Parameters
Map.addLayer(lst, {
min: 18.47, // Minimum LST value
max: 42.86, // Maximum LST value
palette: [
'040274', '040281', '0502a3', '0502b8', '0502ce', '0502e6',
'0602ff', '235cb1', '307ef3', '269db1', '30c8e2', '32d3ef',
'3be285', '3ff38f', '86e26f', '3ae237', 'b5e22e', 'd6e21f',
'fff705', 'ffd611', 'ffb613', 'ff8b13', 'ff6e08', 'ff500d',
'ff0000', 'de0101', 'c21301', 'a71001', '911003'
]}, 'Land Surface Temperature 2024');

You might also like