偶看新闻正宗煲仔饭图片大全:MRT(MODIS?Reprojection?Tool)批处理MOD11A1的代码
来源：百度文库 编辑：偶看新闻 时间：2024/05/06 20:47:26
MRT(MODIS Reprojection Tool)批处理MOD11A1的代码
(2007-10-14 10:46:50)
转载 标签： mrt
批处理
modis
以处理MOD11A1数据为例,处理其他MODIS产品类似:代码:day=2001001while [ $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:\LST\MOD11A1.A2001131.h26v04.004.2003115130411.hdf E:\LST\MOD11A1.A2001131.h26v05.004.2003115134041.hdf E:\LST\MOD11A1.A2001131.h27v04.004.2003115150352.hdf E:\LST\MOD11A1.A2001131.h27v05.004.2003115125635.hdf )
INPUT_FILENAME = E:\LST\TmpMosaic.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:\LST\LST.TIF
RESAMPLING_TYPE = NEAREST_NEIGHBOR
OUTPUT_PROJECTION_TYPE = AEA
OUTPUT_PROJECTION_PARAMETERS = (
DATUM = WGS84
OUTPUT_PIXEL_SIZE = 1000