首页 > 其他分享 >GEE 案例:利用2001-2024年的MODIS数据长时序ndvi指数归一化后的结果分析

GEE 案例:利用2001-2024年的MODIS数据长时序ndvi指数归一化后的结果分析

时间:2024-09-18 12:24:16浏览次数:10  
标签:MODIS gte vci median 2024 归一化 var ndvi where

目录

简介

指数

数据

代码

结果


简介

利用2001-2024年的MODIS数据长时序ndvi指数归一化后的结果分析,并加载时序图。

指数

NDVI指数(Normalized Difference Vegetation Index)是用来评估地表植被覆盖度和健康程度的指数。它通过计算红光和近红外光反射率的差异来衡量植被的光合作用活性和生长状况。NDVI值范围从-1到1,较高的值表示较高的植被覆盖度和健康程度。

VCI指数(Vegetation Condition Index)是衡量植被状况的指数。它是基于NDVI指数计算得出的,用来表征植被的生长和状况。VCI指数通常表达为百分比形式,数值越高表示植被状况越好,数值越低表示植被状况越差。

VCI指数是根据NDVI的变化来计算的,通常将一段时间(如一年)内的NDVI数据与长期平均值进行比较,得出植被状况指数。在农业和生态环境监测中,VCI指数被用来评估植被的退化状况、旱情监测以及预测植被的生长情况。

总的来说,NDVI指数用来衡量植被覆盖程度和健康程度,而VCI指数则是根据NDVI的变化计算得出的,用来评估植被的生长和状况。

数据

MODIS/061/MOD13A2数据是来自于MODIS(Moderate Resolution Imaging Spectroradiometer)遥感传感器的产品之一。MODIS是美国国家航空航天局(NASA)和美国国家环境信息中心(NEIC)联合研发的一款具有高空间分辨率和全球覆盖能力的遥感传感器。

MOD13A2数据是基于MODIS传感器观测到的地表植被指数(Vegetation Index,VI)数据集。它提供了有关地表植被覆盖度和健康状况的信息。MOD13A2数据的空间分辨率为500米,时间分辨率为16天。数据包括以下几个主要的植被指数:

  1. 植被指数(Vegetation Index,VI):包括了常用的NDVI(Normalized Difference Vegetation Index)、EVI(Enhanced Vegetation Index)等。这些指数用于评估植被的光合活性和健康状况,反映了地表植被的覆盖度和生长状态。

  2. 植被覆盖度(Fractional Vegetation Cover,FVC):表示地表被植被覆盖的比例,可用于监测植被覆盖变化和评估土地利用。

  3. 植被生长指标(Vegetation Growth Index,VGI):用于评估地表植被的生长情况,根据植被指数的变化来反映植被的生长速率和季节变化。

MOD13A2数据的应用广泛,包括农业监测、水资源管理、气候变化研究、生态环境监测等。通过分析MOD13A2数据,可以获得关于地表植被的信息,为决策制定和资源管理提供科学依据。

代码


var feature = ee.FeatureCollection("FAO/GAUL/2015/level0")

var point = /* color: #98ff00 */ee.Geometry.Point([-89.8212566433932, 40.15005408150289]);

var roi = feature.filterBounds(point);

Map.centerObject(roi)
Map.addLayer(roi)

// define start and end time

var time_start = '2001', time_end = '2003'

// call modis ndvi images 

var ndvi = ee.ImageCollection("MODIS/061/MOD13A2")
.select(['NDVI'])
.filterDate(time_start, time_end);

// clculate the vci
// vegetation condition index
// vci = ((ndvi - min)/(max - min)) * 100

var ndvi_min = ndvi.min().multiply(0.0001);
var ndvi_max = ndvi.max().multiply(0.0001);

var vci = ndvi.map(function(img){
  var band = img.multiply(0.0001);
  var index = band.expression('((ndvi - min)/(max - min))*100.0',{
    'ndvi': band, 'min': ndvi_min, 'max': ndvi_max
    }).rename('vci');
  return index
  .copyProperties(img, ['system:time_start','system:time_end'])
  });
  
