Ocean Color Forum - Not logged in
Forum Ocean Color Home Help Search Login
Previous Next Up Topic SeaDAS / SeaDAS 6.x - General Questions / Stripe lines problems using SWIR atmospheric correction based on seadas5.4 (locked) (1961 hits)
By Yongjiu Date 2011-12-29 10:23
hi, I came across a problem when using the SWIR atmospheric correction, i found the images were not smooth after correction compared with standard NIR correction using seadas6.2 updated 2, Fedora, here is a similar post (http://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?pid=13531;hl=SWIR%20correction#pid13531), in the post, i know it's because of a bad detector in the 1240nm band, and I tried recommended collection5 which seadas5.4 may work to solve the striple line. however, i cannot process using the same command provided script below:
###################################################################

#!/bin/csh -f
alias rm rm
set dir=/media/disk/jy/
cd $dir
clear
foreach name(*L1B_LAC)
set name=$name
set name2=`echo $name|cut -c 1-14`.GEO``
set name3=`echo $name|cut -c 1-14`.L2``
rm $dir/load_seadas
echo "l2gen,7,ifile='$dir/$name',geofile='$dir/$name2',ofile1='$dir/$name3',slat=34,elat=36,slon=138,elon=140,ctl_pt_incr=1,gas_opt=15,resolution=250,aer_opt=-1,aer_wave_short=1240,aer_wave_long=2130,wait=wait" >> $dir/load_seadas
echo "exit"  >> $dir/load_seadas
echo "inputdata="$name""
seadas -em -b load_seadas
end
alias rm rm -i
###################################################################

here coming the messages:
###################################################################
'nknown option: `-
Usage: csh [ -bcdefilmnqstvVxX ] [ argument ... ].
###################################################################

I guess seadas5.4 may lost sth. so it cannot process the datasets using SWIR correction? but switching to seadas6.2 make the stiple line.
how to solve this?

It was greatly appreciated if anyone give some advice

joe
By gwyn Date 2011-12-29 15:43
Joe -

Try the 6.2 version without the "resolution=250" option.
The newer version of l1bgen_modis* does not interpolate over dead detectors, so any dead detectors will cause striping in high-resolution output.
By Yongjiu Date 2011-12-30 01:36
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]

Previous Next Up Topic SeaDAS / SeaDAS 6.x - General Questions / Stripe lines problems using SWIR atmospheric correction based on seadas5.4 (locked) (1961 hits)



Responsible NASA Official: Gene C. Feldman
Curator: OceanColor Webmaster
Authorized by: Gene C. Feldman
Updated: 27 November 2007
Privacy Policy and Important Notices NASA logo