数据代码分享丨1990-2020年中国省级行政区1km分辨率NDVI数据

原创 等你关注的 地理学人GeoCompany 2023-03-31 10:00 发表于广东

收录于合集

#植被变化3个

#数据分享6个

#代码分享2个

1990-2020年中国省级行政区1km分辨率NDVI数据

时间分辨率:年

空间分辨率:1 km

数据大小:753 MB

数据时间范围:1990-2020年

Dataset

内容概要

归一化植被指数(NDVI, Normalized Difference Vegetation Index)是监测区域和全球植被的有效指标,可以准确反映地表植被覆盖状况,被广泛应用于目前的植被研究。

本文参考徐新良研究员“中国30m逐年NDVI最大值数据集”的计算方法,利用Google Earth Engine(GEE)平台,使用美国陆地卫星Landsat 5/7/8遥感影像计算得到1990-2020年逐年度的中国省级行政区(34个)1km尺度NDVI最大值数据集。数据数值范围为-10000-10000(由-1-1放大10000倍),格式为GeoTIFF。

经合并预处理后的2015年全国NDVI数据

计算代码

利用GEE平台进行NDVI计算,代码如下:

//**欢迎关注"地理学人GeoCompany"公众号**\\//2012年的NDVI计算调用Landsat7影像库//2013年及之后的NDVI计算调用Landsat8影像库var l5 = ee.ImageCollection("LANDSAT/LT05/C01/T1_TOA");var l7 = ee.ImageCollection("LANDSAT/LE07/C01/T1_SR");var l8 = ee.ImageCollection("LANDSAT/LC08/C01/T1_RT_TOA");
var provinces = ["study_area"];//提供研究区var startYear = 1990;var endYear = 2020;
for (var year = startYear; year <= endYear; year++) {  print("Processing year: " + year);
  for (var index = 0; index < provinces.length; index++) {    var qu = provinces[index];    var AOI = ee.FeatureCollection("users/xxx/" + qu);//提供研究区路径    var startyear = year + "-01-01";    var endyear = year + "-12-31";    var landsat;        if (year < 2012) {      landsat = l5        .filterBounds(AOI)        .filterDate(startyear, endyear);    } else if (year == 2012) {      landsat = l7        .filterBounds(AOI)        .filterDate(startyear, endyear);    } else {      landsat = l8        .filterBounds(AOI)        .filterDate(startyear, endyear);    }        //Landsat5/7的NDVI计算利用 B4和B3波段    //Landsat8 的NDVI计算利用 B5和B4波段    if (year < 2013) {      var get_NDVI = function(image) {        var NDVI = image.normalizedDifference(['B4', 'B3']).rename(['NDVI']);        image = image.addBands(NDVI);        return image.select("NDVI");      };    } else {      var get_NDVI = function(image) {        var NDVI = image.normalizedDifference(['B5', 'B4']).rename(['NDVI']);        image = image.addBands(NDVI);        return image.select("NDVI");      };    }
    var NDVI_Collection = ee.ImageCollection(landsat)          .map(get_NDVI);
    var greenest = NDVI_Collection.qualityMosaic('NDVI').clip(AOI).multiply(10000).toInt();
    Export.image.toDrive({      image: greenest,      description: qu,      folder : 'ndvi' + year,      fileNamePrefix: qu,      scale: 1000,      maxPixels: 1e13,      region: AOI.geometry(),      fileFormat: 'GeoTIFF',    });        Map.centerObject(AOI, 5);    Map.addLayer(greenest,{palette: ['black', 'green', 'red']},'greenest');  }}

获取方法

关注并星标"地理学人GeoCompany"公众号,后台回复“NDVI”即可获取数据。

参考文献:

徐新良. 中国30m逐年NDVI最大值数据集.资源环境科学数据注册与出版系统, 2022. DOI: 10.12078/2022030801.

*仅供学习交流使用,如有错误欢迎批评指正*

地理学人GeoCompany

这里是分享文献、数据、代码、模型与奇思妙想的地理梦工厂,请尽情探索吧!

17篇原创内容

公众号

撰稿丨GeoDUO

审核丨Geoffy, &DoraGeomon

公众号运营不易,欢迎关注、点赞、分享

收录于合集 #植被变化

3个

上一篇OneEarth文献速递|全球山地森林流失加速对生物多样性热点地区构成威胁下一篇数据分享丨2000-2020年中国30米年际最大NDVI数据集

阅读 359

分享收藏

81

