
//加载边界
var bj = pie.FeatureCollection("NGCC/CHINA_COUNTY_BOUNDARY")
.filter(pie.Filter.eq("name","雄县")) //筛选 eq是等于的意思
.first().geometry();//first数据集里面的第一个选出来 ;geometry 特征集类型转换
Map.centerObject(bj,8); //定位
Map.addLayer(bj,{fillColor:"00000000"},"边界") //边界图层,填充颜色,在地图上显示的名称
//ndvi的计算
function NDVI(img){
//var B5_SR=img.select(["B5"]).multiply(0.0000275).add(-0.2);
//var B4_SR=img.select(["B4"]).multiply(0.0000275).add(-0.2);
var ndvi=img.normalizedDifference(["B5","B4"]).rename ("ndvi");//第五第四进行计算
return img.addBands(ndvi);
}
for(var i=2000,j=2023;i<=2023;i++){//i开始年份 j结束年份
var startdate=pie.Date.fromYMD(i,01,01)
var enddate=pie.Date.fromYMD(i,12,31)
//把一年的所有影像检索出来,计算NDVI,采用最大合成,获取区域里一年的植被生长情况
var dataset; // 定义dataset变量
if (i < 2013) {
dataset = pie.ImageCollection("LT05/02/SR")
.filterBounds(bj) // 按边界检索
.filterDate(startdate, enddate) // 按日期检索
.map(NDVI) // 计算NDVI
.select(["ndvi"]) // 选择NDVI波段
.max(); // 最大值合成
} else if (i >= 2013) {
dataset = pie.ImageCollection("LC08/02/SR")
.filterBounds(bj) // 按边界检索
.filterDate(startdate, enddate) // 按日期检索
.map(NDVI) // 计算NDVI
.select(["ndvi"]) // 选择NDVI波段
.max(); // 最大值合成
}
if(i==2000){
var ndvi_time=dataset;
print(ndvi_time);
}
else{
var ndvi_time=ndvi_time.addBands(dataset);
print(ndvi_time);
}
}
//数据输出
Export.image({
image: ndvi_time,
description: "ndvi_time_year_new",
//assetId: "24287496921681920123",
region: bj, // 导出矢量范围,这里使用features的几何体
scale: 30
});
//var images = pie.ImageCollection("USGS/MOD13Q1/061") //陆地