thanks so much, the image was greatly improved, however there still exist a few stripings. Do i miss sth.?
here is the script i do not use the 'resolution=250', I try to switch the SWIR and NIR according to the turbidity index which was developed by Wang and shi, 2007. also attach two images with and without resolution=250.
#!/bin/csh -f
###################################################################
# ----SeaDAS command mode bat file---
#
# MODIS L1A & GEO data ---> processing---> L1B
#?
###################################################################
alias rm rm
set dir=/home/joe/joe/1/
cd $dir
clear
echo "This bat.file process MODIS L1A & GEO to MODIS L1B"
foreach name(*.L1B_LAC)
set name=$name
set name2=`echo $name|cut -c 1-15`.GEO``
set name3=`echo $name|cut -c 1-15`.stdproc.L2``
rm $dir/load_seadas
##################################################################
echo "l2gen,ifile='$dir/$name',geofile='$dir/$name2',ofile1='$dir/$name3',slat=34,elat=36,slon=138,elon=140,ctl_pt_incr=1,gas_opt=15,aer_opt=-9,aer_wave_short=1240,aer_wave_long=2130,wait=wait" >> $dir/load_seadas
#echo "l2gen,ifile='$dir/$name',geofile='$dir/$name2',ofile1='$dir/$name3',slat=32.25,elat=34.25,slon=130.75,elon=132.75,ctl_pt_incr=1,resolution=-1,/wait" >> $dir/load_seadas
echo "exit" >> $dir/load_seadas
###################################################################
#
#
echo "inputdata="$name""
seadas -em -b load_seadas
end
alias rm rm -i
[img][/img]