MODIS数据批量裁剪并合成月尺度数据：以MOD13A1为例

最新推荐文章于 2023-04-22 21:36:33 发布
最新推荐文章于 2023-04-22 21:36:33 发布
阅读量5.9k 收藏 81
点赞数 11
分类专栏： RS&GIS 文章标签： 经验分享
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45024495/article/details/121309893
版权

1 excel表里列举全部【栅格计算器】 指令。

 其中MODIS文件名可利用python代码全部输出

代码（注意更改路径）：

import os
import glob
path = "E:/2WD/FVC/1MRT/2002/*.tif"
d=glob.glob(path)
for i in d:
    print(i)

 结果：

 2 批量剪裁

  ArcMap中打开矢量文件及某年份全部MODIS数据

【spatial analyst tool】-【extract】-【extract by mask】-【batch】

 输入指令

3 批量最大值法合成月尺度的NDVI数据

【spatial analyst tool】-【local】-【cell statistics】-【batch】

 输入指令：

 4 结果

丑唐的修炼之旅
关注
  • 11
    点赞
  • 踩
  • 81
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 3
    评论
专栏目录
【ArcGIS遇上Python】Python实现Modis NDVI批量化月最大合成
「 刘一哥与GIS的故事」
08-29 1万+
最大合成法（MVC）可以在Envi中的Band Math中进行，式子是B1>B2,但是无法批量化；本文实现在ArcGIS中利用Python代码批量进行，如下： 用到的Modis NDVI数据是在MRT中进行拼接与转投影后的月数据，一个月有两期，数据格式是.tif，文件名的格式如：20040101.1_km_16_days_NDVI.tif，20040102.1_km_16_days_NDVI.tif代表2004年1月份的两期数据。本次处理的为2004-2013年的10年的数据。 Python代码如
GEE下载MODIS逐月NDVI数据并合成最大值
2301_77925375的博客
06-05 1526
首先，通过定义一个矩形范围来确定感兴趣的地理区域。然后，导入MODIS 250m NDVI数据集，并只选择其中的NDVI波段。函数从整个数据集中提取该月份的数据，然后从所有数据中选出最大值。这种方法可以确保我们得到的是该月份的NDVI最大值。需要注意的是，在实际使用时，需要替换代码中缺失的xmin、ymin、xmax和ymax参数，以确定所需区域的精确范围。函数来过滤出所需时间范围内的数据。函数中对每个月份的数据进行处理。函数在1到12月之间生成一个数值序列，并在。对其进行裁剪，并使用。
MODIS16天数据月合成/季合成
ronvicki的博客
11-11 8976
1.MRT重投影、转换格式 重投影后命名为MOD13Q1.A2019129.250m_16_days_NDVI.tif 2.重命名 5月单独一期，即序号为129的一期。 首先去掉前面的MOD13Q1. 再将001.250替换为0101 将129.250替换为0501 得到A20200501_16_days_NDVI.tif 3.按月合成 ...
GEE下载MODIS MOD13Q1 16天合成NDVI数据
ENVI_GIS_xiaobai的博客
12-23 1555
GEE、MOD13Q1
Python GDAL MODIS ET（MOD16A2GF）8天合成月尺度
gis风轻云淡的专栏
02-23 6462
最近处理一批长时间序列MODIS蒸散发数据，官方数据为8-days，想转成月尺度的数据，想用Python+GDAL批处理，本来觉得应该也不算多复杂的一个工作，想着偷懒在网上搜一下看看有没有现成的代码可以参考，结果中英文网页查完后居然都没有找到合适的代码，感觉不太能理解，现有的几篇博客介绍的处理方法跟我自己设想的差距有点远，都是直接找每个月的大概4景，然后做个平均值或者最大值合成，我个人感觉这样处理不是特别科学，因此零零碎碎花了大概2~3天时间自己写了一个代码，我自己写代码水平一般，很多地方可...
MODIS NDVI下载处理 MOD13A1， win10
balabalawho的专栏
12-23 3283
下载链接： LP DAAC - Search 以前是用MODIS Reprojection Tool (MRT) 进行处理，但是根据USGS网站， 以上工具已经退休： The downloadable MODIS Reprojection Tool (MRT) and MRTSwath tool have been retired. Users are encouraged to use NASA'sHDF-EOS to GeoTIFF (HEG)downloadable conversion...
MODIS数据下载、拼接、裁剪、8天合成月方法以及处理过程中遇到的问题及解决办法（简要整理）
热门推荐
m0_64793686的博客
03-30 1万+
一、MODIS数据下载 ①下载网址：https://ladsweb.modaps.eosdis.nasa.gov/​​​​​​ 建议使用Chrome浏览器，并且需要进行注册账号才可以下载。 ②下载步骤：在官网找到search/Find data，依次选择数据类型、时间、区域，最后进行下载。 ...
月最大值合成法、GEE绘图参数、MODIS数据科普
m0_58127906的博客
04-22 1051
月最大值合成法、GEE绘图参数、MODIS数据科普
【ArcGIS自定义脚本工具】MOD13Q1数据集处理：批量将16-day分辨率的NDVI栅格转换为monthly分辨率
Salierib的博客
02-09 3086
文章目录一、功能介绍二、脚本代码三、工具参数四、工具界面 系列文章目录：ArcGIS自定义脚本编程 一、功能介绍 利用Arcmap中的镶嵌至新栅格工具，镶嵌运算符选择MAXIMUM，将由MOD13Q1数据集预处理得到的16-day、250m分辨率的NDVI栅格批量转为monthly逐月分辨率。 进行镶嵌至新栅格时，平闰年1 ~ 12月份NDVI的输入栅格可参考下图进行选择： 二、脚本代码 #!/usr/bin/python # -*- coding: UTF-8 -*- #!/usr/bin/pyt
modis数据批量处理详细过程xzx.docx
07-23
从MODIS植被指数的下载到批量处理为最大合成法NDVI/EVI数据的全套过程
pyhton批量处理modis数据最大值合成，转投影，重采样转asc
06-22
在modis网站下载全球数据，为日值数据，每一天的数据大概为500M到1G不等，一年数据365天，大概300多G，数据量算比较大的了。原始数据原来分辨率0.05度，现在需要按照半月合成，即得到1年24期数据，结果转投影成1984，重采样成0.5度。 考虑用python进行批量操作，直接得到最终asc结果，可以在matlab进行批量操作。除了最大值合成，还可以平均值合成，代码并且考虑了闰年语非闰年的差异结果举例，如2008年1月最大值合成文件名：mod09cmg200801a一月上旬;mod09cmg200801b一月下旬。整个程序处理过程中显示每处理完成一个月所用的时间。 对基础数据的处理很方便，不需要打开ArcGIS,安装了GIS，自动有python软件，只需要在开始菜单——ARCGIS——PYHON2.7下打开python，菜单下open打开hdf2asc_15day_max_prj_resample.py文件，修改年份，文件夹路径以及目标投影参考影像路径即可运行。
MODIS MOD13A1.zip
03-16
以MODIS 13 A1为例，对MODIS影像进行批量提取、镶嵌、重投影，同时进行质量控制。
subset_via_shp_1.rar_MODIS裁剪_modis批量裁剪_矢量数据_裁剪
09-23
在envi+idl中用矢量数据批量裁剪modis数据
GEE google earth engine批量裁剪下载MODIS数据
01-16
请将代码复制到GEE运行，该代码是读取MODIS数据集的NDVI数据，裁剪、下载，可批量处理，使用简单方便，完美运行。
2000年-2023年安徽省积雪覆盖Modis数据（未处理栅格数据）【MOD10A1/A2】
04-04
2000-02-24到2023-02-10的安徽省积雪覆盖数据（未处理栅格数据）【MOD10A1/A2】 里面为分好类的Modis数据
GEE东拼西凑之MOD16A2_ET 月合成
qq_50291870的博客
03-22 1563
GEE之MOD16A2_ET 月合成尝试
GLASS数据预处理/MRT工具批量处理MODIS数据
qq_45689568的博客
06-05 418
卫星数据批处理
【ArcGIS自定义脚本工具】绘图：MODIS数据编号对应月份图
Salierib的博客
10-01 669
一、功能介绍 该脚本以MODIS数据编号为y轴，以编号对应的月份为x轴，绘制散点图。利用这个脚本可以辅助分析MODIS数据的月变化。 二、脚本代码 运行该脚本需要在Arcgis对应的python目录（默认为C:\Python27\ArcGIS10.2\Lib\site-packages）下安装第三方库Matplotlib，本文使用的版本为matplotlib-2.1.0。 import sys ...
如何用R语言批量下载MODIS13Q1数据
最新发布
04-03
使用R语言批量下载MODIS13Q1数据可以通过以下步骤实现： 1. 安装和加载必要的R包：首先，确保你已经安装了`MODIS`和`raster`这两个R包。如果没有安装，可以使用以下命令进行安装： ```R install.packages("MODIS") install.packages("raster") ``` 然后，加载这两个包： ```R library(MODIS) library(raster) ``` 2. 设置MODIS下载参数：使用`MODISoptions()`函数来设置下载参数，包括数据产品、日期范围、下载路径等。例如，设置下载路径为`"C:/MODIS"`，数据产品为`"MOD13Q1"`，日期范围为从`"2019-01-01"`到`"2019-12-31"`： ```R MODISoptions(localArcPath = "C:/MODIS", outDirPath = "C:/MODIS", product = "MOD13Q1", begin = "2019-01-01", end = "2019-12-31") ``` 3. 批量下载MODIS数据：使用`runGdal() `函数来批量下载MODIS数据。例如，下载NDVI数据： ```R runGdal(product = "MOD13Q1", tileH = 18, tileV = c(4, 5), bands = "250m_16_days_EVI", job = "MOD13Q1_download") ``` 其中，`tileH`和`tileV`参数指定了要下载的地理瓦片，`bands`参数指定了要下载的数据集。 4. 查看下载结果：下载完成后，可以使用`listMosaicFiles()`函数来查看下载的文件列表： ```R listMosaicFiles(product = "MOD13Q1", tileH = 18, tileV = c(4, 5), bands = "250m_16_days_EVI") ```

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

丑唐的修炼之旅

创作不易、你的鼓励、我的动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值