博文
MRT(MODIS Reprojection Tool)批处理MOD11A1的代码（转载）
||
转自http://blog.sina.com.cn/s/blog_4edba2dd01000bjj.html~type=v5_one&label=rela_nextarticle
MRT(MODIS Reprojection Tool)批处理MOD11A1的代码
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
#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 = ( E:LSTMOD11A1.A2001131.h26v04.004.2003115130411.hdf E:LSTMOD11A1.A2001131.h26v05.004.2003115134041.hdf E:LSTMOD11A1.A2001131.h27v04.004.2003115150352.hdf E:LSTMOD11A1.A2001131.h27v05.004.2003115125635.hdf )
INPUT_FILENAME = E:LSTTmpMosaic.hdf
SPECTRAL_SUBSET = ( 1 )
#ORIG_SPECTRAL_SUBSET = ( 1 0 0 0 0 0 0 0 0 0 0 0 )
SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG
SPATIAL_SUBSET_UL_CORNER = ( 50.0 124.457906149 )
SPATIAL_SUBSET_LR_CORNER = ( 30.0 115.470053838 )
OUTPUT_FILENAME = E:LSTLST.TIF
RESAMPLING_TYPE = NEAREST_NEIGHBOR
OUTPUT_PROJECTION_TYPE = AEA
OUTPUT_PROJECTION_PARAMETERS = (
DATUM = WGS84
OUTPUT_PIXEL_SIZE = 1000
https://wap.sciencenet.cn/blog-331307-277287.html
上一篇：MRT(Modis Reprojection Tool)批处理MODIS影象(以海河流域为例)（转载）
下一篇：Cygwin使用指南
