SPATIAL_SUBSET_LR_CORNER = ( 29.0 116.0 )两行 如果需要不同的投影,如经纬度投影,可以修改 OUTPUT_PROJECTION_TYPE = UTM 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 ) DATUM = WGS84 UTM_ZONE = 50
如果拼接时只选择了1个波段,可以修改 SPECTRAL_SUBSET = ( 1 ) 当然这句
INPUT_FILENAME = G:\\TmpMosaic.hdf批处理是必须改掉的了 关于像素分辨率,一般留空,也就是默认不改变影像的像素大小。 2关于DOS批处理裁切MOD产品 单次处理
Dir g:\\mod11*A%1*.hdf /s /b >> g:\\n%1.txt
f:\\mrt\\modis\\bin\\mrtmosaic -i \1 0 0 0 0 \
f:\\mrt\\modis\\bin\\resample -p \del g:\\prm%1.prm /q /f del g:\\n%1.txt /q /f del g:\\temp_%1.hdf /q /f 批量处理
for /f %%i in (e:\\lis.txt) do @单次处理文件名 %%i
如果双击运行,用%%i号表示变量。在命令行中%i表示变量。 其中e:\\lis.txt是年份及天数 2008001
2008009 2008017 2008025 2008033 2008041 2008049
G:\\n%1.txt在运行时为g:\\n2008001.txt,其中包括2008年第1天的多景MODIS产品影像名。 Dir g:\\mod11*A%1*.hdf /s /b >> g:\\n%1.txt
用来列出g盘(数据存储盘)以mod11打头的指定日期的的hdf格式文件,即MODIS影像。 这里要说明一下,MODIS产品命名的规则
MYD11A2.A2002185.h27v05.005.2007221162948.hdf
从左至右,前7位表示产品类型MYD11指下午星(MYD)第11号产品(land surface temperature),A2表示处理级别。A2002185表示过境时间是2002年第185天,h27v05是MODIS产品在全球的轨道行列号,中国一般是h27v05、h27v06、h28v05、h28v06,005表示HDF5文件格式,早期MODIS采用HDF4格式存储,所以有时也可见到004的产品。2007等一串数字表示影像处理的年月日时分秒,hdf表示存储格式为EOS-HDF。
了解了MODIS产品命名规则之后,可以用DOS dir命令来挑选符合条件的MODIS产品路径了。 >> g:\\n%1.txt 表示把结果转存为g:\\n%1.txt %1表示接收的第一个参数,即2008001 3Envi生成时间列表 2008001 2008009 2008017 2008025
这样的列表文件可以用在Envi中定义一个pro来生成: pro daylist,startdate,enddate,result_txt_path
;example daylist,2008001,2008013,'example.txt' ; ;
a=(startdate mod 1000) /8 *8+(startdate /1000)*1000+1
b=((enddate mod 1000) /8 +1)*8+(enddate /1000)*1000+1 help,a,b step=8
openw,lun,result_txt_path,/get_lun res=a;
for i=a,b,step do begin
printf,lun,string(res,format='(%\
res=res+step endfor free_lun,lun print,'done'; close,/all end
在命令行窗口敲入
daylist,2009001,2009033,'f:\\g4.txt' 打开f:\\g4.txt如下: 2009001 2009009 2009017 2009025 2009033 2009041
4Envi平台下整合前面工作 留给读者思考 提示:spawn函数