以处理MOD11A1数据为例,处理其他MODIS产品类似:
代码:
day=2001001
while [ $day -le 2001365 ] do for trl in MOD11A1.A$day.*.hdf do if [ -f $trl ] then HDFFILES=$(ls MOD11A1.A$day.*.hdf) echo $HDFFILES > MOSAICINPUT.TXT MRTMOSAIC -i MOSAICINPUT.TXT -s "1 0 0 0 0 0 0 0 0 0 0" -o MOSAIC_TMP_$day.hdf RESAMPLE -p LST.PRM -i MOSAIC_TMP_$day.hdf -o MOD11A1_$day.tif rm MOSAIC_TMP_$day.hdf rm MOD11A1.A$day.*.hdf let "day=day+1" break else let "day=day+1" break fi done done LST.PRM参数文件:
#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.
INPUT_FILENAME = E:\LST\TmpMosaic.hdf SPECTRAL_SUBSET = ( 1 ) SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG SPATIAL_SUBSET_UL_CORNER = ( 50.0 124.457906149 ) OUTPUT_FILENAME = E:\LST\LST.TIF RESAMPLING_TYPE = NEAREST_NEIGHBOR OUTPUT_PROJECTION_TYPE = AEA OUTPUT_PROJECTION_PARAMETERS = ( DATUM = WGS84 OUTPUT_PIXEL_SIZE = 1000 |
|