《遥感云计算与科学分析——应用与实践》中常用算法与实例教程代码整理

一、常用算法

1.去云

Landsat TOA

https://code.earthengine.google.com/feb1d236038ef0f50ecd710ce5e9556a

Landsat SR

https://code.earthengine.google.com/5841963e413265c72741525f1b6607bd

Sentinel-2

https://code.earthengine.google.com/6c3d05daed20dd5dc6e60e911771ff0c

2.时间序列数据合成

合成年度数据

https://code.earthengine.google.com/4be04f02b8326cb5804a82bb7b3bb08c

合成月度数据——会出现空波段影像

https://code.earthengine.google.com/cf237797cf5f19060ba4413b6f82cc01

Filter时间序列合成(月) ——不出现空波段影像

https://code.earthengine.google.com/cf70daeb2f03e063b651ee389d1c2008

Join时间序列合成(月)——不出现空波段影像

https://code.earthengine.google.com/a3b1ba75364a079d7359437c0a52b134

3.形态学变化操作

开操作：先腐蚀后膨胀（消除细小物体）

闭操作：先膨胀后腐蚀（填充物体内细小空洞）

腐蚀focal_min，膨胀focal_max

https://code.earthengine.google.com/0fb2f77ef6ad0c00d3f189f17341b7b3

4.统计分析操作

统计某地影像覆盖度

.reduce(ee.Reducer.count())计算每个像素有效的像素数

https://code.earthengine.google.com/c6d7e8bdf4b177a7bc77fda4fc3c95bf

统计某个区域内最大值或最小值的经纬度

ee.Reducer.max()或ee.Reducer.min()计算极值

ee.Image.pixelLonLat()获取经纬度

https://code.earthengine.google.com/7c87e592ff46aa75bdb028392baf8803

5.缨帽变换

https://code.earthengine.google.com/b8a328b3e7d16ed2c737c8ccd114cec4

6.影像集合数据导出

https://code.earthengine.google.com/2603384e06c7e0fdc27b3343f450d475

二、实例教程

1.动态变化监测

黄河入海口30年变化

https://code.earthengine.google.com/51fbcea132cd65c7b4068696b40d6215

2.遥感地物分类

ee.Clusterer：非监督分类，包括K-mean等

ee.Classifier：监督分类，包括svm、cart、decisionTree等

ee.ConfusionMatrix：混淆矩阵，做精度评价，提取精度、Kappa系数

2.1非监督分类

非监督分类

https://code.earthengine.google.com/df35058a7620f442b177ec7a65eb8d79

添加图例后的非监督分类

https://code.earthengine.google.com/b201d170e0a76ca2579787a09323cf7f

2.2监督分类

利用已有的分类结果随机采样（MODIS）

https://code.earthengine.google.com/87789ca41109750ec94b04ff068e1209

监督分类

https://code.earthengine.google.com/72fcf9641dcd9e859b0177ed314be828

2.3面向对象分类

先分割再分类ee.Algorithms.Images.Segmention

https://code.earthengine.google.com/fb946692b5ed34830d42977a4fd31e94