|||
此处的产品为MOD02HKM,我只是从官网上下了一部分影像,在批处理时,没有顺序,所以就把每个文件名中的天数全部自己写出来,因为影像相对较少,所以可行,当文件较多时,可以再修改文件的读取方式
(另外:L1B产品的处理还需要相应的Geo文件,即MOD03)
(直接复制粘贴的代码$都变了,贴张图说明一下)
day=(254 255 256 257 258 258 259 260 260)
time=(130 215 120 200 240 245 325 325 230)
i=0
while [ $i -lt 9 ]
do
for trl in MOD02HKM.A2016${day[i]}.0${time[i]}.006.*.hdf
do
for geo in MOD03.A2016${day[i]}.0${time[i]}.006.*.hdf
do
if [ -f $trl ]
then
swath2grid -if=$trl -gf=$geo -of=MOD02HKM.A2016${day[i]}.0${time[i]} -off=GEOTIFF_FMT -sds="EV_500_RefSB, 1, 1; EV_250_Aggr500_RefSB, 1" -kk=NN -oproj=GEO -oprm=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 -osp=8
fi
done
let "i=i+1"
break
done
done
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 04:54
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社