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'
  });
}