Hi,
Thank you in advance!It's really took me a lot of time to study the Seadas all the way by myself and I really got stuck for a long time but asked for your help.Firstly I want to check is there anything wrong in my data processing?
I had installed the seadas6.4 on my linux system and just wanted to process the the area of Poyang Lake,namely the latitude is 28.3 to 29.8N and the longitude is between 115.7 and 116.8E with the method SWIR/NIR proposed by Wang(2007). Here is my processing procedure.
Firstly I directly downloaded the L0 file from the
http://oceancolor.gsfc.nasa.gov/cgi/browse.pl?sen=am website without any email and chose the Terra satellite with the data of 2002.08.20.Secondly I processed the L0 date from the seadas GUI and didn't change any default.Then I set the L1A_LAC as my input file and got the L1B file and finally I chose the chlor_a,CDOM,all bands of the Rrs and nLw as my output and switched the Aerosol mode option to 2-band model selection with NIR/SWIR,is there anything wrong in my processing?And the follow was my .par file
ifile=/home/junjun/T2002232035000.L1B_LAC
geofile=/home/junjun/T2002232035000.GEO
ofile1=/home/junjun/T2002232035000.L2_LAC
l2prod1=chlor_a,Kd_490,cdom_index,Rrs_412,Rrs_443,Rrs_469,Rrs_488,Rrs_531,Rrs_547,Rrs_555,Rrs_645,Rrs_667,Rrs_678,Rrs_748,Rrs_859,Rrs_869,Rrs_1240,Rrs_1640,Rrs_2130,nLw_412,nLw_443,nLw_469,nLw_488,nLw_531,nLw_547,nLw_555,nLw_645,nLw_667,nLw_678,nLw_748,nLw_859,nLw_869,nLw_1240,nLw_1640,nLw_2130
spixl= 1
epixl= -1
dpixl= 1
sline= 1
eline= -1
dline= 1
ctl_pt_incr= 1
proc_ocean= 1
atmocor= 1
proc_land= 1
proc_sst= 1
resolution= 1000
gas_opt= 15
pol_opt= 3
aer_opt= -9
aermodmin= -1
aermodmax= -1
aermodrat= 0.00000
mumm_alpha= 1.94500
mumm_gamma= 1.00000
mumm_epsilon= 1.00000
aer_rrs_short= -1.00000
aer_rrs_long= -1.00000
aer_swir_short=1240
aer_swir_long=2130
aer_wave_short= 748
aer_wave_long= 869
aer_iter_max= 10
brdf_opt= 7
iop_opt= 0
glint_opt= 1
outband_opt= 2
filter_opt= 1
filter_file=$OCDATAROOT/modisa/msl12_filter.dat
met1=/home/junjun/seadas6.4/run/var/anc/2002/232/N200223200_MET_NCEPN_6h.hdf
met2=/home/junjun/seadas6.4/run/var/anc/2002/232/N200223206_MET_NCEPN_6h.hdf
met3=/home/junjun/seadas6.4/run/var/anc/2002/232/N200223206_MET_NCEPN_6h.hdf
ozone1=/home/junjun/seadas6.4/run/var/anc/2002/231/N200223100_O3_TOMSOMI_24h.hdf
ozone2=/home/junjun/seadas6.4/run/var/anc/2002/232/N200223200_O3_TOMSOMI_24h.hdf
ozone3=/home/junjun/seadas6.4/run/var/anc/2002/232/N200223200_O3_TOMSOMI_24h.hdf
no2file=$OCDATAROOT/common/no2_climatology.hdf
sstfile=/home/junjun/seadas6.4/run/var/anc/2002/232/N2002232_SST_OIV2AVAM_24h.nc
icefile=/home/junjun/seadas6.4/run/var/anc/2002/232/N200223200_SEAICE_NSIDC_24h.hdf
land=$OCDATAROOT/common/landmask.dat
water=$OCDATAROOT/common/watermask.dat
gain=[0.9999,0.9975,1.0033,1.0036,0.9893,0.9901,0.9905,1.0241,0.9909,0.9975,0.9976,1.0161,1.00,1.00,1.00,1.00]
offset=[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,0.0]
albedo= 0.0270000
rhoamin= 0.000200000
qaa_adg_s= 0.0150000
chloc2_wave=[488,547]
chloc2_coef=[0.2500,-2.4752,1.4061,-2.8237, 0.5405]
chloc3_wave=[443,488,547]
chloc3_coef=[0.2424,-2.7423,1.8017, 0.0015,-1.2280]
chloc4_wave=[]
chloc4_coef=[]
chlclark_wave=[443,488,547]
chlclark_coef=[0.789273,-3.925523,11.637764,-27.157997,27.936958,-10.398587]
nlwmin= 0.150000
wsmax= 12.0000
tauamax= 0.300000
epsmin= 0.800000
epsmax= 1.35000
glint= 0.00500000
windspeed= -1000
windangle= -1000
pressure= -1000
ozone= -1000
watervapor= -1000
relhumid= -1000
sunzen= 70.0000
satzen= 60.0000
maskland= 1
maskcloud= 1
maskglint= 0
maskbath= 0
masksunzen= 0
masksatzen= 0
maskhilt= 1
maskstlight= 0
giop_wave=[412,443,469,488,531,547,555]
giop_fit_opt= 1
giop_maxiter= 500
giop_aph_opt= 2
giop_adg_opt= 1
giop_bbp_opt= 3
giop_aph_file=$OCDATAROOT/common/aph_default.txt
giop_adg_file=$OCDATAROOT/common/adg_default.txt
giop_bbp_file=$OCDATAROOT/common/bbp_default.txt
giop_aph_s= 0.500000
giop_adg_s= 0.0145000
giop_bbp_s= 1.03373
First I want you to help me to confirm was this the right way of data processing.Second I got confused when I displayed the L2 product,the district of lake was divided into pieces and how can I change them into a continuous body of water with color image? Any suggestion would be appreciated! Thank you very much!
