MRT批处理 MODIS影像

任务：

现在要把2019-2022年8天一期的MOD17A2H转换投影并镶嵌，有4景范围 h25/26, v05/06. 

方法：

1.在MRT的界面上设置好参数可以保存成prm文件

1.1导出指定波段：注意打开prm文件看一下，MOD17A2H有GPP,PsnNet,Psn_QC三个波段，如果导出PsnNet波段，那么prm中这个参数应该是

ORIG_SPECTRAL_SUBSET = ( 0 1 0 )

1.2多景镶嵌：注意INPUT_FILENAMES中是四景，并把INPUT_FILENAME（不带S）注释掉。这一条似乎用INPUT_FILENAME也可以，不修改也没影响，还是四景合一

INPUT_FILENAMES = ( D:\NPP\MOD17A2H.A2019001.h25v05.006.2019079210742.hdf D:\NPP\MOD17A2H.A2019001.h25v06.006.2019079210823.hdf D:\NPP\MOD17A2H.A2019001.h26v05.006.2019079210745.hdf D:\NPP\MOD17A2H.A2019001.h26v06.006.2019079210838.hdf )

1.3 WGS1984设置：在edit Projection Parameter中设置下datum为WGS1984
其实以上所有设置都可以在prm文件中手动修改（文末附上prm）

2.在MRT安装目录/bin 中打开cmd，然后敲入

java -jar MRTBatch.jar -d D:\npp -p Convert2WGS1984.prm -o D:\npp2

其中，D:\npp是输入目录，D:\npp2是输出目录，Convert2WGS1984.prm在bin文件夹中，储存相应参数。

3.再敲入

MRTBatch.bat

就可以执行了。

最后输出的就是四景镶嵌好了的，且投影改成WGS84了

#The "INPUT_FILENAMES" field would be commented. If you want to load multiple input files please uncomment the "INPUT_FILENAMES" field and comment the"INPUT_FILENAME" field. 
#Also the "ORIG_SPECTRAL_SUBSET" field needs to be uncommented and changed to "SPECTRAL_SUBSET". The initial "SPECTRAL_SUBSET" field should be deleted.

INPUT_FILENAMES = ( D:\NPP\MOD17A2H.A2019001.h25v05.006.2019079210742.hdf D:\NPP\MOD17A2H.A2019001.h25v06.006.2019079210823.hdf D:\NPP\MOD17A2H.A2019001.h26v05.006.2019079210745.hdf D:\NPP\MOD17A2H.A2019001.h26v06.006.2019079210838.hdf )

#INPUT_FILENAME = D:\TmpMosaic.hdf

#SPECTRAL_SUBSET = ( 1 )
ORIG_SPECTRAL_SUBSET = ( 0 0 1 )

SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG

SPATIAL_SUBSET_UL_CORNER = ( 39.999999996 91.37851024 )
SPATIAL_SUBSET_LR_CORNER = ( 19.999999998 95.775999513 )

OUTPUT_FILENAME = D:\MOD17A2H2019001.tif

RESAMPLING_TYPE = NEAREST_NEIGHBOR

OUTPUT_PROJECTION_TYPE = GEO

OUTPUT_PROJECTION_PARAMETERS = ( 
 0.0 0.0 0.0
 0.0 0.0 0.0
 0.0 0.0 0.0
 0.0 0.0 0.0
 0.0 0.0 0.0 )

DATUM = WGS84