利用MRT进行多年LAI数据（MOD15A2）拼接和投影转换

最近需要处理连续10年的MODIS LAI产品数据（MOD15A2）。根据自己的需求，写了一个MS-DOS批处理小程序，主要解决的问题是：YEAR和DOY两重循环，并将YEAR和DOY连接成YEAR+DOY的形式（如2001001）用以搜索符合条件的待拼接文件。

需要注意的地方（代码中红色标注处）：

1、用set设置数据文件夹环境变量时，等号的前后都不能出现空格。如果等号前面有空格，则环境变量名称变为“MRTDATADIR ”，如果等号后有空格，则数据文件夹路径前面多出空格会导致MRT工具找不到相关数据文件。

2、DOY的设置采用了“曲线救国”方案，即先设置为10001的形式，然后截取最后三个字符作为标准DOY。

3、DOS下的字符截取操作：%DOY：~-3%表示截取DOY变量的最后三个字符，等效于%DOY：~4,3%，表示从DOY变量的第5个字符开始截取3个字符。

4、进行YEAR的循环时需要对DOY进行初始化。

代码如下：

@echo off

rem **This is a MS-DOS batch file to do multi-yearly mosaicing and resampling preprocessing by using MRT

rem **author: pengbin, institute of remote sensing and digital earth, CAS

rem **Memo: developed at May 7th, 2013

rem **Set the MRTDATADIR environmental var to the MRT data directory.

set MRTDATADIR=D:\MOD15A2_h24-26_v4-5

rem **set batch data start time

set /a YEAR = 2001

set /a DOY = 10001

rem **set batch data end time

set /a ENDYEAR = 2010

set /a ENDDOY = 10365

:start

if %YEAR% leq %ENDYEAR% (goto ORDER1) else exit

:ORDER1

if %DOY% leq %ENDDOY% (goto ORDER2) else (

set /a YEAR = %YEAR% +1

set /a DOY = 10001

goto start

)

:ORDER2

set /a DAY = %YEAR%%DOY:~-3%

rem **save the file name into a notepad

dir %MRTDATADIR%\*%DAY%.*.hdf/a/b/s > %MRTDATADIR%\MOSAICINPUT.TXT

rem **execute mosaic

mrtmosaic.exe -i %MRTDATADIR%\MOSAICINPUT.TXT -o %MRTDATADIR%\MOSAIC_TMP.hdf

rem **execute resample

resample.exe -p %MRTDATADIR%\MOD15A2.prm -i %MRTDATADIR%\MOSAIC_TMP.hdf -o %MRTDATADIR%\%DAY%.tif

set /a DOY = %DOY% + 8

goto ORDER1

转载本文请联系原作者获取授权，同时请注明本文来自彭彬科学网博客。
链接地址： http://blog.sciencenet.cn/blog-430991-687486.html

利用MRT进行多年LAI数据（MOD15A2）拼接和投影转换相关推荐

利用MRT进行Modis NDVI数据（MOD13Q1）投影变换格式转换操作图文教程
本实例以Modis NDVI(MOD13Q1,空间分辨率为250m)一景影像数据为例,演示利用MRT进行Modis NDVI影像变换,主要内容包括:将.hdf格式转为.tif格式,将坐标系转为Albe ...
利用MRT对MODIS数据进行批量重投影+批量波段合成
写在前面: 官方渠道已经下载不到MRT了,为什么呢?退休了呗. LP DAAC - The downloadable MODIS Reprojection Tool (MRT) and MRTSwat ...
Windows10下VB6.0开发——利用PictureBox控件实现数据点实时绘图
前言:VB中可以利用PictureBox控件实现数据点实时绘图功能. 1. 控件PictureBox 下面是PictureBox的图标: 将该控件拖到目标位置后调整它的形状大小,修改它的背景色 ...
第二、三代基因组测序数据混合拼接软件综述
第二.三代基因组测序数据混合拼接软件综述王昊1,陈挺2* (1.清华大学计算机科学与技术系,北京 10084: 2.清华大学计算机系.清华大学人工智能研究院,北京 10084) 摘要 ...
利用python进入数据分析之数据规整化：清理、转换、合并、重塑（二）
数据转换移除重复数据 In [106]: data = DataFrame({'k1': ['one'] * 3 + ['two'] * 4,'k2': [1, 1, 2, 3, 3, 4, 4]} ...
利用Python进行数据分析：数据规整（基于DataFrame）
利用Python进行数据分析:数据规整在许多应用中,数据可能分散在许多文件或数据库中,存储的形式也不利于分析.本部分关注可以聚合.合并.重塑数据的方法. 文章目录利用Python进行数据分析:数据 ...
最新数据处理之最新DSM（AW3D30）数据批量拼接--文末附数据获取方式
本文主要针对JAXA发布的最新版本3.1&3.2(AW3D30)数据进行数据拼接处理. JAXA利用ALOS的PRISM立体模式制作了全球的DSM数据,并免费分布了30米版本产品.详细的数据介 ...
长时间序列遥感数据植被物候提取/遥感数据产品分析暨MODIS NDVILAI多年产品数据批处理分析/Python长时间序列遥感数据处理及在全球变化、物候提取、植被变绿与固碳分析、生物量估算与趋势分析
基于MATLAB长时间序列遥感数据植被物候提取与分析 1.本课程基于matlab语言 2.提供所有代码 3.以实践案例为课程内容主线,原理与操作相结合 4.根据讲解内容,布置作业,巩固所学内容及拓展在 ...
利用SQl对数据库实行数据拆分与组合
利用SQl对数据库实行数据拆分与组合实现提供以下几种方案: 方法一: WITH CTE AS ( SELECT A.Id,A.[Uid],UserName FROM ( SELECT A.[id], ...

利用MRT进行多年LAI数据（MOD15A2）拼接和投影转换

利用MRT进行多年LAI数据（MOD15A2）拼接和投影转换相关推荐

最新文章

热门文章