SCRIPT URBAN HEAT ISLAND KABUPATEN KUNINGAN

 SCRIPT URBAN HEAT ISLAND KABUPATEN KUNINGAN

var maskL8 = function(image) {

  var qa = image.select('BQA');

  var mask = qa.bitwiseAnd(1 << 4).eq(0);

  return image.updateMask(mask);

}

var dataset = ee.ImageCollection("LANDSAT/LC08/C01/T1")

 .filterDate('2019-01-01', '2019-12-31')

 .filterBounds(geometry)

  .map(maskL8)

 .select ("B10");

 

var first = dataset.first();

var K1 = first.get('K1_CONSTANT_BAND_10');

var K2 = first.get ('K2_CONSTANT_BAND_10');

var A = first.get ('RADIANCE_ADD_BAND_10');

var M = first.get ('RADIANCE_MULT_BAND_10');

print (A);

print (M);

print (K1);

print (K2);

var Celcius = dataset.map(function(img){

  var id= img.id();

  return img.expression ('((1321.08/(log(774.89/((TIR*0.0003342)+0.1)+1)))-272.15)'

,{'TIR':img})

  .rename('B10')

  .copyProperties (img, ['system:time_start'])

});

var celciusmean = Celcius.mean().clip(geometry);

print (celciusmean)

var param = {min:10, max:33, palette:['blue','limegreen','yellow','darkorange', 'red']};

Map.addLayer(celciusmean, param, 'SUHU RATA-RATA KUNINGAN 2019 Januari-Desember'

);

// Load a landsat image and select three bands.

var landsat = ee.Image('LANDSAT/LC08/C01/T1_TOA/LC08_123032_20140515')

  .select(['B4', 'B3', 'B2']);

// Create a geometry representing an export region.

//var geometry = ee.Geometry.Rectangle([116.2621, 39.8412, 116.4849, 40.01236]);

// Export the image, specifying scale and region.

Export.image.toDrive({

  image: celciusmean,

  description: 'KNG19',

  scale: 60,

  region: geometry,

});