first script

1 - NO NEED for field data, nor for atmospheric correction
2 - this is demonstrated in this website, using a variety of hyper/multi spectral data
Requirements are
1 - homogeneous water body and atmosphere
2 - some coverage of optically deep water
3 - some coverage of dry land
Problems are
1 - the precision on estimated depth is found wanting, because the noise-equivalent change in radiance  of accessible data is too high for shallow water column correction work 
2 - radiance data should be preprocessed by the provider at level 1 in order to improve S/N ratio
3 - exponential decay: the deeper/darker the bottom, the poorer the performances
I keep digging
until suitable data
become available
nice -20  4SM.4.08 -Process/06_JUL_2011  -Origin/Tarawa_Kiribati                            \
-DB/tarawa-subset/27_1_0_0/3_1/650_600/718.079_169.221/1_1                            \
-Mis/Tarawa/Kiribati/Spot/XS/NA/UTM_19_008/0.02_0.02/01_JAN_1986                   \
-MakePIX                                                                                                                                     \
-Import/data*tarawa-full/dbnc_3_0_0_0/R1630_L1860/Origin_709.049_182.231_0_0/chIn1,3/chOut1,3/import.shp    \

Modify Channel Descriptor    file="tarawa-subset.pix";

cm01=" 1  1_Raw_550 nm";                    dboc=1  ; r mcd
cm01=" 2  2_Raw_650 nm";                    dboc=2  ; r mcd
cm01=" 3  3_Raw_850 nm";                    dboc=3  ; r mcd
cm01=" 4 Special_Effect_mask; use pctmSE";            dboc=4  ; r mcd
cm01=" 5 wZ: use pctwZ";                    dboc=5  ; r mcd
cm01=" 6 wR: use pctwR";                    dboc=6  ; r mcd
cm01=" 7 Z0->250 in dm: use pctZ250";                dboc=7  ; r mcd
cm01=" 8 Z: use pctZ20";                    dboc=8  ; r mcd
cm01=" 9 B: use pctB";                        dboc=9  ; r mcd
cm01="10 Bottom Typing: use pctClassified";            dboc=10 ; r mcd
cm01="11 LSBCC_blue ";                        dboc=11 ; r mcd
cm01="12 LSBCC_green";                        dboc=12 ; r mcd
cm01="13 LSBCC_red  ";                        dboc=13 ; r mcd
cm01="14 dZED: use pctZ250";                    dboc=14 ; r mcd
cm01="15 BPL_12  ";                        dboc=15 ; r mcd
cm01="16 BPL_13  ";                        dboc=16 ; r mcd
cm01="17 BPL_23  ";                        dboc=17 ; r mcd
cm01="18 void";                            dboc=18 ; r mcd
cm01="19 void";                            dboc=19 ; r mcd
cm01="20 void";                            dboc=20 ; r mcd
cm01="21 void";                            dboc=21 ; r mcd

cm01="22  1_Ldeglinted_550 nm";                    dboc=22 ; r mcd
cm01="23  2_Ldeglinted_650 nm";                    dboc=23 ; r mcd
cm01="24  3_Ldeglinted_850 nm";                    dboc=24 ; r mcd

cm01="25  1_LSB_550 nm";                    dboc=25 ; r mcd
cm01="26  2_LSB_650 nm";                    dboc=26 ; r mcd
cm01="27  3_LSB_850 nm";                    dboc=27 ; r mcd

cm01="28_ZG: Zshallow in centimeters, others at -1";        dboc=28 ; r mcd

WLMin[ 1]=540.0 WL[ 1]=550.0 WLMax[ 1]=560.0   
WLMin[ 2]=640.0 WL[ 2]=650.0 WLMax[ 2]=660.0   
WLMin[ 3]=840.0 WL[ 3]=850.0 WLMax[ 3]=860.0