GEE开发之Modis_NPP数据分析和获取
前言:主要介绍使用GEE获取Modis_NPP数据。
1.NPP
- 净初级生产量(NPP, Net Primary Production)指在初级生产过程中,植物光合作用固定的能量中扣除植物呼吸作用消耗掉的那部分,剩下的可用于植物的生长和生殖的能量。生产量通常用每年每平方米所生产的有机物质干重[g/(m2·a)]或每年每平方米所固定能量[J/(m2·a)]表示。
- MOD17A3H V6 产品以 5 亿像素分辨率提供有关年度净初级生产力 (NPP) 的信息。年度 NPP 源自给定年份的 45 个 8 天净光合作用 (PSN) 产物 (MOD17A2H) 的总和。PSN 值是 GPP 与维持呼吸 (MR) (GPP-MR) 之差。
- MOD17A3HGF V6 产品以 5 亿像素分辨率提供有关年度净初级生产力 (NPP) 的信息。年度 NPP 源自给定年份的所有 8 天净光合作用 (PSN) 产品 (MOD17A2H) 的总和。PSN 值是总初级生产力 (GPP) 和维持呼吸 (MR) (GPP-MR) 之差。
- 官网链接https://developers.google.cn/earth-engine/datasets/catalog/MODIS_006_MOD17A3HGF
- https://developers.google.com/earth-engine/datasets/catalog/MODIS_061_MOD17A3HGF
2.MOD系列和MYD系列区别
- Terra卫星的观测数据为MOD数据,Terra是于当地的上午10:30过境。
- aqua卫星观测数据为MYD数据,auqa是于当地下午1:30过境。
- A3是NPP数据,每一年进行合成。
- 006_A3系列只有NPP数据,061_A3系列NPP和GPP都有
3.MOD17A3H_v006(500m/1年)
var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset = ee.ImageCollection('MODIS/006/MOD17A3H').filter(ee.Filter.date('2014-01-01', '2014-12-31'));
print(dataset)
var npp = dataset.select('Npp');
var nppVis = {
min: 0.0,
max: 19000.0,
palette: ['bbe029', '0a9501', '074b03'],
};
Map.addLayer(npp.mean().clip(geometry), nppVis, 'NPP');
4.MOD17A3HGF_v006(500m/1年)
var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset = ee.ImageCollection('MODIS/006/MOD17A3HGF').filter(ee.Filter.date('2014-01-01', '2014-12-31'));
print(dataset)
var npp = dataset.select('Npp');
var nppVis = {
min: 0.0,
max: 19000.0,
palette: ['bbe029', '0a9501', '074b03'],
};
Map.addLayer(npp.mean().clip(geometry), nppVis, 'NPP');
5.MOD17A3HGF_v061(500m/1年)
- 注意:指数包括NPP、GPP、NPP_QC
var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset = ee.ImageCollection('MODIS/061/MOD17A3HGF').filter(ee.Filter.date('2014-01-01', '2014-12-31'));
print(dataset)
var npp = dataset.select('Npp');
var nppVis = {
min: 0.0,
max: 19000.0,
palette: ['bbe029', '0a9501', '074b03'],
};
Map.addLayer(npp.mean().clip(geometry), nppVis, 'NPP');
6.MYD17A3H_v006(500m/1年)
6.1 MYD17A3H下的指数
var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset = ee.ImageCollection('MODIS/006/MYD17A3H').filter(ee.Filter.date('2014-01-01', '2014-12-31'));
print(dataset)
6.2 遥感影像查看
var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset = ee.ImageCollection('MODIS/006/MYD17A3H').filter(ee.Filter.date('2014-01-01', '2014-12-31'));
print(dataset)
var npp = dataset.select('Npp');
var nppVis = {
min: 0.0,
max: 19000.0,
palette: ['bbe029', '0a9501', '074b03'],
};
Map.addLayer(npp.mean().clip(geometry), nppVis, 'NPP');
7.MYD17A3HGF_v061(500m/1年)
- 注意:有NPP、GPP、NPP_QC指数
var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset = ee.ImageCollection('MODIS/061/MYD17A3HGF').filter(ee.Filter.date('2014-01-01', '2014-12-31'));
print(dataset)
var npp = dataset.select('Npp');
var nppVis = {
min: 0.0,
max: 19000.0,
palette: ['bbe029', '0a9501', '074b03'],
};
Map.addLayer(npp.mean().clip(geometry), nppVis, 'NPP');
8.MYD17A3HGF_v006(500m/1年)
var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset = ee.ImageCollection('MODIS/006/MYD17A3HGF').filter(ee.Filter.date('2014-01-01', '2014-12-31'));
print(dataset)
var npp = dataset.select('Npp');
var nppVis = {
min: 0.0,
max: 19000.0,
palette: ['bbe029', '0a9501', '074b03'],
};
Map.addLayer(npp.mean().clip(geometry), nppVis, 'NPP');
9.NPP年数据下载(以MOD17A3HGF为例)
var geometry = ee.FeatureCollection('users/www1573979951/HHLYt');
var dataset = ee.ImageCollection("MODIS/061/MOD17A3HGF");
for(var i=2021;i<=2023;i++){
var data_collection = dataset.filterDate(i+'-01-01',i+'-12-31').select('Npp');
print(data_collection)
var YR_collection = data_collection.mean().multiply(0.0001).clip(geometry);//单位换算
Export.image.toDrive({
image: YR_collection,
description: "HHLY_NPP_"+i,
fileNamePrefix: "HHLY_NPP_"+i,
scale: 1000,
crs: "EPSG:4326",
region: geometry,
maxPixels: 1e13,
folder: 'Modis'
});
}