subset_stitch_grid.exe

1.TIF 转换

使用方法为File，然后选择第二项，导入数据处理即可

重采样类型：Resampling Type：

• Bilinear：双线性
• Nearest Neighbor 就近
• Cubic Convolution 三次卷积

点击 batch run 开始执行

批量处理时不能使用默认的区域参数，手动设置

2.拼接

2.1 工具拼接

tool->switch/subset 切换工具

2.2 python 脚本拼接

保存后的是一个txt文件，Unix编码（这一点很重要）的txt文件，可以用记事本打开这里记录了处理参数，然后咱们需要做的就是批量更改参数并调用HEG来处理。需要修改的有这几个参量

NUMBER_INPUTFILES 文件个数
INPUT_FILENAMES 文件名（多个文件名之间以‘|’分隔）
OUTPUT_FILENAME 输出文件名
OUTPUT_STITCHED_FILENAME 这一行属于过程文件的输出路径，虽然选了NO，但是批量调用的时候它还是会给你输出，建议这一行删掉（不影响程序运行）

NUM_RUNS = 1

BEGIN
NUMBER_INPUTFILES = 1
INPUT_FILENAMES = E:\MOD13Q1.A2021305-China\MOD13Q1.A2021305.h21v03.061.2021322012318.hdf
OBJECT_NAME = MODIS_Grid_16DAY_250m_500m_VI|
FIELD_NAME = 250m 16 days NDVI|
BAND_NUMBER = 1
SPATIAL_SUBSET_UL_CORNER = ( 59.999999995 46.671714797 )
SPATIAL_SUBSET_LR_CORNER = ( 49.999999996 79.99999998 )
OUTPUT_OBJECT_NAME = MODIS_Grid_16DAY_250m_500m_VI|
OUTGRID_X_PIXELSIZE = 0.0020833366514795377
OUTGRID_Y_PIXELSIZE = 0.0020833366514795377
RESAMPLING_TYPE = BI
OUTPUT_PROJECTION_TYPE = GEO
ELLIPSOID_CODE = WGS84
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  )
OUTPUT_FILENAME = E:\MOD13Q1.A2021305-China\MOD13Q1.A2021305.h21v03.061.2021322012318_MODIS_Grid_16DAY_250m_500m_VI.tif
SAVE_STITCHED_FILE = NO
OUTPUT_STITCHED_FILENAME = E:\MOD13Q1.A2021305-China\MOD13Q1.A2021305.h21v03.061.2021322012318_MODIS_Grid_16DAY_250m_500m_VI_stitched_.hdf
OUTPUT_TYPE = GEO
END