GEE（ccdc-1）——利用Landsat系列影像获取研究区范围的多波段影像（包含ccdc中已经定义的多波段）

//加载一个点用于定位
var geometry = /* color: #d63000 */ee.Geometry.Point([-69.68833790720782, -11.422826263974468]);
//加载两个定义好的函数库
var utils = require('users/parevalo_bu/gee-ccdc-tools:ccdcUtilities/inputs.js')
var ccdc_utils = require('users/parevalo_bu/gee-ccdc-tools:ccdcUtilities/ccdc.js')
var composite = require("users/google/toolkits:landcover/impl/composites.js").Composites
/**
 * UI app等定义可视化参数和样式。
 * 
 * 第一个字典包含标题参数，第二个字典包含 "运行 CCDC "按钮的样式。
 * 第二个字典包含 "运行 CCDC "按钮的样式。
 */ 
var visLabels = {
  fontWeight: 'bold', 
  fontSize: '14px', 
  width: '590px',
  padding: '4px 4px 4px 4px',
  border: '1px solid black',
  color: 'white',
  backgroundColor: 'black',
  textAlign: 'left'
  }
var runStyle = {
  fontWeight: 'bold', 
  fontSize: '14px', 
  width: '590px',
  padding: '4px 4px 4px 4px',
  border: '1px solid black',
  color: 'black',
  backgroundColor: 'white',
  textAlign: 'center'
  }  
  
var sensor = 'landsat'
/**
 * 全球 GLANCE 陆地网格。
 * 
 * 该特征集包含所有大陆与陆地重叠区域的
 * 与陆地重叠区域的所有大陆的网格系统。
 * 有关 GLANCE 和网格的更多信息，请点击此处： 
 *  http://sites.bu.edu/measures/
 */ 
var grids = ee.FeatureCollection('projects/GLANCE/GRIDS/GEOG_LAND/GLANCE_v01_GLOBAL_TILE_LAND')
/**
 * 用于放置 ui.Widgets 的面板。
 * 
 * 有一个放置所有内容的主面板，以及一个名为geoPanel 的面板，用于放置选择输出范围的选项。
 */ 
var mainPanel = ui.Panel({style: {width: '600px'}})
.add(ui.Label('Continuous Change Detection and Classification (CCDC)',visLabels))
.add(ui.Label('User interface for submitting CCDC'))
.add(ui.Label('Please see instructions'))
.add(ui.Label('Input Parameters',visLabels))
ui.root.add(mainPanel)
var geoPanel = ui.Panel()
/**
 * 定义全局变量 
 * 
 * 这样方便控件之间的交互变得更容易。
 */ 
//定义变量名称和常量变量
var outGeo, outGeos, outGeosSize, tool, drive, asset
var h = 0
var v = 0
var r = 0
var geoType = 0
var regionList = ['Select Method','Draw on Map','Single Tile','Tile Intersecting Point','Multiple Tiles','Draw Multiple Tiles on Map']
/**
 * 在地图上画出研究区的引用
 * 
 * Written by Justin Braaten (I think): https://emapr.github.io/LT-GEE/
 * https://code.earthengine.google.com/82b08b69bd596ada4747cb4bb7ea9526
 * 
 * 此功能可让您在地图上点击，点击 5 次后，这些位置就会连接起来，形成一个多边形。该多边形可用作 分析的地理范围。
 */ 
var DrawAreaTool = function() {
  var map = Map
  var drawingToolLayer = ui.Map.Layer({name: 'Area Selection Tool', visParams: {palette:'#4A8BF4', color:'#4A8BF4' }});
  this.map = map;
  this.selection = null;
  this.active = false;
  this.points = [];
  this.area = null;
  
  this.listeners = [];
  var tool = this;
  
  this.initialize = function() {
    this.map.onClick(this.onMouseClick);
    map.layers().reset()
    map.layers().set(0, drawingToolLayer);
  };
  
  this.startDrawing = function() {
    this.active = true;
    this.points = [];
    this.map.style().set('cursor', 'crosshair');
    drawingToolLayer.setShown(true);
  };
  
  this.stopDrawing = function() {
    tool.active = false;
    tool.map.style().set('cursor', 'hand');
    if(tool.points.length < 2) {
      return;
    }
    var closedPoints = tool.points.slice(0,-1);
    tool.area = ee.Geometry.Polygon(closedPoints)///.bounds();
    
    var empty = ee.Image().byte();
    var test = empty.paint({
      featureCollection: ee.FeatureCollection(tool.area),
      color: 1,
      width: 4
    });
  
    drawingToolLayer.setEeObject(test);
    tool.listeners.map(function(listener) {
      listener(tool.area);
    });
    runButton.widgets().get(0).setDisabled(false)
    runButton.widgets().get(0).style().set('backgroundColor','#5ab4ac')
  };
  
  this.onMouseClick = function(coords) {
    if(!tool.active) {
      return;
    }
    
    tool.points.push([coords.lon, coords.lat]);
    var geom = tool.points.length > 1 ? ee.Geometry.LineString(tool.points) : ee.Geometry.Point(tool.points[0]);
    drawingToolLayer.setEeObject(geom);
    if(tool.points.length > 4) {
      tool.stopDrawing();
    }
  };
  
  this.onFinished = function(listener) {
    tool.listeners.push(listener);
  };
  
  this.initialize();
};
/**
 * 绘制输出范围的辅助函数。
 * 
 * 使用DrawAreaTool() 函数绘制输出范围。
 */ 
function doAreaTool() {
  geoPanel.add(ui.Label('Slowly click five points on the map and the application will generate a rectangle for the output extent geometry.'))
  tool = new DrawAreaTool();
  tool.startDrawing();
  tool.onFinished(function(geometry) {
    outGeo = ee.Feature(geometry);
  });
  
}
/**
 * 绘制输出范围和查找重叠瓦片的辅助函数
 * 
 * 使用DrawAreaTool() 函数绘制输出范围。
 */ 
function drawMultipleTiles() {
  geoPanel.add(ui.Label('Slowly click five points on the map and the application will generate a rectangle for the output extent geometry.'))
  tool = new DrawAreaTool();
  tool.startDrawing();
  tool.onFinished(function(geometry) {
    var tempGeo = grids.filterBounds(geometry);
    tempGeo.size().evaluate(function(val) {
      if (val > 0) {
        outGeos = ee.FeatureCollection(tempGeo)
        Map.addLayer(outGeos, {},'Output Geometry')
        Map.centerObject(outGeos)
        runButton.widgets().get(0).setDisabled(false)
        runButton.widgets().get(0).style().set('backgroundColor','#5ab4ac')
        outGeosSize