【转载学习】数据代码分享丨1990-2020年中国省级行政区1km分辨率NDVI数据相关推荐

  1. 数据代码分享|PYTHON用NLP自然语言处理LSTM神经网络TWITTER推特灾难文本数据、词云可视化...

    全文下载链接:http://tecdat.cn/?p=28877 作者:Yunfan Zhang Twitter是一家美国社交网络及微博客服务的网站,致力于服务公众对话.迄今为止,Twitter的日活 ...

  2. php熊掌号怎么设置json-ld,织梦DEDECMS熊掌号JSON LD结构化数据代码分享

    织梦DEDECMS熊掌号JSON LD结构化数据代码分享 效果: 注意:@id里面的地址必须是自适应网站或者手机站,pc端链接是不支持! 百度官方JSON-LD示例: { "@context ...

  3. HMC5883L读取数据代码分享

    文章目录 前言 一.硬件.软件准备 二.代码内容 1.寄存器地址的宏定义 2.IIC功能实现函数 3.寄存器A,B与模式寄存器配置(初始化) 4.读取xyz寄存器中的数据 三.总结 前言 最近初入门学 ...

  4. 年处理两万亿条环保数据,埃睿迪要做中国最懂环保的大数据公司

    随着国家<"十三五"生态环境保护规划>明确提出打好大气.水.土壤污染防治的三大战役,生态环境保护已经成为国家重点政策.今天,追求良好的环保生态已成为社会各界的共识,国家 ...

  5. 数据代码分享|Python用NLP自然语言处理LSTM神经网络Twitter推特灾难文本数据、词云可视化与SVM,KNN,多层感知器,朴素贝叶斯,随机森林,GBDT对比

    作者:Yunfan Zhang Twitter是一家美国社交网络及微博客服务的网站,致力于服务公众对话.迄今为止,Twitter的日活跃用户达1.86亿.与此同时,Twitter也已成为突发紧急情况时 ...

  6. 墨言教育分享丨大开眼界!跟着中国9大湿地美景学配色

    看风景,一样可以学配色哦!墨言教育提醒喜欢到处走走看看的设计师朋友,中国这9大湿地美景千万别错过,大自然的颜色真是巧夺天工呀,看到这样的色彩也很难不被吸引吧?一起来欣赏大美中国顺便学配色~

  7. php数组 城市代码,PHP数组,包含了中国省级行政区划名称及代码

    return array( 'BEIJING'=>array( 'code'=>110000, 'name'=>'北京市' ), 'TIANJIN'=>array( 'code ...

  8. 第五届“泰迪杯”数据分析技能赛 经验代码分享

    第五届"泰迪杯"数据分析技能赛 经验/代码分享 品牌:"泰迪杯"数据分析技能赛 组织单位:泰迪杯数据分析技能赛组织委员会.广东泰迪智能科技股份有限公司.广东省工 ...

  9. 【数据分享】1901-2021年1km分辨率逐月平均气温栅格数据(全国/分省/免费获取)

    气温数据是我们最常用的气象指标之一,之前我们给大家分享过1950-2022年0.1° x 0.1°精度的逐月平均气温栅格数据和逐年平均气温栅格数据(均可查询之前的文章获悉详情)! 本次我们分享的是精度 ...

最新文章

  1. 解决IE下jquery ajax无法获得最新数据的问题(IE缓存)
  2. python购物记录程序_python ATM购物程序
  3. add结果 bigdecimal_java.math.BigDecimal.add(BigDecimal augend)方法
  4. 《Node.js区块链开发》
  5. 关于DevExpress的心得
  6. httpclient+Jsoup总结
  7. npm ERR! Unexpected end of JSON input while parsing near '...Comment: https://open'
  8. 15K薪资轻松到手,要低调~
  9. 拓端tecdat|使用R语言进行Metroplis-in-Gibbs采样和MCMC运行分析
  10. SAP License:ERP项目升级评标会偶遇以前同事感慨
  11. GoLang爬取花瓣网美女图片
  12. Vivado使用系列:使用自定义BD
  13. 一元三次方程求解(洛谷)c语言
  14. OpenCV基础(7)OpenCV中的颜色空间
  15. mac最方便的文件管理器--Qspace
  16. 1、每天学习一点点之 contenteditable 属性
  17. MySQL数据库常见错误与解决方法总结
  18. 解决:NoClassDefFoundError: javax/activation/DataSource
  19. UltraEdit 15.10 注册码
  20. Hangzhou Invitation Day1

热门文章

  1. A Game of Thrones(68)
  2. 微软 Windows 10 1903 九月更新官方 ISO 镜像
  3. 如何使用手机登录邮箱?
  4. 用python编写一个点餐程序_用python做一个点餐系统,热敏打印机怎么连?
  5. 服务器又宕机了,怎么办?
  6. 原生js+div+css,动态画结构图/流程图
  7. python就业发展前景_如何参与蓬勃发展的Python就业市场
  8. 2009年中国电脑病毒疫情及互联网安全报告(图)
  9. 汽车金融-贷款分期购车业务系统建设
  10. 读文章笔记-自适应和加速GB生成算法