// drought mapping using vci

var vci_median = vci.median();

Map.addLayer(vci_median.clip(roi),[],'vci_median',false)

// histogram

print(
  ui.Chart.image.histogram(vci_median, roi, 1000)
  )

// vci classification

// making a constant image 

var cons = ee.Image.constant(0);

// define classes

var extreme = cons.where(vci_median.gte(0).and(vci_median.lt(10)), 1);
var severe = extreme.where(vci_median.gte(10).and(vci_median.lt(20)), 2);
var moderate = severe.where(vci_median.gte(20).and(vci_median.lt(30)), 3);
var light = moderate.where(vci_median.gte(30).and(vci_median.lt(40)), 4);
var no1 = light.where(vci_median.gte(40).and(vci_median.lt(60)), 5);
var no2 = no1.where(vci_median.gte(60).and(vci_median.lt(80)), 6);
var no3 = no2.where(vci_median.gte(80), 7);

Map.addLayer(moderate.clip(roi),{min:1, max: 7},'drought_map',false)


// vci classification for each date


var time_start = '2001', time_end = '2024'

var ndvi2 = ee.ImageCollection("MODIS/061/MOD13A2")
.select(['NDVI'])
.filterDate(time_start, time_end);

var ndvi_min2 = ndvi2.min().multiply(0.0001);
var ndvi_max2 = ndvi2.max().multiply(0.0001);

var vci2 = ndvi2.map(function(img){
  var band = img.multiply(0.0001);
  var index = band.expression('((ndvi - min)/(max - min))*100.0',{
    'ndvi': band, 'min': ndvi_min2, 'max': ndvi_max2
    }).rename('vci');
  return index
  .copyProperties(img, ['system:time_start','system:time_end'])
  });
  
var vci_class = vci2.map(function(img){
  var cons = ee.Image.constant(0);
  var extreme = cons.where(vci_median.gte(0).and(vci_median.lt(10)), 1);
  var severe = extreme.where(vci_median.gte(10).and(vci_median.lt(20)), 2);
  var moderate = severe.where(vci_median.gte(20).and(vci_median.lt(30)), 3);
  var light = moderate.where(vci_median.gte(30).and(vci_median.lt(40)), 4);
  var no1 = light.where(vci_median.gte(40).and(vci_median.lt(60)), 5);
  var no2 = no1.where(vci_median.gte(60).and(vci_median.lt(80)), 6);
  var no3 = no2.where(vci_median.gte(80), 7);
  return no3
  .copyProperties(img, img.propertyNames())
  });
  
var vci_map = vci_class.mode();

Map.addLayer(vci_map.clip(roi),{palette: ['black','brown','red','orange','yellow', 'lightgreen','darkgreen'],
        min: 1, max: 7},'vci_mode', false)

Export.image.toDrive({
  image: vci_map.clip(roi),
  description: 'vci_map',
  region: roi,
  maxPixels: 1e13,
  crs: 'EPSG:4326',
  folder: 'drought',
  scale: 1000
  })


// drought area calculation

var drought_area = (ee.Image.pixelArea().divide(1e6)).addBands(vci_map);

print(
  ui.Chart.image.byClass(drought_area, 'constant', 
  roi, ee.Reducer.sum(), 1000, ['nan','extreme drought','severe drought','moderate drought','light drought','no drought1','nodrought2','no drought3'])
  )

结果

标签:MODIS,gte,vci,median,2024,归一化,var,ndvi,where
From: https://blog.csdn.net/qq_31988139/article/details/142141912

