GEE(CCDC-3)——根据CCDC segment分割后的影像进行地类变化统计和绘制土地覆被变化地图

简介: GEE(CCDC-3)——根据CCDC segment分割后的影像进行地类变化统计和绘制土地覆被变化地图

简介

这一部分是最后一部分,因为这里我们已经根据前两个部分获取了分割后的影像,这里就可以选择指定的时间范围来查看其指定时间段的土地分类结果,同时可以导出研究区发生变化的时间结果。

之前两个部分请查看以下两篇博客

GEE(ccdc-1)——利用Landsat系列影像获取研究区范围的多波段影像(包含ccdc中已经定义的多波段)-CSDN博客 GEE(CCDC-2)——根据以获取的研究区CCDC系数进行土地覆被分类分析-CSDN博客

绘制土地覆被和土地覆被变化地图APP

https://code.earthengine.google.com/?scriptPath=users%2Fparevalo_bu%2Fgee-ccdc-tools%3AAPPS%2Flandcover_app

输入分割影像

这里是根据已经加载后的影像输入路径,然后点击load就可以进行后续参数的设定

分类制图

这里我们首先可以选择我们需要进行查看土地分类指定时间节点的影像,这个时间节点取决于你影像的时间范围,也就是CCDC-1教程中获取的时间影像的时间范围。如果你获取的仅仅是2000年的之后的影像,那么就无法获取1990年之间的影像。这里按照标准格式进行添加,可以选择添加多个时间段的影像,然后查看不同时间段的土地分类结果。

第二个选项我们可以查看在某个时间段内地类发生变化的影像,而分类选项我们根据我们分类对应的数字查看从某一类到某一类地类的转化过程,同时也支持多选,也就是你剋同时查看草地,林地转化成建筑物的结果,最后结果按照红色的栅格所展示。

不同时间段的土地分类结果

地类发生转变的部分

添加地类后的地图

所有代码

/**
 * @license
 * Copyright 2019 Boston University
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *    https://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *
 * @fileoverview User interface more making land cover maps with CCDC.
 * @author Eric L. Bullock, PhD
 */
//............. UTILITIES ................//
var utils = require('users/parevalo_bu/gee-ccdc-tools:ccdcUtilities/api')
var uiUtils = require('users/parevalo_bu/gee-ccdc-tools:ccdcUtilities/ui') 
 
//............. VARIABLES ................//
var properties = function() {
  this.coefs = [
    "INTP", "SLP", 
    "COS", "SIN", 
    "COS2", "SIN2",
    "COS3", "SIN3", 
    "RMSE"]
    
  this.segs = ["S1", "S2", "S3","S4", "S5", "S6"]
 
  /**
   * 可视化参数
   */ 
  this.visLabels = {
    fontWeight: 'bold', 
    fontSize: '14px', 
    padding: '4px 4px 4px 4px',
    border: '1px solid black',
    color: 'white',
    backgroundColor: 'black',
    textAlign: 'center',
    stretch: 'horizontal'
  }
  this.visLabelsSub = {
    fontWeight: 'bold', 
    fontSize: '14px', 
    padding: '4px 4px 4px 4px',
    color: 'black',
    backgroundColor: 'white',
    textAlign: 'center',
    stretch: 'horizontal'
  }
}
var horizontalStyle = {
  stretch: 'horizontal', 
  width: '100%'
}
var PROPS
var HELPER
var WIDGS = {}
var ccdParams = {}
var runParams = {}
var vizParams = {}
var globGeom = []
var INDICES = ['NDVI', 'NBR', 'EVI', 'EVI2', 'NDFI', 'GREENNESS', 'BRIGHTNESS', 'WETNESS']
var BANDS = ['BLUE','GREEN','RED', 'NIR', 'SWIR1', 'SWIR2'] 
var FULLBANDS = BANDS.concat(INDICES)
var BPBANDS = ['GREEN', 'RED', 'NIR', 'SWIR1',
相关文章
|
7月前
|
定位技术
ArcGIS地形起伏度+地形粗糙度+地表切割深度+高程变异系数提取
ArcGIS地形起伏度+地形粗糙度+地表切割深度+高程变异系数提取
2137 0
|
3月前
GEE图表——趋势线图表的加载和展示包含纵坐标间隔的设定(以某区域年均降水总量为例)
GEE图表——趋势线图表的加载和展示包含纵坐标间隔的设定(以某区域年均降水总量为例)
36 0
|
3月前
|
图形学 计算机视觉
GEE错误——如何将原有矢量将维度转化为地理坐标系,重投影坐标坐标无法实现?
GEE错误——如何将原有矢量将维度转化为地理坐标系,重投影坐标坐标无法实现?
34 0
|
4月前
|
数据处理
分面中嵌入趋势线
分面中嵌入趋势线
31 1
|
5月前
|
定位技术
基于ENVI实现栅格遥感影像按图层行列号与像元数量划定矩形研究区域并裁剪
基于ENVI实现栅格遥感影像按图层行列号与像元数量划定矩形研究区域并裁剪
|
5月前
ArcGIS矢量面要素中零碎小面积空洞区域补全与单独部分区域分离并剔除
ArcGIS矢量面要素中零碎小面积空洞区域补全与单独部分区域分离并剔除
|
10月前
|
编解码 定位技术
GIS空间分析 数字地形分析2 基本地形因子的提取
在本文中,你将学习到如何利用DEM数据在ArcGIS中提取出坡度、坡向等地形因子。
183 0
|
10月前
|
定位技术
ArcGIS:(模型构建器)如何获取各个区域下不同地物的面积?
ArcGIS:(模型构建器)如何获取各个区域下不同地物的面积?
146 0
|
7月前
|
存储 编解码 定位技术
案例!从天地图中提取全市的建筑物矢量轮廓-以苏州市为例
案例!从天地图中提取全市的建筑物矢量轮廓-以苏州市为例
116 1
|
10月前
|
数据挖掘
ArcGIS:如何进行栅格数据的拼接和裁剪、坡度坡向的提取、地形透视图的建立、等高线的提取、剖面图的创建?
ArcGIS:如何进行栅格数据的拼接和裁剪、坡度坡向的提取、地形透视图的建立、等高线的提取、剖面图的创建?
267 0