L4578ndvi分析+m-k分析

var geometry = 
    /* color: #d63000 */
    /* shown: false */
    /* displayProperties: [
      {
        "type": "rectangle"
      }
    ] */
    ee.Geometry.Polygon(
        [[[-67.62548053231656, -25.849835182304385],
          [-67.62548053231656, -26.114016902313434],
          [-67.3211433157106, -26.114016902313434],
          [-67.3211433157106, -25.849835182304385]]], null, false);
//Center Map 
Map.centerObject(geometry, 10);
 
 
 
var months=ee.Filter.calendarRange(9, 12, "month")
 
//Cloud filter based on Cloud_cover of the whole scene
var cloudcover= 80;
 
 
//
///TEMPORAL SERIES OF LANDSAT/// 
//
 
 
//-----------------------------------------------------------------------------------------------------------------------------
/// Landsat - 4  //
//-----------------------------------------------------------------------------------------------------------------------------
 
var Landsat4 = ee.ImageCollection('LANDSAT/LT04/C02/T1_L2')  
               .filterBounds(geometry) 
               .filterDate("1982-08-22", "1993-06-24") //periodo de L5  "2012-06-01"
               .filterMetadata("CLOUD_COVER", 'less_than',cloudcover)
               .filter(months);
 
// Scaling factors. according to https://developers.google.com/earth-engine/datasets/catalog/LANDSAT_LT05_C02_T1_L2
function applyScaleFactorsL5(image) {
  var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
  var thermalBand = image.select('ST_B6').multiply(0.00341802).add(149.0);
  return image.addBands(opticalBands, null, true)
              .addBands(thermalBand, null, true);
}
 
 Masks for NDVI 
function maskClouds(image) {             
  var cloudShadowBitMask1 = (1 << 4);    
  var cloudsBitMask1 = (1 << 3);     
  var highcloud1 = (1 << 9);
  var snow = (1 << 5);
  var qa1 = image.select('QA_PIXEL').rename('Cloud');    
  var mask1 = qa1.bitwiseAnd(cloudsBitMask1).eq(0)               
               .and(qa1.bitwiseAnd(highcloud1).eq(0)            
               .or(qa1.bitwiseAnd(cloudShadowBitMask1).eq(0)))
               .and(qa1.bitwiseAnd(snow).eq(0));
  var maskedComposite = image.updateMask(mask1);
return maskedComposite.addBands(mask1);            
}
 
 
 
// Apply scale factors
Landsat4 = Landsat4.map(applyScaleFactorsL5);
 
// Apply the mask clouds
var collectionL4 = Landsat4.map(maskClouds);
 
// I select bands for the calculations of the indices
var col_L4 = collectionL4.select(['SR_B1','SR_B2','SR_B3','SR_B4','SR_B5', 'SR_B7', 'ST_B6']);
 
 
//--------------------------------------------------------------------------------------