相关文章

  • 20240918_114105 mysql 认识索引
    关于索引MySQL的索引是数据库管理系统中用于提高数据检索效率的一种数据结构。MySQL支持多种类型的索引,每种索引都有其特定的用途和优化方式。以下是MySQL中常见的几种索引类型:1.主键索引(PrimaryKeyIndex)定义:主键索引是一种特殊的唯一索引,它不允许有NULL值,且表中每一行数据......
  • 海南2024下半年软考准考证打印时间11月4日开始
    根据2024下半年海南软考考务通知的说明,2024下半年海南软考准考证打印相关事项如下:一、2024下半年海南软考准考证打印时间2024年11月4日-11月12日。二、2024下半年海南软考准考证打印入口网址海南2024年下半年软考准考证打印入口:中国计算机技术职业资格网(https://www.ruankao.org.......
  • 2024.9.18训练记录
    订正昨天早上的模拟赛T1还没做,dp写法好像要记录什么的感觉好麻烦T2考试没做出来,其实是挺裸的dp状态记pair<int,int>\(f[i][j][k]\)表示前\(i\)个物品,拉出来\(j\)个\(1\),\(k\)个\(2\)所需要的\({背包数,最后一个背包剩的空间}\)。可以分讨最后这一位是否被拉出......
  • 职场人该如何学习使用AI大模型,都2024年还不会用AI办公的你真的out了!
    【写在开篇:这是一篇针对非技术背景的职场人,学习和使用AI大模型的完全攻略。】非技术背景的职场人想要学习和使用AI大模型,可以遵循以下步骤:1.基础学习:首先,需要掌握人工智能的基础知识,包括但不限于机器学习、深度学习等领域。可以通过阅读《ArtificialIntelligence:AMod......
  • 江西2024下半年软考准考证打印时间11月4日9点后开始
    根据2024下半年江西软考考务通知的说明,2024下半年江西软考准考证打印相关事项如下:一、2024下半年江西软考准考证打印时间2024年11月4日9点-11月12日8点30分。二、2024下半年江西软考准考证打印入口网址江西2024年下半年软考准考证打印入口:中国计算机技术职业资格网(https://www.r......
  • 山东2024下半年软考准考证打印时间11月5日9点后开始
    根据2024下半年山东软考考务通知的说明,2024下半年山东软考准考证打印相关事项如下:一、2024下半年山东软考准考证打印时间2024年11月5日9:00。二、2024下半年山东软考准考证打印入口网址山东2024年下半年软考准考证打印入口:中国计算机技术职业资格网(https://www.ruankao.org.cn/)......
  • 【视频教程】手把手AppWizard轻松制作一个emWin滑动主界面控制框架,任意跳转控制(2024-0
    现在的新版AppWizard已经比较好用,用户可以轻松的创建各种项目常规界面。比如早期创建一个支持滑动的主界面框架,并且可以跳转各种子界面,仅仅界面布局和各种图片格式转换都要花不少时间,而现在使用AppWizard,可以说轻轻松松,毫不费力。用户唯一要做的就是根据自己的芯片性能做一定的速度......
  • 河南2024下半年软考准考证打印时间11月2日开始
    根据2024下半年河南软考考务通知的说明,2024下半年河南软考准考证打印相关事项如下:一、2024下半年河南软考准考证打印时间考前一周。二、2024下半年河南软考准考证打印入口网址河南2024年下半年软考准考证打印入口:中国计算机技术职业资格网(https://www.ruankao.org.cn/)。点击进入......
  • 实操触发器的使用 mysql 20240918_102020
    需求新建日志表用于记录老师表的数据化情况起个名字teacher_log需要的列idoperationmsg建老师日志表CREATETABLEteacher_log( idINTPRIMARYKEYAUTO_INCREMENT, operationVARCHAR(11)NOTNULL, msgVARCHAR(200)NOTNULL);定义添加触发器如果往老师表tea......
  • 2024百元蓝牙耳机哪个好?2024四款优选爆款百元机型推荐!
    随着科技的不断进步和消费者需求的多样化,蓝牙耳机市场在近年来迎来了爆发式的增长,2024年,百元级别的蓝牙耳机以其亲民的价格和日益提升的性能,成为了众多消费者的新宠,那么2024百元蓝牙耳机哪个好?在这个价位段,各大品牌纷纷推出了各具特色的产品,旨在满足不同用户的使用场景和个性化......