Annexes
rho = b4.multiply(math.pi)
TPM = landsat_season.expression('A*(rho/(1 - (rho/C)))', {
'A' : A,
'rho' : rho,
'C' : C
}).clip(roi);
task2 = ee.batch.Export.image.toDrive(image=TPM,
folder=folder,
region=roi,
scale=30,
fileNamePrefix='TPM_%d' %(m))
task2.start()
#--------- Sea surface Temperature algorithme ------b10 = landsat_season.select('B10').multiply(0.1).subtract(273.5).c
lip(roi)
b11 = landsat_season.select('B11').multiply(0.1).subtract(273.5).c
lip(roi)
SST = landsat_season.expression('b10 + (2.946*(b10 - b11))-0.038',
{
'b10' : b10,
'b11' : b11
}).clip(roi);
task3 = ee.batch.Export.image.toDrive(image=SST,
folder=folder,
region=roi,
scale=30,
fileNamePrefix='SST_%d' %(m))
task3.start()
Initialisation de l'extraction
oceanL8(landsat, coast_geo, 'Landsat_images')
Aperçu des résultats
algiers_map.addLayer(chl, {min:0,max:2,'palette':['green','yellow','orange
','red','blue']})
algiers_map