【GEE笔记2】数据整合Reducer/遍历循环.map

原创已于 2022-08-06 21:15:35 修改 · 2.6k 阅读

26 ·

CC 4.0 BY-SA版权

文章标签：

#GEE #遥感 #JS

于 2022-07-31 22:26:15 首次发布

遥感同时被 2 个专栏收录

38 篇文章

订阅专栏

GEE

15 篇文章

订阅专栏

该博客介绍了如何使用Google Earth Engine（GEE）进行遥感数据分析，具体操作包括：通过imageCollection筛选特定区域和时间范围的卫星影像，计算NDVI并将其作为新波段添加到影像中，然后利用reduceRegion方法计算每个影像在ROI上的NDVI平均值，并将结果作为属性字段。整个过程展示了GEE中的map函数在并行处理大量影像数据时的高效性。

Reducer是GEE中按照时间、空间、波段等进行整合数据的方法；类似于从部分到整体。

ee.Reducer()可以指定数据整合的方式，如：max\min\mean\etal。

时间：imageCollection.reduce()
空间：image.reduceRegion() ; image.reduceNeighborhood()
波段：image.reduce()
Feature Collection 的属性：Feature Collection.reduceColumns()

以计算NDVI数值并添加至影像合集为例

具体需求：
为一个影像集合添加NDVI波段，同时计算每一景影像上roi的NDVI均值

var mean = ndvi.reduceRegion({  
   reducer: ee.Reducer.mean(),  
   geometry: roi,   
   scale: 30  
});  
print("reduceRegion value is: ", mean);

map遍历循环

map：用的最多，执行效率最快，相当于“并行执行任务”，能用 map,就用 map，但是函数里一般不能用 print,Map,toDrive 等方法
for、while：就是基本的循环结构语句，但是在GEE里很少用，基本用 map 替代，因为是顺序执行代码，效率太慢
iterate：“迭代”，顺序执行，每次执行都在前面执行的结果上继续执行，很遗憾，这次的例子没有体现“迭代”的意思
evaluate：主要用来异步导出和显示，加快速度，一般没有返回值，结合 for 和 map 一起用，这里的 map 里可以有 print,Map,toDrive 方法

完整代码：
roi是我自己画的
在这里插入图片描述

var roi = ee.FeatureCollection(geometry);  
Map.addLayer(roi, {color: "red"}, "roi"); // 显示roi 
Map.centerObject(roi, 9); // 以roi为中心，缩放等级9
// 筛选影像 
var Col = ee.ImageCollection("LANDSAT/LC09/C02/T1_TOA")  
              .filterBounds(roi)  
              .filterDate("2022-1-1", "2022-7-1")  
              // 循环遍历影像集合，为每景影像计算并添加一个NDVI波段
              .map(function(image){  
                var ndvi = image.normalizedDifference(["B5", "B4"])  
                                .rename("NDVI");  
                return image.addBands(ndvi);  
              })  
              // 选择NDVI波段用于后续计算
              .select("NDVI")  
              // 循环遍历影像集合，计算每景影像roi的NDVI均值，并赋为新的属性字段“ndvi”
              .map(function(image) {  
                var dict = image.reduceRegion({  
                 reducer: ee.Reducer.mean(),  
                 geometry: roi,  
                 scale: 30  
                });
                var ndvi = ee.Number(dict.get("NDVI"));  
                image = image.set("ndvi", ndvi);  
                return image;  
              });  
// 可视化参数，加载显示新的影像集合，使用first()，返回集合的第一个非空元素，即第一景影像
print("Col", Col);  
var visParam = {  
 min: -0.2,   
 max: 0.8,  
 palette: ["FFFFFF", "CE7E45", "DF923D", "F1B555", "FCD163",   
            "99B718", "74A901", "66A000", "529400", "3E8601",   
            "207401", "056201", "004C00", "023B01", "012E01",   
            "011D01", "011301"]  
};  
Map.addLayer(Col.first().clip(roi), visParam, "Col");