PIE-Engine利用modis计算ndvi

借助国产遥感计算云服务平台PIE-Engine，利用MOD09A1这套数据获取长时序年平均NDVI，因为PIE才刚刚起步，所以找了好久都没有特别适合的代码，以下是我摸索的成果，希望帮助到有需要的人，如果发现问题，也欢迎指证。

var map1 = pie.FeatureCollection('user/xxx/xxx_SHP');
var roi = map1.first() //转化为对象.geometry(); //转化为几何形状// 创建一个函数用来去除有云的像元
// 参考文章链接：https://blog.csdn.net/little00bee/article/details/105801593
// 参考GEE去云函数：https://blog.csdn.net/qq_41664159/article/details/124130777?spm=1001.2101.3001.6661.1&utm_medium=distribute.pc_relevant_t0.none-task-blog-2%7Edefault%7ECTRLIST%7ERate-1-124130777-blog-128925130.235%5Ev32%5Epc_relevant_default_base3&depth_1-utm_source=distribute.pc_relevant_t0.none-task-blog-2%7Edefault%7ECTRLIST%7ERate-1-124130777-blog-128925130.235%5Ev32%5Epc_relevant_default_base3
// GEE中如何使用 QA 波段和位掩码：https://blog.csdn.net/qq_31988139/article/details/120315322?spm=1001.2101.3001.6661.1&utm_medium=distribute.pc_relevant_t0.none-task-blog-2%7Edefault%7ECTRLIST%7EPayColumn-1-120315322-blog-124130777.235%5Ev32%5Epc_relevant_default_base3&depth_1-utm_source=distribute.pc_relevant_t0.none-task-blog-2%7Edefault%7ECTRLIST%7EPayColumn-1-120315322-blog-124130777.235%5Ev32%5Epc_relevant_default_base3&utm_relevant_index=1
var maskClouds = function(image) {// 选择质量评估波段var QA = image.select('sur_refl_state_500m')// 1<<10表示二进制第10位，1表示第10位表示有云，0表示第10位无云。var bitMask = 1 << 10;// 使得检测出含云像元置为0，进行掩膜去除含云var cloudMask = QA.bitwiseAnd(bitMask).eq(0);return image.updateMask(cloudMask);
}//加载Terra星全球500m地表反射率8天合成产品(MOD09A1 V6)
var imgs = pie.ImageCollection('USGS/MOD09A1/006').filterDate('2000-01-01', '2000-12-31').select(['sur_refl_b01', 'sur_refl_b02','sur_refl_state_500m']).filterBounds(roi).map(maskClouds).mean()// .mosaic() //掩膜;       function NDVI_new(image){var b1 = image.select("sur_refl_b01").toFloat();var b2 = image.select("sur_refl_b02").toFloat();var ndvi = (b2.subtract(b1)).divide(b2.add(b1));return ndvi;
}var ndvi = NDVI_new(imgs);print(ndvi);                 //地图显示中心
Map.centerObject(roi, 5);//NDVI绘制样式
var visParamNDVI = {min: -0.2,max: 0.8,palette: 'CA7A41, CE7E45, DF923D, F1B555, FCD163, 99B718, '+'74A901, 66A000, 529400,3E8601, 207401, 056201, 004C00,'+'023B01, 012E01, 011D01, 011301'
};
//加载影像集合的最大值
// Map.addLayer(img.mean().clip(roi), visParamNDVI, "NDVI");
Map.addLayer(ndvi.clip(roi), visParamNDVI, "NDVI");//加载显示中国边界图
Map.addLayer(roi, {color:"ffff00ff", fillColor:"00000000"}, "研究区");/*** 导出影像，参数分别为：*  image：导出影像*  description：导出任务描述也就是任务名称*  assetId：导出的《资源》中数据集合路径*  region：导出范围*  scale：分辨率30米* */
Export.image({image:ndvi.clip(roi),description:'ndvi2000',region:roi,crs:'EPSG:4326',scale:500,
});

PIE-Engine利用modis计算ndvi相关推荐

Google Earth Engine（GEE）计算NDVI、EVI、NDBI等指数并相关性分析
Google Earth Engine(GEE)计算NDVI.EVI.NDBI等指数并相关性分析 NDVI为归一化植被指数.EVI为增强植被指数,据GEE官方介绍NDTI为归一化热差指数,用于区分城市 ...
Google Earth Engine（GEE）——计算NDVI\EVI\RVI\DVI\SAVI计算并下载（2）
有粉丝还问到关于Landsat8如何计算,上一次的文章提到了关于各类植被指数的下载和加载:Google Earth Engine(GEE)--计算NDVI\EVI\RVI\DVI\SAVI归一化植被指 ...
利用IDL计算植被覆盖度(VFC)
0. 前言正巧IDL实验课考核的作业是利用4个Function和主Pro过程写一个遥感图像处理的代码,要求是前一个方法的输出是另一个方法的输入.以前一直想着能不能计算NDVI和植被覆盖度(VFC ...
ArcGIS计算NDVI为什么只有1和-1及0值
在遥感影像的应用中,归一化植被指数是非常常用的一个指标,是反映农作物长势和营养信息的重要参数之一,在专业处理遥感影像的应用软件(ENVI.ERDAS等)中都有直接的工具可以来计算NDVI.NDVI的计 ...
Python批量计算NDVI
PYTHON批量计算NDVI 做了少量修改,剔除了异常值,执行代价时需要更换影像对应波段及文件存储位置 import os import numpy as np from osgeo import g ...
python计算csv文件内的数据_Python利用pandas计算多个CSV文件数据值的实例
功能:扫描当前目录下所有CSV文件并对其中文件进行统计,输出统计值到CSV文件 pip install pandas import pandas as pd import glob,os,sys in ...
微带线特性阻抗计算公式_利用HFSS计算微带线的特性阻抗
系别电子工程系课程名称 <微波技术与天线>实验班级实验名称利用 HFSS 计算微带线的特性阻抗姓名实验时间学号指导教师报告内容一.实验目的和任务 1 .了解微带 ...
python ks值计算_利用Python计算KS的实例详解
在金融领域中,我们的y值和预测得到的违约概率刚好是两个分布未知的两个分布.好的信用风控模型一般从准确性.稳定性和可解释性来评估模型.sOf免费资源网一般来说.好人样本的分布同坏人样本的分布应该是有很 ...
10分钟上线 - 利用函数计算构建微信小程序的Server端
摘要: 阿里云函数计算是一个事件驱动的全托管计算服务.通过函数计算,您无需管理服务器等基础设施,只需编写代码并上传.微信小程序是一种不需要下载安装即可使用的应用,它可以在微信内被便捷地获取和传播. 当 ...

PIE-Engine利用modis计算ndvi

PIE-Engine利用modis计算ndvi相关推荐

最新文章

热门文章