Hello Mr.Taher Mostafapour

This is Durga Prasad Jakka from India

This is the code for lst calculation for landsat images using Google Earth engine

Just remove visualisation parameters from Map.addLayer() function which produce error for you and you add your own visualisation parameters

//Importing Landsat Images

var Images = ee.ImageCollection("LANDSAT/LT05/C02/T1_L2")

.filterBounds(ROI)

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

.filterMetadata('CLOUD_COVER', 'less_than', 5)

.mean();

Map.addLayer(Images.clip(ROI),imageVisParam,'Landsat5');

print(Images);

//median (ndvi)

{

var ndvi = Images.normalizedDifference(['SR_B4', 'SR_B3']).rename('NDVI');

var ndvi_clip = ndvi.clip(ROI)

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

print(ndvi,'ndvi');

Map.addLayer(ndvi_clip, ndviParams, 'ndvi');

}

//select thermal band 10(with brightness tempereature), no calculation

var thermal= Images.select('ST_B6').multiply(0.00341802).add(149.0);

var thermal_clip = thermal.clip(ROI)

var b10Params = {min: 291.918, max: 302.382, palette: ['blue', 'white', 'green']};

Map.addLayer(thermal_clip, b10Params , 'thermal');

// find the min and max of NDVI

{

var min = ee.Number(ndvi_clip.reduceRegion({

reducer: ee.Reducer.min(),

geometry: ROI,

scale: 30,

maxPixels: 1e9

}).values().get(0));

print(min, 'min');

var max = ee.Number(ndvi_clip.reduceRegion({

reducer: ee.Reducer.max(),

geometry: ROI,

scale: 30,

maxPixels: 1e9

}).values().get(0));

print(max, 'max')

}

//fractional vegetation

{

var fv =(ndvi_clip.subtract(min).divide(max.subtract(min))).pow(ee.Number(2)).rename('FV');

print(fv, 'fv');

Map.addLayer(fv);

}

//Emissivity

var a= ee.Number(0.004);

var b= ee.Number(0.986);

var EM=fv.multiply(a).add(b).rename('EMM');

var imageVisParam3 = {min: 0.9865619146722164, max:0.989699971371314};

Map.addLayer(EM, imageVisParam3,'EMM');

//LST in Celsius Degree bring -273.15

//NB: In Kelvin don't bring -273.15

var LST = thermal.expression(

'(Tb/(1 + (0.00115* (Tb / 1.438))*log(Ep)))-273.15', {

'Tb': thermal.select('ST_B6'),

'Ep': EM.select('EMM')

}).rename('LST');

//min max LST

{

var min = ee.Number(LST.reduceRegion({

reducer: ee.Reducer.min(),

geometry: ROI,

scale: 30,

maxPixels: 1e9

}).values().get(0));

print(min, 'minLST');

var max = ee.Number(LST.reduceRegion({

reducer: ee.Reducer.max(),

geometry: ROI,

scale: 30,

maxPixels: 1e9

}).values().get(0));

print(max, 'maxLST')

}

Map.addLayer(LST, {min: 12.19726940544291, max:41.6701111498399, 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'

]},'LST');

Export.image.toDrive({

image: LST,

description: 'LST',

folder: "image EE",

scale: 30,

region: ROI,

fileFormat: 'GeoTIFF',

formatOptions: {

cloudOptimized: true

}

});

Which software you want to use? Please be more specific in your question.

Arc Gs pro or Envi or Matlab.