基于GEE平台采用Landsat8L2级的数据进行反照率的计算

案例：采用GEE平台基于Landsat8L2级的数据进行反照率的计算
代码如下:

// 定义点，坐标，坐标系
var p1 = ee.Geometry.Point([-88.241,41.8406], "EPSG:4326")
// 将点集成列表，可以写入多个点
var pts = ee.FeatureCollection(ee.List([ee.Feature(p1).set('name','p1')]))// 去云函数，以Landsat8为例
function maskL8sr(image){// Bits 4 and 3 are cloud shadow and cloud, respectively.var cloudShadowBitMask = (1 << 4);var cloudsBitMask = (1 << 3);// Get the pixel QA band.var qa = image.select('QA_PIXEL');// Both flags should be set to zero, indicating clear conditions.var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0).and(qa.bitwiseAnd(cloudsBitMask).eq(0));return image.updateMask(mask);
}// 定义NDVI计算函数，以Landsat8为例
function NDVI_LS(image){var ndvi = image.expression("0.13*B1+0.115*B2+0.143*B3+0.180*B4+0.281*B5+0.108*B6*0.042*B7",{"B1": image.select("SR_B1").multiply(0.0000275).add(-0.2),"B2": image.select("SR_B2").multiply(0.0000275).add(-0.2),
"B3": image.select("SR_B3").multiply(0.0000275).add(-0.2),
"B4": image.select("SR_B4").multiply(0.0000275).add(-0.2),
"B5": image.select("SR_B5").multiply(0.0000275).add(-0.2),
"B6": image.select("SR_B6").multiply(0.0000275).add(-0.2),
"B7": image.select("SR_B7").multiply(0.0000275).add(-0.2)}).rename('NDVI');return image.addBands(ndvi)
}// 产生数据集
var dataset = ee.ImageCollection("LANDSAT/LC08/C02/T1_L2").filterDate('2013-01-01','2013-12-31').filterBounds(pts)                .map(maskL8sr).map(NDVI_LS).select('NDVI')var ft = ee.FeatureCollection(ee.List([]))                
var fill = function(img, ini) {var inift = ee.FeatureCollection(ini)var ft2 = img.sampleRegions({collection:pts,properties:ee.List(['name']),scale:30});var date = img.date().format()var ft3 = ft2.map(function(f){return f.set("date", date)})return inift.merge(ft3)
}// 迭代获取影像集中每一幅影像的值
var newft = ee.FeatureCollection(dataset.iterate(fill, ft))// 结果输出
Export.table.toDrive({collection: newft,description: 'FSL_US_IB2',fileFormat: 'CSV'
});

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！