PIE-Engine 哨兵2FAI指数

最新推荐文章于 2024-07-23 13:00:00 发布

clearlove2015

最新推荐文章于 2024-07-23 13:00:00 发布

阅读量197

点赞数

文章标签：经验分享

本文链接：https://blog.csdn.net/clearlove2015/article/details/131674209

版权

start = "2018-04-03"
end = "2018-07-03"
var roi = pie.Geometry.Polygon([
    [
        [
            119.76725495049737,
            31.59154939469161
        ],
        [
            120.68907544090496,
            31.59154939469161
        ],
        [
            120.68907544090496,
            30.88772538310468
        ],
        [
            119.76725495049737,
            30.88772538310468
        ],
        [
            119.76725495049737,
            31.59154939469161
        ]
    ]
], null);
var faiprm = {
    palette: [
'#edf8e9','#c7e9c0','#a1d99b','#74c476','#41ab5d','#238b45','#005a32'
    ]
}
//哨兵数据
var image = pie.ImageCollection("S2/L2A")
                                .filterBounds(roi)
                                .filterDate(start,end)//作物生长
                                .filter(pie.Filter.lt("cloudyPixelPercentage", 10))//云盖量属性小于5
                                .select(["B1","B2","B3","B4","B5","B6","B7","B8","B9","B10","B11","B12"])//选择需要的波段
                                .mosaic(["B1","B2","B3","B4","B5","B6","B7","B8","B9","B10","B11","B12"])//选择需要的image波段并拼接影像
                                .clip(roi);//根据研究区范围裁剪影像
print(image);
Map.addLayer(image.select(["B5","B4","B3"]),{min:0,max:3000},"影像",false);
// 提取FAI
function fai(img){
    var nir1 = img.select(["B4","B8","B11"]).expression(
        'B4+((B11-B4)*(B8-B4)/(B11-B4))',{
            B4: img.select("B4"),
            B11:img.select("B11"),
            B8: img.select("B8"),
        }).rename('nir1')
    var nir = img.select('B11')
    var fai = nir.subtract(nir1).rename('fai')
    return fai
}
var fai1 = fai(image)
print(fai1)

var water = pie.Algorithm.Sentinel2.MNDWI(image).gt(0.3)
var ndvi = pie.Algorithm.Sentinel2.NDVI(image).updateMask(water)
var fai2 = fai1.updateMask(water)
print(fai2)

Map.addLayer(fai2,faiprm,"蓝藻")

Export.imageToCloud({
  image:fai2,
  description: '',
  region:roi,
  scale:10,
})